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This  dissertation  focuses  on  two  main  aspects.  One,  developing  new  tools  for  the 
robustness  analysis  of  uncertain  linear  systems.  Two,  the  synthesis  of  robust  controllers 
for  uncertain  delay  systems. 

Traditional  methods  to  analyze  the  robust  stability  of  linear  systems  have  depended 
on  structured  singular  value  bounds.  However,  this  approach  neglects  phase  dependence 
of  the  uncertainty  and  results  in  undue  conservatism.  In  this  dissertation,  the  concept  of 
the  Nyquist  robust-stability  margin,  is  introduced  for  characterizing  the  closed-loop 
stability  of  uncertain  systems.  The  approach  makes  direct  use  of  Nyquist  domain 
arguments  and  is  based  on  the  analysis  of  the  perturbed  eigenvalue  loci,  hence  avoiding 
undue  conservatism  through  the  use  of  singular-value  upper  bounds.  A  key  element  in 
the  new  approach  is  Critical  Direction  Theory  applied  to  uncertainties  in  the  Nyquist 
plane.  The  critical  direction  method  is  based  on  recognizing  that,  at  any  given  frequency 
on  the  Nyquist  plane,  there  is  only  one  direction  of  perturbation  of  relevance  to  the 
stability  analysis.    This  allows  the  characterization  of  robust  stability  margins  for 
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uncertain  systems  characterized  by  irregular  perturbation  templates,  a  problem  that  poses 
significant  challenges  to  other  analysis  methods.  Examples  of  practical  relevance  are 
given  to  illustrate  the  application  of  the  new  theory.  Using  the  new  approach,  the 
problem  of  assessing  robust  stability  and  computing  stability  margins  for  SISO  systems 
with  affine  complex  parametric  uncertainties  is  tackled  successfully.  Exact  analytical 
results  are  derived  for  geometrically  simple  uncertainty  sets  such  as  ellipses  and 
rectangles. 

In  the  later  part  of  the  dissertation,  the  synthesis  of  robust  controllers  for  state- 
delayed  and  input-delayed  systems  is  considered.  In  particular.  Sliding  Mode  Control  is 
chosen  as  the  technique  of  choice,  as  it  possesses  the  combination  of  robustness  and 
performance  guarantees  that  one  seeks  in  a  control  system.  Robust  stability  to  the 
chosen  perturbation  characterization  is  rigorously  proven.  Practical  difficulties  in 
implementation  are  pointed  out  and  ways  to  overcome  these  hurdles  are  presented.  Some 
open  questions  in  the  literature  are  brought  out  and  theoretical  analysis  and  analytical 
answers  are  presented.  Finally,  ideas  and  interesting  formulations  are  presented  for 
future  work. 


IX 


CHAPTER  1 
INTRODUCTION 

Modeling  of  any  practical  system  for  control  design  invariably  involves  uncertainty. 
Since  the  exact  model  of  a  process  may  be  difficult  or  impossible  to  determine,  the 
logical  approach  is  to  design  a  control  strategy  based  on  a  suitable  nominal  (i.e.,  most 
likely)  model.  Once  a  system  has  been  identified  and  a  nominal  model  and  associated 
uncertainty  description  is  available,  there  are  two  main  tasks.  The  first  task  is  the 
analysis  of  the  uncertain  system.  The  following  issues  arise  naturally:  The  first,  is  the 
question  of  robustness  analysis.  What  is  the  measure  of  the  maximum  uncertainty  that 
the  system  can  tolerate  and  still  sustain  stability  and  performance?  Alternatively,  given 
the  size  of  the  uncertainty,  one  might  wish  to  study  the  effect  of  the  uncertainty  on  the 
stability  and  performance  of  the  system.  A  system  is  said  to  be  robustly  stable  if  it  can 
sustain  stability  in  spite  of  suffering  from  uncertainty. 

The  second  major  issue  is  the  synthesis  of  robust  controller  to  ensure  stability  and 
performance  of  the  physical  system.  Over  the  past  two  decades,  a  plethora  of  robust 
controller  methods  have  been  investigated  and  implemented  in  a  wide  array  of  industries. 
Although,  robust  analysis  and  synthesis  is  an  intellectually  stimulating  and 
mathematically  challenging  task,  there  are  some  valid  criticisms  that  can  be  leveled 
against  this  approach.  Is  it  a  good  idea  to  design  a  system  for  worst  case  uncertainty?  If 
the  probability  of  worst  case  occurring  is  very  low,  why  should  we  sacrifice  the 
performance  of  the  system  by  designing  for  the  worst  case?  Should  the  uncertainties  be 
simply  bounded  by  norms  (or)  should  there  be  probabilistic  ideas  inherently  associated 
with  the  uncertainty  descriptions?  These  are  areas  of  intense  interest  and  importance  and 
the  interested  reader  is  referred  to  the  research  monograph  by  Smith  and  Dahleh  [1994] 


and  the  references  therein. 

In  this  dissertation  attention  will  be  focused  on  both  the  analysis  and  the  synthesis 
aspects  of  robust  control.  The  thesis  is  organized  as  follows:  Chapters  2-3  address  the 
robustness  analysis  problem,  while  Chapters  4-5  deal  with  controller  design  for  robustly 
stabilizing  a  class  of  uncertain  delay  systems.  In  Chapter  6,  future  research  directions  are 
presented. 

In  Chapter  2,  the  concept  of  the  Nyquist  robust-stability  margin,  for  characterizing 
the  closed-loop  stability  of  uncertain  systems  is  introduced.  The  predominant  robustness 
analysis  methods  typically  rely  on  arguments  based  on  real  singular  value  or  magnitude 
bounds.    Such  is  the  case  of  well  known  multivariable  stability  margin  A:,„  [Safonov, 

1982]  and  the  structured  singular  value  ji  [Doyle,  1982].  In  contrast,  the  proposed  new 
approach  makes  direct  use  of  Nyquist  arguments  and  is  based  on  the  analysis  of  the 
perturbed  eigenvalue  loci,  hence  avoiding  undue  conservatism  that  may  result  through  the 
use  of  singular-value  upper  bounds.  A  key  element  in  the  new  approach  is  Critical 
Direction  Theory  applied  to  uncertainties  in  the  Nyquist  plane.  The  critical  direction 
method  is  based  on  recognizing  that,  at  any  given  frequency  on  the  Nyquist  plane,  there  is 
only  one  direction  of  perturbation  of  relevance  to  the  stability  analysis.  This  critical 
direction  is  defined  by  the  oriented  line  that  has  its  origin  at  the  location  of  the  nominal 
eigenvalue  of  the  unperturbed  system  and  passes  through  the  critical  point  -1+jO.  In  fact, 
from  a  stability  point  of  view  all  points  on  the  uncertainty  eigentemplate  that  do  not  lie  on 
the  critical  hne  can  be  ignored.  This  allows  the  characterization  of  robust  stability 
margins  for  uncertain  systems  characterized  by  irregular  perturbation  templates,  a 
problem  that  poses  significant  challenges  to  other  analysis  methods.  Other  attractive 
features  of  the  Nyquist  robust-stability  margin  theory  are  that  its  MIMO  (Mukiple  Input- 
Multiple  Output)  form  is  a  natural  extension  of  the  SISO  (Single  Input-Single  Output) 
case,  and  that  its  versions  for  continuous  and  discrete  systems  are  formally  identical.  A 
SISO  and  a  MIMO  example  are  given  illustrating  the  concepts. 


Chapter  3  focuses  on  the  complex  parametric  uncertainty  problem  from  the 
perspective  of  the  critical  direction  paradigm.  Considerable  attention  has  been  given  to 
the  problem  of  assessing  robust  stability  and  computing  stability  margins  for  SISO 
systems  with  parametric  uncertainties.  While  the  general  robust  stability  margin 
computation  problem  is  known  to  be  NP-hard  [Braatz  et  al.,  1994],  its  rank-one 
counterpart  results  in  a  convex  optimization  problem  that  can  be  readily  computed,  and 
can  be  often  derived  as  an  explicit  analytical  expression.  It  is  now  well  known  that  affine 
parametric  uncertainty  problems  result  in  a  rank-one  /i  formulation.  This  chapter 
addresses  the  case  of  element-wise  complex  parametric  uncertainties  where  each 
individual  uncertainty  lies  in  a  highly  structured  domain  that  is  convex  along  any  line 
from  the  origin.  This  is  a  mild  convexity  condition,  since  the  actual  domain  could  be 
highly  non-convex.  Exact  analytical  results  are  derived  for  geometrically  simple 
uncertainty  sets  such  as  ellipses  and  rectangles.  This  approach  however,  is  applicable  to  a 
much  wider  class  of  systems  than  the  ones  considered  here.  The  only  information  that  is 
needed  is  the  phase  dependent  magnitude  of  the  uncertainty.  The  approach  presented 
here  is  simple,  intuitive  and  yet,  mathematically  rigorous.  For  the  special  case  of  disk- 
bounded  affine  uncertainties,  the  approach  recovers,  in  a  transparent  manner,  relevant 
results  by  previous  methods. 

Chapter  4  and  Chapter  5  deal  with  the  issue  of  stabilizing  an  uncertain  time-delay 
systems  (both  state-delay  and  input-delay  are  considered).  Time  delay  exists  in  various 
branches  of  engineering.  From  biological,  optical,  electric  networks,  to  chemical  reaction 
systems,  time  delay  occurs  and  affects  the  stabiUty  and  performance  of  a  wide  variety  of 
systems  [See  Schell  and  Ross,  1986,  Inamdar  et  ai,  1991].  For  example,  input  delays 
occur  as  one  of  the  common  sources  of  time  delay,  such  as  in  chemical  processes, 
transportation  lags  and  measurement  delays,  etc.  It  is  well  known  that  the  existence  of 
time  delay  degrades  the  control  performances  and  makes  the  closed-loop  stabilization 
difficult.   Hence,  the  study  of  time  delay  systems  has  been  paid  considerable  attention 


over  the  past  years.  The  description  of  time-delay  systems  lead  to  differential-difference 
equations,  which  require  the  past  values  of  the  system  variables.  There  are  no  known 
methods  to  get  necessary  and  sufficient  conditions  for  the  robust  stability  of  such 
systems.  There  are  many  approaches  that  yield  sufficient  conditions  with  varying  degrees 
of  sophistication. 

In  Chapter  4,  a  Variable  Structure  Control  (VSC)  design  methodology  is  proposed 
to  robustly  stabilize  an  uncertain  state-delay  system  with  nonlinear  and  possibly 
mismatched  uncertainties,  utilizing  the  concept  of  perturbation  compensation.  The  study 
of  uncertain  state-delay  systems  has  received  much  attention  by  researchers  in  the  past 
decade.  One  reason  for  this  interest  is  that  time  delay  is  encountered  in  various 
engineering  systems  and  can  be  the  cause  of  instability.  Another  reason  is  that  practical 
control  systems  unavoidably  include  uncertainty  or  disturbances  due  to  modeling  or 
measurement  errors  and  other  approximations.  The  sliding  mode  in  VSC  possesses  well- 
known  features  that  make  it  very  attractive  for  control  systems.  These  include  fast 
response,  insensitivity  to  parameter  variations,  and  decoupling  design  procedure,  among 
others.  Sufficient  robust  stability  conditions  are  derived  which  are  independent  of  the 
size  of  the  delay;  furthermore  no  matching  conditions  on  the  uncertainty  are  assumed. 
The  main  advantages  of  the  proposed  approach  are  that  i)  a  realizable  control  law  is 
obtained,  ii)  chattering  reduction/elimination  is  achieved,  and  in)  guidelines  are  given  for 
the  design  of  the  control  matrix  involved  in  the  specification  of  the  switching  function. 
Finally  the  results  are  illustrated  with  an  open-loop  unstable  example. 

Chapter  5  proposes  an  approach  for  uncertain  input-delayed  systems  with  linear  and 
possibly  mismatching  uncertainties  via  Sliding  Mode  Control.  Sufficient  robust  stability 
conditions  are  derived.  We  consider  a  transformation  to  convert  the  original  system  into 
an  equivalent  system  without  delay  [Artstein,  1982].  This  reduction  enables  the  use  of 
known  stabilizing  methods  for  linear  systems  without  input  delay.  The  main  advantages 
of  the  proposed  approach  are  that  (i)  a  realizable  control  law  is  obtained,  (ii)  chattering 


reduction  is  acliieved  since  a  continuous  control  law  is  used,  (iii)  guidelines  are  given  for 
tiie  design  of  the  switching  function  to  ensure  the  asymptotic  stability  of  the  normal 
system  in  the  sliding  mode.  Certain  open  issues  in  the  sliding  mode  literature  are  brought 
out  and  possible  solutions  are  discussed. 

Finally,  future  research  directions  and  possible  extensions  to  the  presented  work  are 
discussed  in  chapter  6. 


CHAPTER  2 

THE  NYQUIST  ROBUST  STABILITY  MARGIN--A  NEW  METRIC  FOR  THE 
STABILITY  OF  UNCERTAIN  SYSTEMS 


2.1    Introduction 

Many  of  the  now  well-known  results  on  robust  stability  can  trace  their  origins  to 
the  Nyquist  stability  criterion  (SISO  case),  or  the  Generalized  Nyquist  stability  criterion 
(MIMO  case)  [MacFarlane,  1980).  Such  is  the  case  of  the  well  known  multivariable 
stability  margin  /:,„  [Safonov,  1982]  and  the  structured  singular  value  jii  [Doyle,  1982]. 

Although  the  generalized  Nyquist  criterion  addresses  the  issue  of  stability  using  complex 
variable  (eigenvalue  or  transfer  function)  arguments,  the  predominant  robustness  analysis 
methods  typically  rely  on  arguments  based  on  real  singular  value  or  magnitude  bounds. 
One  reason  for  this  choice  was  clearly  the  more  advantageous  conditioning  of  singular 
values  for  numerical  calculations.  However,  another  important  factor  was  the  historical 
fact  that  uncertainty  measurements  were  initially  proposed  for  MIMO  systems  using 
singular  value  bounds,  which  collapse  to  magnitude  (disk)  bounds  for  SISO  systems. 

In  this  chapter  we  introduce  the  concept  of  the  Nyquist  robust-stability  margin, 
^yy,  for  characterizing  the  closed-loop  stability  of  uncertain  systems.    The  approach 

makes  direct  use  of  Nyquist  arguments  and  is  based  on  the  analysis  of  the  perturbed 
eigenvalue  loci,  hence  avoiding  undue  conservatism  that  may  result  through  the  use  of 
singular- value  upper  bounds. 

A  key  element  in  the  new  approach  is  Critical  Direction  Theory  applied  to 
uncertainties  in  the  Nyquist  plane.  A  basic  version  of  the  critical  direction  concept  was 
first  presented  in  [Guzzella,  et  ai,  1985]  as  a  tool  for  analyzing  the  robust  stability  of 
polynomials  with  ellipsoidal  uncertainties.  In  [Latchman  and  Crisalle,  1995]  the  critical 


direction  notion  is  redefined  directly  on  the  Nyquist  plane  to  facilitate  the  derivation  of 
necessary  and  sufficient  robust  stability  conditions  for  transfer  functions  subject  to 
arbitrary  perturbations.  In  this  paper  we  further  extend  the  critical  direction  theory  for  the 
analysis  of  MIMO  systems  by  formulating  the  concept  in  the  generalized  Nyquist  plane, 
introducing  new  nomenclature  to  facilitate  the  seamless  extension  from  the  SISO  to  the 
MIMO  case,  and  defining  the  new  concept  of  the  Nyquist  robust  stability  margin. 

The  critical  direction  method  is  based  on  recognizing  that,  at  any  given  frequency 
on  the  Nyquist  plane,  there  is  only  one  direction  of  perturbation  of  relevance  to  the 
stability  analysis.  This  critical  direction  is  defined  by  the  oriented  line  that  has  its  origin 
at  the  location  of  the  nominal  eigenvalue  of  the  unperturbed  system  and  passes  through 
the  critical  point  -1+jO.  In  fact,  from  a  stability  point  of  view  all  points  on  the  uncertainty 
eigentemplate  that  do  not  lie  on  the  critical  line  can  be  ignored.  This  allows  the 
characterization  of  robust  stability  margins  for  uncertain  systems  characterized  by 
irregular  perturbation  templates,  a  problem  that  poses  significant  challenges  to  other 
analysis  methods.  Other  attractive  features  of  the  Nyquist  robust-stability  margin  theory 
are  that  its  MIMO  form  is  a  natural  extension  of  the  SISO  case,  and  that  its  versions  for 
continuous  and  discrete  systems  are  formally  identical. 

The  chapter  is  organized  as  follows.  Section  2.2  gives  a  succinct  review  of 
robust-stability  margin  results  for  the  SISO  and  MIMO  cases.  Section  2.3  formulates  the 
critical  direction  theory  for  analyzing  the  robustness  of  uncertain  SISO  systems,  and 
Section  2.4  gives  the  details  of  the  new  critical  direction  theory  for  the  MIMO  case.  A 
SISO  and  a  MIMO  example  are  given  in  Section  2.5,  followed  by  final  remarks  and 
conclusions  in  Section  2.6. 

2.2     Frequency-Domain  Approach  to  Robust  Stability  Analysis 

The  new  Nyquist  robust  stabiUty  margin  is  an  object  derived  using  Nyquist 
arguments.   In  order  to  provide  a  contextual  background  for  the  ensuing  discussions,  in 
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this  section  we  review  relevant  SISO  and  MIMO  results  on  robust  stability  analysis, 
placing  particular  attention  on  the  Nyquist  arguments  that  also  lie  at  the  root  of  the  new 
developments  presented  here. 

The  transfer  function  g{s)  of  a  SISO  system  can  be  written  in  terms  of  an  additive 
perturbation  5{s)  about  a  nominal  transfer  function  ^o(^) 

8{s)^go{s)  +  S{s),  5{s)ed  (2.I) 

where  d  represents  the  set  of  allowed  perturbations.  When  all  8{s)  e  d  are  considered, 
then  at  each  frequency  the  map   g{jco)  =  gQ{j(o)  +  5{jco)  defines  a  region  denoted 

uncertainty  template  or  value  set  [Barmish,  1993].  It  is  then  possible  to  visualize  the 
transfer-function  uncertainty  directly  in  terms  of  the  classical  Nyquist  plot  as  uncertainty 
templates  about  the  locus  of  the  nominal  system  gQ{jco). 

The  mathematical  description  for  the  set  of  perturbations  may  be  formulated  in 

different  spaces.  For  example,  the  Laplace-function  domain  d  is  often  defined  implicitly 

in  terms  of  a  frequency-domain  description.    Such  is  the  case  where  the  uncertainty 

templates  are  specified  to  be  circular  at  each  frequency.    This  makes  the  analysis  of 

robust  stability  particularly  straightforward.    On  the  other  hand,  the  templates  for  a 

number  of  useful  uncertainty  descriptions  are  not  circular  but  rather  irregularly  shaped, 

exhibiting  highly  directional  features.  This  is  the  case  for  parametric  uncertainty 
descriptions  of  the  form  S{s)  =  g{s,p)-g[s,pQ),  where  g{s,pQ)  and  pQ,  respectively, 

denote  the  nominal  transfer  function  and  the  nominal  parameter  set,  and  where  g{s,p)  is 
characterized  by  m  real,  uncertain  parameters  p  edp  a  %'" .  In  this  framework  the  set  d 
is  defined  implicitly  through  the  Euclidean  domain  dp.  Typical  uncertainty 
representations  for  the  parameter  uncertainty  set  dp  are  polytopes  and  ellipsoids.  The 
special  case  where  the  uncertainty  templates  are  ellipses  is  an  interesting  example  of 
directional  uncertainty  templates  which  also  has  a  tractable  mathematical  representation. 

The  representation  (2.1)  is  general  and  adequately  encompasses  additive 
perturbations  as  well  as  uncertainties  that  may  appear  as  multiplicative  perturbations  of 
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the  nominal  system,  as  parametric  variations  on  pole/zero  locations,  and/or  variations  of 
the  coefficients  of  the  numerator  and  denominator  polynomials  in  the  transfer  function. 
The  difference  in  each  case  is  the  shape  and  orientation  of  the  uncertainty  templates 
associated  with  5{s). 

Of  particular  interest  are  three  types  of  SISO  uncertainty  descriptions. 
Specifically,  we  denote  the  set  of  uncertainties  giving  rise  to  unstructured  (circular), 
elliptical,  and  structured  (arbitrary)  templates,  respectively,  by  the  notation  du,  d^,  and 

ds- 

In  the  case  of  uncertain  MIMO  systems,  the  matrix  transfer  function  G{s) can  be 
written  as  an  additive  perturbation  A{s)  about  a  known  nominal  transfer  matrix  Gg{s) 

G{s)  =  G,{s)  +  A{s),  A{s)eD  (2.2) 

where  D  is  the  set  of  allowed  perturbations.  The  uncertainty  set  D  is  described  in  terms 
of  available  information  about  the  modeling  errors.    We  use  the  designation  D^  to 

represent  Unstructured  uncertainties  with  a  single  norm  bound,  D^  to  denote  block- 
diagonal  norm-bounded  uncertainties,  Dq  to  represent  element-wise  circular  (disk- 
bounded)  uncertainties  and  D^  to  represent  element-by-element  elliptical  uncertainties. 
Other  uncertainty  descriptions  are  defined  in  terms  of  errors  in  the  elements  of  the  state- 
space  matrices,  or  errors  in  the  real  parameters  of  the  matrix  transfer  functions.  In  all 
these  cases  the  uncertainty  description  still  can  be  mapped  (via  eigenvalue  inclusion 
regions)  on  to  frequency-response  uncertainty  templates,  in  the  Nyquist  plane  which  will 
then  form  the  basis  for  defining  the  Nyquist  robust  stabihty  margin. 

Classical  SISO  Robust-Stability  Results 

The  SISO  robust  stability  analysis  is  concerned  with  the  stability  of  the  system 
that  results  when  the  system  g{s)  in  (2.1)  is  arranged  in  a  unity  negative-feedback 
configuration.   It  is  normally  assumed  that  (/)  the  nominal  system  is  stable  under  unity 
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negative-feedback,  and  that  (ii)  the  nominal  and  uncertain  system  have  the  same  number 

of  open-loop  unstable  poles.  These  assumptions  are  adopted  throughout  this  paper. 

The  most  studied  case  is  perhaps  that  where  the  uncertainty  description  is 
unstructured,  5{s)edu  i.e.,  when  \5{j(o)\<\W{jco)\,  where  W{s)  is  a  known  function 

whose  frequency-response  magnitude  defines  the  radius  of  the  circular  uncertainty 
templates. 

Invoking  the  Nyquist  stability  criterion  leads  to  the  necessary  and  sufficient 
stability  requirement  that  all  uncertainty  templates  exclude  the  critical  point  -1-hjO,  i.e., 

8{j(o)  =  go{jco)  +  5{jco)^-l  Vd{jco)  and  \/q)  ^2.3) 

This  condition  can  be  expressed  through  the  inequality  [Doyle  et  ai,  1992] 

Wjjco) 


or  equivalently, 

W{j(0) 


<1  Vg) 


(2.4) 


<1 


(2.5) 


l  +  ^oW 

The  robustness  analysis  for  uncertainty  descriptions  with  non  circular  templates  is 
more  challenging.  For  this  reason  it  is  quite  common  for  elliptical  or  arbitrarily  shaped 
uncertainty  templates,  arising  for  example  from  parametric  uncertainties,  to  be 
circumscribed  by  an  appropriate  circle  [Bhattachrya  et  ai,  1993].  Although  this  approach 
yields  sufficient  conditions  for  robust  stability,  it  is  nevertheless  inherently  conservative. 
Classical  MIMO  Robust-Stability  Results 

As  in  the  SISO  case,  the  robustness  analysis  for  MIMO  systems  makes  use  of  the 
generalized  Nyquist  stability  criterion  [MacFarlane,  1980],  where  avoidance  of  the 
critical  point  1-f-jO  is  at  the  center  of  interest  from  the  point  of  view  of  absolute  stability 
assessment.  This  has  been  the  major  focus  of  much  research  interest  in  the  development 
of  robust  multivariable  stability  margins.  Consider  the  uncertain  transfer  matrix  G(s)  in 


11 


(2.2)  with  nominal  model  G^is)  affected  by  an  uncertainty  A(s).  Again  it  is  assumed 
that  (0  the  nominal  system  is  stable  under  unity  negative-feedback,  and  that  (ii)  the 
nominal  and  uncertain  system  have  the  same  number  of  open-loop  unstable  poles. 

The    conditions    under   which    no    eigenvalue    of   the    uncertain    system 
Gg{jQ})  +  A{jo))  is  equal  to  -1+jO,  i.e., 

X{Go{jco)  +  A{jco))^-l  Vzi(7ffl)  and  Vo)  ^2.6) 

is  readily  shown  to  lead  to  the  determinantal  condition 

det(l  +  M{jco)A{jco))  ^  0  Vzi(;ffi>)  and  \/co  ^2  i\ 

where  M{jco)  :=  (l  +  Gq  (jffl))"^ . 
Condition  (2.7)  is  of  course  equivalent  to  the  eigenvalue  condition 

^M{jco)A{jco))^-l  VA{joj)  and  Vco  ,281 

The  multivariable  stability  margin  k„j  and  the  structured  singular  value  ji  were 

independently  defined  by  Safonov  [1982]  and  Doyle  [1982],  respectively,  in  terms  of  the 

determinantal  stability  condition  (2.7),  rather  than  the  eigenvalue  conditions  (2.6)  or 

(2.8).  For  an  uncertainty  of  a  given  description  D,  the  multivariable  stability  margin  k„^ 

is  defined  as  the  matrix  2-norm  of  the  smallest  destabilizing  uncertainty  in  the  given 
class,  namely: 

k„t{o})  =  mmla[A{jQ))) :   det(l  -f-  M{ja))A{jco))  =  o} 

This  definition  corresponds  to  the  reciprocal  of  the  robustness  measure  /j.  defined  in 
[Doyle,  1982],  which  is  defined  as  the  inverse  of  the  norm  of  the  smallest  destabilizing 
uncertainty  in  the  class,  namely: 

Li(co)  =  nnn\a(A(jco)) :  det(l  +  M(j(o)A(jco))  =  O] 

AeD^    V    vv    ^;  V  \J    J    \J    I)       J  ^2.10) 
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Note  that  if  no  uncertainty  in  the  allowable  class  destabilizes  the  system  (i.e.,,  makes  an 
eigenvalue  of  Go{jco)  +  A{jcd)eqml  to  -1  at  any  frequency)  then  k„^  =  00  and  /i  =  0. 

These  margins  provide  a  measure  of  tolerable  uncertainty  size. 

A  number  of  powerful  robust-stability  results  can  be  obtained  in  terms  of  the 
robust  stability  margins  (2.9)and  (2.10).  Consider,  for  example,  the  following  case  of 
interest.  An  uncertainty  description  D  is  said  to  be  dosed  under  contraction  and  rotation 
if  for  any  A{s)  g  D,  then  ye^^A{s)  e  Dior  all  0  <  7  <  1  and  for  all  0<9<  2k..  Starting 
from  the  eigenvalue  stability  condition  (2.8),  we  can  state  the  following  lemma. 

Lemma  2.1 

Consider  an  uncertainty  description  D  that  is  closed  under  contraction  and 
rotation.   Then  the  uncertain  system  G{s)  =  G^{s)  +  A{s)  is  stable  for  allA{s)  €  D  if  and 

only  if 

sup  p[M{jco)A{jco))  <  1  Vo) 

AbD  (2.11) 

Proof:     Since  the  eigenvalue  condition(2.8)  is  necessary  and  sufficient  for  robust 

stability,  it  is  clear  that  (2.11)  immediately  constitutes  a  sufficient  stability  condition. 

That  (2.11)  is  also  necessary  is  established  by  contraposition  as  follows:  Consider  an 
uncertainty  Ag{s)sD  such  that   p(M2l^)  =  [7]~  >1,  then  due  to  the  closure  under 

contraction  we  can  always  find  an  uncertainty  Ai  e  Dwith  A^  =yA^\  0<7<1,  such 
thatp(M/4] )  =  1 .  Since  the  phase  of  the  uncertainty  class  is  arbitrary,  due  to  the  closure 
under  rotation  we  can  assign  a  scalar  phase  multiplier  to  A^^  to  get  Ar^  such  that 
X{MAqf)  =  -\,  which  of  course  implies  instability.     Thus,  it  follows  that  for  the 

description  D  specified  in  the  Lemma,  condition  (2.11)  is  both  necessary  and  sufficient 

for  robust  stability.  V 

For  the  classes  of  uncertainties  covered  by  Lemma  2.1,  the  spectral  radius 

condition  (2.1 1)  is  in  fact  equivalent  (after  appropriate  normalizations)  to  the  definitions 
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of  fi  and  /:,„  and  thus  it  naturally  follows  that  /i  is  equal  to  the  left  hand  side  of  (2.1 1). 
Hence  the  conclusion  of  Lemm.a  2. 1  can  be  equivalently  formulated  as 

K^>lVco  (2.12) 

or 

^<l\/co  (2.13) 

The  A;„j  and  {J.  stability  margins  serve  to  provide  valuable  characterization  of 
destabilizing  uncertainties.  Specific  procedures  and  results  have  been  developed  for 
computing  these  stability  margins  (or  at  least  good  upper  and  lower  bounds)  for  various 
useful  classes  of  uncertainties.  Additionally,  the  jj.  problem  formulation  has  produced  an 
extensive  set  of  results  associated  with  the  class  of  block  diagonal  bounded  uncertainties, 
and  has  thus  attained  wide  acceptance  for  robustness  analysis,  and  more  recently,  also  for 
controller  synthesis  using  the  ^-synthesis  method.  In  all  these  applications  the 
computation  of  the  stability  margin  is  effected  by  calculating  singular-value  upper  bounds 
which  give  sufficient,  and  in  some  cases  necessary  and  sufficient,  stability  conditions. 

Unfortunately,  however,  the  use  of  norm-bounds,  especially  singular  value 
bounds,  in  the  characterization  of  destabilizing  uncertainties  often  causes  the  rich 
structural  properties,  including  phase  and  directionality  to  be  ignored.  On  the  other  hand, 
as  we  show  below  significant  advantages  can  be  gained  by  studying  the  robustness 
problem  directly  from  an  eigenvalue  point  of  view,  especially  since  necessary  and 
sufficient  stability  conditions  correspond  precisely  to  an  eigenvalue  condition  via  the 
generalized  Nyquist  criterion. 

2.3    The  Critical  Direction  Method  for  SISO  Systems 

In  this  section  we  present  an  exact  robust  stability  result  for  the  case  of  SISO 
systems  using  the  critical  direction  theory.  We  adopt  a  nomenclature  that  permits  the 
extension  of  the  critical-direction  concepts  to  the  MIMO  case.   We  also  introduce  the 
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definition  of  the  Nyquist  robust  stability  margin  for  SISO  systems  in  a  form  which 
readily  extends  to  the  MIMO  case. 

Figure  2.1  shows  a  typical  Nyquist  diagram  for  a  SISO  system  illustrating  the 
nominal  frequency  response  go{jOi)  and  an  irregularly  shaped  uncertainty  template.  We 

define  the  critical  line  at  a  given  frequency  co  as  the  directed  line  which  originates  at  the 
nominal  point  go{jco]  and  passes  through  the  critical  point  -1+jO.    The  figure  is  also 

useful  for  identifying  the  entities  defined  below: 

1.  The  crtitical  direction 

J,.    X  l  +  ^„(;ft)) 

|1  +  ^«U«)|  (2.14) 

which  may  be  interpreted  as  the  unit  vector  that  defines  the  direction  of  the  critical  line. 

2.  The  uncertainty  template 

1(0))  :=  {gijco)  I  gijco)  =  g^{jco)  +d(jQ}),  6(s)  e  d}  ^2.15) 

3.  The  critical  template 

%{co)  :=  {z  €  T(ffl)  I  z  =  goijco)  +0(d{jw)  for  some  aeSi^  ] 

namely,  the  set  of  perturbed  frequency-responses  lying  along  the  critical  direction 
d{jco). 

4.  The  critical  perturbation  radius 

Pr((o)  •=  max  (a  U  =  gJjco)  +cxd(j(o)  g  T^(co)] 

As  an  illustration  of  the  previous  definitions  note  that  the  critical  line  at  frequency 
CO  =  Q}i  is  readily  identified  in  Figure  2.1  as  the  directed  line  with  origin  at  gg[jco^)  and 

passing  through  the  point  -1+jO.   The  critical  direction  d[ja)\)  is  simply  the  unit-length 
vector  that  characterizes  the  direction  of  the  critical  line.   In  addition,  the  critical  radius 
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pc(ft>i)  is  shown  as  the  distance  between  the  nominal  Nyquist  point  go{j(Oi)  and  the 

point  where  the  critical  line  intersects  with  the  boundary  of  the  template.  Finally,  the 
critical  template  %{o))  is  also  readily  characterized  in  Figure  2.1  as  the  subset  of  the 
uncertainty  template  1(0))  that  intersects  with  the  critical  line. 

By  definition  %{co)  is  a  subset  of  the  critical  line,  thus  it  follows  that  %{co)  is 
either  a  single  straight-line  segment,  or  the  union  of  such  segments.  The  critical  template 
may  also  contain  isolated  points  should  the  boundary  of  the  template  T(ft))  be  tangent  to 
the  critical  direction.  Figure  2.1  shows  the  case  where  the  critical  template  is  a 
continuous  segment,  and  hence  %{q))  is  a  convex  set  even  though  the  entire  template  is 

highly  non-convex.  The  critical  template  is  discontinuous  when  it  is  made  up  of  the 
union  of  distinct  segments.    In  this  case  1(.{co)  is  not  convex,  but  each  of  its  member 

segments  is  a  convex  subset. 

With  these  definitions  in  hand  we  can  proceed  to  state  the  following  robust- 
stability  theorem.  For  ease  of  exposition  we  assume  that  the  critical  template  %ico)  is  a 

convex  set  for  all  frequencies.  This  restriction  can  be  relaxed  through  obvious 
modifications  to  account  for  each  of  the  convex  segments  of  %{o)). 

Theorem  2.1 

Let  the  nominal  SISO  system  gg{s)be  subject  to  an  uncertainty  5{s)  e  d.    Then  the 

uncertain  closed-loop  system  remains  stable  under  unity  feedback  if  and  only  if 

I /~^ — VJ<^   ^^ 

\^  +  M8o{m)\  (2.18) 

The  proof  of  the  theorem  is  omitted  since  all  details  are  given  in  [Latchman  and 
Crisalle,  1995].  V 

Motivated  by  (2.18)  we  now  propose  the  following  definition  of  the  Nyquist 
Robust-Stability  Margin  for  SISO  systems: 
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From  Theorem  2. 1  it  follows  that 


kf^{co)<l  Va  ^2.20) 

is  a  necessary  and  sufficient  condition  for  robust  closed-loop  stability.  Furthermore,  the 
quantity  y  =  [k;^{a))]     specifies  the  amount  by  which  the  uncertainty  template  should  be 

increased  or  decreased  to  attain  the  limiting  case  of  stability. 

The  SISO  critical  direction  result  enables  an  exact  assessment  of  stability  for 
systems  with  directional  uncertainty  templates,  without  having  to  circumscribe  the 
templates  with  larger  circular  uncertainties.  The  calculation  of  the  critical  perturbation 
radius   pc{(o)  involves  determining  the  intersection  of  a  straight  line  (the  critical 

direction)  with  a  curve  (the  boundary  of  the  template  'T{co)).    Section  2.5  shows  an 
example  where  the  critical  radius  can  be  calculated  analytically. 

It  is  also  worth  noting  that  the  critical  perturbation  radius  also  gives  a  systematic 
methodology  for  determining  the  uncertainty  weights  which  can  be  used  as  one  of  the 
inputs  for  various  robust  synthesis  methods.  Of  particular  interest  in  this  regard  are  cases 
where  there  exists  a  Critical  Weighting  Function,  W^{s),  such  that  its  magnitude  satisfies 

the  interpolation  condition 

\%{4  =  Pcioi)  \/(D  ^2.21) 

When  (2.21)  can  be  exacdy  satisfied.  Theorem  2.1  leads  to  the  famiUar  H^  condition 

W,{s) 


<1 

(2.22) 


which  is  of  the  form  (2.4).  If  (2.22)  cannot  be  satisfied  exactly  but  can  be  approximated 
using  standard  frequency-domain  regression  methods,  then  the  resulting  approximate 
weight  can  still  be  used  as  the  basis  for  practical  robust-synthesis  design  for  systems  with 


highly  structured  templates.  | 

t 
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It  should  be  noted  that  Theorem  2.1  recovers  as  a  special  case  the  situation  where 
the  uncertainties  are  circular  or  disk-bounded.  The  interested  reader  is  referred  to 
[Latchman  and  Crisalle,  1995]  for  further  details  on  this  point. 

In  the  following  section  we  show  how  the  critical  direction  approach  and  the 
Nyquist  robust-stability  margin  may  be  generalized  to  the  case  of  uncertain  MIMO 
systems. 

2.4    The  Critical  Direction  Metliod  for  MIMO  Systems 

For  MIMO  systems  we  focus  attention  on  the  effect  of  uncertainties  on  the 
eigenvalue  loci  along  the  critical  direction.  This  treatment  has  several  desirable 
properties.  First  of  all  we  derive  necessary  and  sufficient  robust  stability  conditions 
based  upon  eigenvalue  relationships  obtained  directly  from  eigenvalue  uncertainty 
templates.  Thus,  in  principle  this  method  yields  a  stability  margin  measure  even  in  cases 
where  singular  value  conditions  fail  to  give  necessary  and  sufficient  stability  conditions, 
provided  that  a  method  can  be  found  to  define  tight  inclusion  regions  for  the  uncertain 
eigenvalues.  Furthermore  by  considering  only  the  subset  of  the  eigenvalue  inclusion 
region  that  lies  along  the  critical  direction,  we  significantly  reduce  the  computation 
involved  in  using  the  boundary  of  the  eigentemplates  (See  E-contours  [Daniel  and 
Kouvaritakis,  1985] )  for  stability  assessment. 

Consider  the  generalized  Nyquist  plots  in  Figure  2.2  which  (for  illustrative  i 

I 

purposes  only)  show  two  irregularly  shaped  eigenvalue  inclusion  regions  (templates) 
about  two  of  the  n  eigenvalues  Xi[Gg{jco)),  i  =1,  2,...,n.  In  analogy  with  the  SISO  case, 

for  each  of  the  n  eigenloci  we  define  at  each  frequency  an  associated  critical  line,  defined 
as  the  directed  line  originating  at  the  location  of  nominal  eigenvalues  X^[G^{j(X)))  and 

passing  through  the  critical  point  -1+jO.  Also  in  analogy  with  the  SISO  case  discussed  in 
Section  2.3  we  define  the  following  entities: 
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1.  The  crtitical  directions 

•^'    '        |l  +  Ai(G„(»)|  (2-23) 

which  may  be  interpreted  as  the  unit  vectors  which  define  the  direction  of  each  critical 
line. 

2.  The  uncertainty  eigentemplates 

r^{co)  :=  {^G{jco))\  ^G{jco))  =  A,(G,(»  +  A{jco)),  A(s)  e  d]    ^2.24) 

3.  The  critical  eigentemplates 

Tci{a>)'={zeTi{cd)\z  =  ^G^{jco))  +  ad^{jco)  for  some  aeSi^  }(2.25) 

4.  The  critical  perturbation  radii 

p, .  {co)  :=  max  (a,-  \\z  =  X,  (G,  {jco))  +  a,^-  {jco)  g  %■  {(o)]  .3  26) 

The  entities  defined  above  have  identical  interpretations  to  their  SISO 
counterparts  defined  in  Section  2.2.  Hence,  their  geometrical  interpretation  can  be 
obtained  directly  from  the  SISO  Nyquist  diagram  given  in  Figure  2.1,  provided  that  the 
nominal  frequency-response  plot  goijco)  is  substituted  by  an  eigenvalue  plot  Gg{ja). 

Obviously,  the  tight  E-contour  templates  for  unstructured  and  structured  uncertainty 
descriptions  utilized  in  [Daniel  and  Kouvaritakis,  1985;  Kouvaritakis  and  Latchman, 
1985]  are  respectively  equivalent  to  the  uncertainty  templates  1^(tt»). 

Using  these  definitions  we  can  now  state  the  following  robust  stability  theorem. 

As  in  the  SISO  case,  for  simplicity  of  exposition  we  assume  that  the  critical  uncertainty 
eigentemplates  T^^ico)  are  convex  sets  at  all  frequencies. 
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Theorem  2.2 

Let  the  nominal  MIMO  system  G^{s)  be  subject  to  an  uncertainty  A{s)  e  D.   Then 
the  uncertain  closed-loop  system  remains  stable  under  unity  feedback  if  and  only  if 

Pc,i(^)  p,(ft)) 

max    I —/ — - — -,  =  I ^-fH^ — :7T  <  1  Vct)  /^  ^-yx 

/=i,2,...n|l  +  l,(G,(ya;))|     |l  +  A,(G,(»)|  (2.27) 

where  p^{.)and\{.)  are  respectively  used  to  denote  the  critical  perturbation  radius  and 

the  eigenvalue  associated  with  the  eigenloci  resulting  from  the  maximization  over  all 
i  =  1,2,- --n  in  (2.27). 

Proof.  Assuming  nominal  closed  loop  stability  and  that  the  nominal  and  perturbed 
open  loop  systems  have  the  same  number  of  open  loop  unstable  poles,  the  Generalized 
Nyquist  stability  criterion  guarantees  that  the  uncertain  closed  loop  system  is  stable  if  and 
only  if 

Xi(G^[jco)  + A[ja))^  ^  -1  V/  and  \/co  r?  28) 

Using  the  defining  equation  for  the  eigenvalues  of  Gq{j(o)  +  A{jco),  namely 

dti(G^[ico)-zJ  +  A{i(o))  =  Q  ^2.28) 

at  each  frequency  co  we  can  parametrize  the  uncertain  eigenvalues  by 

z  =  ^/fe(;«))  +  P,-e^^'  (2.29) 

where  d^  =  9i{co)  varies  in  the  range  0<0i<2k,  and  for  each  value  of  9^  the  scalar 
Pi  =  pi{co)  varies  in  the  range  0  <  p ■  <  p ■  where  2,i[G^{jco))  +  p-e-'  '  is  an  eigenvalue 
corresponding  to  the  boundary  of  the  eigentemplate  Ti{cl)).  Then  condition  (2.28)  can  be 
written  as 

l  +  A,(G,(;fo))  +  AV^'    ^   0  V/  (2.30) 

Since  the  term  1  + A,- (0^(70}))  is  fixed,  the  only  possibility  for  violating  the  stability 
condition  (2.30)  is  for  p-e^  '  to  be  oriented  along  the  critical  direction,  namely 
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p^e^^' ^  ai{co)di{jco)  (2.31) 

where  0  <  o:;(g))  <  p^,..  Thus,  for  robust  stability  we  have  the  necessary  and  sufficient 
stability  condition 

l  +  ?ii[G^{jco))  +  ai{a))di{jO))^  0  V/  and  \/co  ('2  32) 

Invoking  definition  (2.23)  it  follows  that  a  sufficient  stability  condition  is  given  by 

\^^UG,{jco))\  (2.33a) 

Furthermore,  since  along  the  critical  direction  a^{co)  <  p^^{(o)  stability  is  ensured  if 

Pc  (^) 

I 7 -7   <  1  \/i  and  \/co  ^  qqu\ 

\l  +  ^G,{jco))\  (2.33b) 

The  proof  is  completed  by  taking  the  maximum  over  all  i  and  noting  that,  because  of  the 
convexity  assumption  on  the  critical  templates  T^^{co)  conditions  (2.33a)-(2.33b)  are  also 

necessary  for  stability.  V 

In  complete  parallel  with  the  SISO  case 

|l  +  A,(G,(j«))|  (2.34) 

defines  the  MIMO  Nyquist  Robust-Stability  Margin,  and  the  scalar  7  =  [^y^ (&))]"    again 

specifies  the  amount  by  which  the  uncertainty  eigentemplate  should  be  increased  or 
decreased  to  attain  the  limiting  case  of  stability.  Consequently,  from  Theorem  2.2  it 
follows  that  a  necessary  and  sufficient  condition  for  robust  stability  is  given  by  the 
Nyquist-derived  constraint 

k^{co)<\  \fco  (235) 

The  Nyquist  robust  stability  condition  (2.35)  provides  an  exact  answer  to  the  robust 
stability  margin  problem,  much  in  the  same  vein  as  the  margins  /:„j('^)  ^nd  li{co). 


21 

A  most  pleasing  feature  of  the  new  stability  measure  is  that  the  SISO  version 
(2.20)  and  the  MIMO  version  (2.35)  of  the  Nyquist  robust  stability  margin  are  formally 
identical.     Much  more  importantly,  the  kf^{(a)  formalism  adheres  to  the  classical 

generalization  of  SISO  system  properties  to  the  MIMO  case  via  the  vehicle  of  the 
eigenvalues  of  the  frequency  response  matrix. 
Relationship  between  A:,„(fi)),/i(ft)),  and  ki^{(D) 

In  light  of  Theorem  2.2  and  the  remarks  above,  it  is  clear  that  there  are  strong 
equivalences  between  the  Nyquist  stability  margin  k!^{co)  and  the  margins  A:,„(ft>))  and 
/J.{co).  To  see  this  recall  that  the  definition  (2.9)  for  k„^{a))  is  equivalent  to  finding  at 
each  frequency  the  smallest  destabilizing  A  such  that 

l  +  ^M{j(o)A{jco))  ^0  V/  (2.36) 

which  is  precisely  the  condition  exploited  in  Theorem  2.2  to  get  the  necessary  and 
sufficient  stability  condition 

l  +  Xi[G^{jco))  +  ai{co)di{jco)^  0  Vr  and  Va  Q  37) 

For  uncertainties  which  satisfy  the  condition  of  Lemma  2.1,  from  (2.12)-(2.13)  it  follows 
that 

ki^{co)  <  1  o  /i(ft))  <  1  <=>  k,„{co)  >  1  <^  sup  p{MA)  <  1 

4eD  (2.38) 

In  fact,  even  for  uncertainties  for  which  Lemma  2. 1  does  not  hold,  it  is  still  true  that  for 

all  A(s)  e  D 

l  +  ^M{jco)A{jco))  ^0^k^{o))<\  (2.39) 

The  practical  utility  of  Theorems  2.1  and  2.2  requires  the  computation  of  the 
critical  perturbation  radius  pdco).  Several  methods  are  presently  being  developed  to 
exploit  the  critical  direction  theory  to  compute  pdo))  and  hence  kj^{a))  for  the  case  of 
affine  and  multi-affine  parametric  uncertainties  in  SISO  and  MIMO  systems.  The  results 
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thus  far  are  very  encouraging,  and  have  shown  significant  computational  and  algorithmic 
simplifications  due  to  the  focus  on  the  critical  direction.  For  the  case  of  structured  and 
unstructured  MIMO  uncertainties,  tight  eigenvalue  inclusion  regions  may  be  obtained 
using  the  singular-value  based  E-contours  method  [Daniel  and  Kouvaritakis,  1993; 
Kouvaritakis  and  Latchman,  1985],  with  similarity  or  non-similarity  scaling  deployed  to 
remove  or  reduce  conservatism. 

2.5     Examples 
SISO  Example 

This  section  presents  an  example  of  a  SISO  system  with  an  uncertainty 
description  motivated  from  statistical  parameter-estimation  techniques.  The  uncertainty 
considered  is  an  ellipsoidal  parameter-space  model  which,  in  addition  to  its  intrinsic 
merits,  facilitates  the  direct  analysis  in  the  Nyquist  plane  permitting  an  explicit 
characterization  of  the  Nyquist  robust-stability  margin  ki^{co), thus  illustrating  the 

application  of  the  critical  direction  theory  proposed  in  this  paper.  Furthermore,  a 
discrete-time  representation  of  the  model  is  deliberately  chosen  to  emphasize  the  fact  that 
the  critical  direction  method  and  its  attendant  Nyquist  robust-stability  margin  concept  are 
applicable  to  both  the  continuous  and  discrete  domains. 

Let  us  consider  a  general  SISO  system  model  given  by  the  discrete-time  transfer 
function 

H(z;p)  =    hkz'^ 

k  =  1  (2.40) 

with  q  parameters  defined  by  h  =  \hl,1^2,■■■h  \  .  The  nominal  system,  denoted  H^(z),  is 
obtained  when  h  assumes  the  nominal  values  h^=  hi,h^,---h^     ,  so  that  the  real 

system  is  modeled  as  H{z)  =  H^{z)  +  5H{z).    Let  h  in  '^^  be  such  that  h  =  h^  +  dh. 
Then  we  define  the  uncertain  ellipsoidal  parametric  uncertainty  description 
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Sh  =  h  -  h^    G   Df^  (2.41) 


D,,   =  [Sh  eSi"^  s.t  5h''Qj;'5h;Q,=Qj;>0] 


(2.42) 

We  argue  that,  in  addition  to  the  remarkable  mathematical  tractability  which  we 

shall  show  later,  the  ellipsoidal  parametric  description  of  uncertainties  is  quite  natural  in 

many  applications  and  offers,  in  contrast  to  hyper-rectangular  descriptions,  the  further 

advantage  of  allowing  the  dependence  among  various  system  parameters  to  be  taken  into 

explicit  account.  Ellipsoidal  models  often  arise  quite  naturally,  as  for  example,  whenever 

linear  regression  or  least-squares  analysis  is  used  in  model  estimation  [Guzzella,  et  al, 

1991;Kosut6'ra/.,  1992]. 

Clearly  the  nominal  parameter  vector  h^  defines  the  nominal  transfer  function 

Hoiz)  -     J^Kz'^ .  Consider     now     the     frequency     response      H{e^^) 

k=  1 


"         -k 


.0    ,0         ,  0 


where /f(z)  =     t.h^       contains  parameters  h^^  hi,h2,---h^       belonging  to  the 

k=  1 

parameter  ellipsoid.  Under  these  conditions  the  following  lemma  shows  that  the 
parameter  space  ellipsoid  maps  precisely  to  ellipses  at  all  frequencies  except  for  co  =  0 
and  co  =  7T.  Let  Re  (z)  and  Im  (z)  represent  respectively  the  real  and  imaginary 
components  of  a  complex  number  z. 

Lemma  2.2 

The  parameter  space  ellipsoid  defined  by  (2.42)  maps  to  the  elliptical  uncertainty 
template 


T{(o)  =  lz{co)  =  X,  (co)  I X2  (co)  =  0,1  |xi  {(d)  -  x'}{co)\  <  ^vi^{cofQi,v;^{co),co  =  0    co  =  k 


(2.43) 


where 


X(£0):=[x,(ffl)  X2{(0)]    e^'_  X^{co)-[x'^{oj)  x^i 


-0^-^    ^0(«; 


%3i2 


such  that 


X{co)- 


XoH 


Re  H(e^'^ 
Im  H  e^'^ 


Re  nJeJ'' 


Im  Hje^^ 


=  V{co)h 


=  V{aj)h, 


Q^    :=   F(«)(2,F(«f  e  9^2x2 


V{0))  :-- 


COSG)     coslco     •••     cosqO) 
sinft)     sin2a)     ■•■     sinqo) 
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A?  CO  =  0    a/trf  Q)  =  7r,  matrix  V{co)  is  singular  and  the  uncertainty  template  is  entirely 
real  and  given  by 


1{co)  =  lz{co)  =  x^ico)  I  X2ico)  =  0,1  |xi(ffl)  -  x^{(a)\  <  ^vi^{a)f  QhVn{(o),CD  =  0    a^n 


n 

where  v^(a))  :=  [cos  co    cos  2(0    ■■■    cos  qty] 


Proof.  The  proof  makes  use  of  a  result  given  in  [Guzzella,  et  al,  1985].  Details 
are  omitted.  V 

Our  intention  is  to  characterize  the  Nyquist  robust-stability  margin  for  the  unity 
feedback  system  comprised  of  the  nominal  discrete-time  system  and  its  associated 
parametric  uncertainty  ellipsoid  (2.42).  Note  that  the  robust  stability  of  the  closed-loop 
system  cannot  be  analyzed  using  the  /"  results  in  [Tsypkin  and  Polyak,  1991]  because 
their  method  is  applicable  only  to  continuous-time  systems  and  to  ellipsoidal  uncertainty 
descriptions  where  the  matrix  g/,  is  diagonal  {i.e.,  the  principal  directions  are  aligned 
with  the  coordinate  axes).  In  fact,  none  of  the  conventional  methods  for  analyzing  robust 
stability  of  SISO  systems  appears  capable  of  treating  in  a  systematic  fashion  the  case  of 
ellipsoidal  uncertainties  D/,  considered  in  this  example. 
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Theorem  2.3 

Under  the  assumptions  of  nominal  closed-loop  stability  and  that  the  nominal  and 
perturbed  systems  share  the  same  number  of  open-loop  unstable  poles,  the  unity  negative 
feedback  system  with  open-loop  transfer  function  (2.40)  and  uncertain  parameters  (2.42) 
is  stable  for  all  %  e   Di^  if  and  only  if 


where 


k^ico)   <   I  V  CO  ^2.44) 


k^io))  =     ,  _  for  CO  e  {0,7c) 

K{a)Qafd^{co)  (2.45) 


,^(,)  =  ,r^jff  f;7y  for  CO  =  M 

\     dc  («K(«)  (2.46) 


and 


d^{co)  :-- 


-I 
0 


-  V{cD)h, 


(2.47) 

Proof.  Inequality  (2.44)  is  a  direct  result  of  applying  the  necessary  and 
sufficient  condition  (2.18)  of  Theorem  2.1  in  the  Nyquist  robustness  margin  form  (2.20). 
It  then  suffices  to  prove  (2.45)  and  (2.46).   First  consider  the  case  where  co   e  {0,7t). 

Using  two-dimensional  vector  analysis  and  the  definitions  given  in  Section  3,  it  is  clear 
that  the  components  of  d^{Q})  given  by  (2.47)  are,  respectively,  the  unnormahzed  real 

and  imaginary  parts  of  the  critical  direction  d{jco)  as  defined  in  (4.1).  Now  from  Lemma 
2,  at  each  co  the  frequency-response  Hie-'^]  is  located  inside  an  elliptical  template  with 
center  at  the  point  hJc-^^).   Clearly,  the  elliptical  template  is  a  convex  set  because  it  is 

continuous  along  any  ray.  A  fundamental  result  from  two-dimensional  coordinate 
geometry  gives  the  length  of  the  line  joining  the  center  of  an  ellipse  to  the  point  of 
intersection  along  the  vector  d^{Q))  as 
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PcM  =  . 


(2.48) 


Recognizing  that 


l  +  Hje^'' 


dJ{Q))d,{co) 


(2.49) 


it  readily  follows  that 


k^{co)  = 


p,{Q)) 


l  +  Hje^'' 


(2.50) 


is  of  the  form  (2.45).   Analogously,  the  proof  for  the  singular  cases  «  =  0,;r  is  derived 
noting     that     from     Lemma     2.2     it     follows     that     at     these     frequencies 


Pc{o))  =  --^VRico)  QijVj^{oj)  from  which  (2.46)  is  readily  established. 

In  summary,  for  the  case  of  ellipsoidal  parametric  uncertainties,  the  critical 
direction  theory  permits  the  exact  characterization  of  the  Nyquist  robust-stability  margin, 
and  the  derivation  of  the  exact  necessary  and  sufficient  condition  (2.44)  for  robust 
stability.  V 

MIMO  Example 


Consider  the  MIMO  nominal  transfer  function 


G.(s)  =  0.5  (-lX^  +  2)(s  +  3)     (s  +  4) 
and  the  associated  element-by-element  elliptical  uncertainty  description. 
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D^    =  {A(s)  s.t.AikO'o))   e   E;^{jco)  V«} 


\ 


(2.51) 


where  the  boundary  5E,'^,  (jo))  of  each  elliptical  domain  E^^{jcl))  is  given  by  the  map 


J<t>ik 


dE,,{jco)  =   2A.,,.^^*cos0,,  +  j  IB^'^me,,,  i,k  =  1,2         ^2.52) 
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and  where  the  coefficients  A-^  and  B-j.  are  the  elements  of  the  frequency-dependent 
matrices  A{q))  and  5(a)),  and  the  major-axis  orientation  parameters  0,-^,  are  the  elements 
of  the  matrix  0(a)).  The  objective  is  to  illustrate  the  calculation  of  the  Nyquist  robust 
stabihty  margin  k^{co)  at  the  frequency  cOi  =1.21,  where 


A(a),) 


0.1264  0.0359 
0.0680  0.3185 


5(0)1 ) 


0.0246  0.0278' 
0.0537  0.2707 


0(a)i) 


2.1799    4.5759' 
1.4906     3.8865 


Note  that  formulations  based  on  the  multivariable  stability  margin  k^^ico)  or  the 

structured  singular  value  iJ.{a))  cannot  currently  offer  an  obvious  approach  for  calculating 

the  stability  margins  for  the  uncertainty  description  Dg. 

At  each  frequency  of  interest  we  propose  the  following  calculation  procedure. 
First,  the  nominal  eigenvalues  /l,(G„(ja)))  and  the  critical  directions ti; (o)), /  =  1,2  are 

computed. 

Second,  for  each  eigenvalue,  z  satisfying  Zi  =  \  +  pidi  is  found.  In  the  above 

expression,  D  is  a  diagonal  similarity  scaling  matrix,  P  is  the  diagonal  matrix  where 
each  diagonal  element  is  the  radius  of  the  circle  circumscribing  the  elliptical 
uncertainties.  E\  and  E2  are  derived  from  the  diagonalization  of  the  uncertainty 
[Kouvaritakis  and  Latchman,  1985]  and  A,-  represents  the  eigenvalues  of  G^.  This  yields 
the  upper  bound  for  A,(Go  -I-  zl)  in  the  critical  direction.  Then  an  optimization  is  carried 
out  over  all  Q^y.  minimizing  the  function 

2 


/=zf-A,(G,  +  zl) 


-H 


d, 


1  +  A,(G,+Zi) 
|l  +  A,(G,+zi)| 


(2.53) 

Equation  (2.53)  gives  the  eigenvalue  of  (G^  +  4)  in  the  critical  direction,  which  is  closest 
to  the  upper  bound  z° .  The  basic  principle  is  similar  to  that  described  in  [Kouvaritakis  et 

al,  1991]. 

At  the  selected  frequency  0)]  =1.21,  the  nominal  eigenvalues  for  the  system 
considered     in     this      example      are      \(Gg[ia}{)]  =  -0.2Q<d5- jLl^U,      and 
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?i2[Gg{j(jOifl=  1.9332-72.4525  and  have  the  associated  critical  directions 
^j(;a)i)  = -0.5446 +  70.8387  and  J2(M)  = -0.7672  +  70.6414.  Then  using  the 
structured  E-contour  method  for  the  fixed  critical  directions  d][jcOi)  and  d2{jC0\)  yields 
the  critical  radii  p^^  («,)  =  0.3140  and  p,,^  {co{)  =  0.6445.  Finally,  from  (2.27)  and  (2.34) 

it  is  found  that 


k^^  {(o)  =  max< 


PcM) Pcii^l) 


=  max{0.2139,0.1686}(2  54) 


^  l  +  Ai(Go(jft),))|'l  +  A2(Go(j«i)) 
Hence,  the  Nyquist  robust  stability  margin  is  k'^^[cOi)<  1,  and  it  is  concluded  that  the 
system  satisfies  the  necessary  and  sufficient  condition  for  robust  stability  at  this  particular 
frequency. 

Figure  2.3  shows  a  Nyquist  diagram  with  the  eigen-plots  for  the  nominal  system 
considered  in  this  example.  The  eigentemplates  Ti((Di)  and  7^2  (0)1)  are  not  shown  in  the 

figure  because  for  the  purpose  of  determining  the  Nyquist  robust-stability  margin  it  is  not 
necessary  to  calculate  the  entire  templates.  However,  for  reference  the  figure  shows  the 
Optimal  D-scaling  bounds  obtained  using  the  scalings  given  in  [Kouvaritakis  and 
Latchman,  1985].  The  actual  eigentemplates  are  bounded  by  these  Optimal  D-scaling 
contours.  The  points  denoted  z, (7(1)1)  and  Z2{jC0\)  are  on  the  boundary  of  the  critical 
templates  %  (ft),)  and  %^_{cO]),  respectively,  and  were  obtained  using  the  E-contour 

method.  Note  that  for  this  example  the  Optimal  D-scaling  bounds  provide  a  good 
estimate  for  the  boundary  point  Zi{jO)i),  but  yield  a  poor  estimate  for  22(^^1)  •  Note  also 
that  the  critical  eigentemplates t;  [co^]  and  T^^_{cO\)  can  be  readily  identified  from  the 
figure  as  the  straight-line  segments  joining  each  nominal  eigenvalue  /L,' (0^(70)1))  with  its 
corresponding  boundary  point  Zi{jC0i). 

2.6     Conclusions 

In  this  paper  we  have  proposed  the  Nyquist  Stability  Margin,  k^  (ft))  as  a  new 

metric  for  robustness  analysis  of  SISO  and  MIMO  systems.    The  definition  of  new 


29 


stability  margin  is  based  on  the  Critical  Direction  Theory  which  provides  a  single 
framework  for  robustness  analysis  for  SISO  and  MIMO  systems.  The  analysis 
methodology  makes  direct  use  of  the  generalized  Nyquist  diagram,  and  in  contrast  to  the 
prevalent  approaches  which  emphasize  singular-value  perturbations,  it  focuses  attention 
on  eigenvalue  perturbations. 

The  main  advantage  of  the  critical-direction  theory  is  that  it  provides  necessary 
and  sufficient  conditions  for  robust  stability  in  the  presence  of  highly  structured 
uncertainties  with  phase  and  directionality  constraints.  Other  approaches  to  these 
problems  either  do  not  have  the  inherent  capability  to  deal  with  these  structural  details,  or 
the  directionality  and  phase  constraints  are  deliberately  ignored,  giving  rise  in  either  case 
to  sufficient-only  conditions  such  as  those  associated  with  singular-value  theory.  On  the 
other  hand  the  new  method  explicitly  exploits  the  detailed  directionality  and  phase 
constraints  of  the  uncertainties  as  these  are  manifested  in  the  frequency  domain 
uncertainty  templates.  Thus,  the  new  method  is  applicable  to  a  number  of  uncertainty 
descriptions  for  which  other  methods  fail,  such  as  the  case  of  element-by-element 
ellipsoidal  uncertainties  in  the  transfer-function  matrix,  and  other  uncertainty  descriptions 
with  highly  directional  frequency-domain  templates. 

The  new  critical-direction  technique  opens  up  new  avenues  for  robustness 
analysis  and  could  lead  to  novel  approaches  for  robust  control  synthesis.  There  is 
significant  promise  for  fruitful  new  results  in  this  area  where  the  computational  efforts 
are  concentrated  on  a  single  and  well-defined  frequency-dependent  directed  line. 
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Nomenclature  for  Chapter  2 


s 

Laplace  variable 

d,{jQ)) 

direction  in  the  complex  plane 

SoiJOi) 

Nominal  SISO  transfer  function 

gijo)) 

Uncertain  SISO  transfer  function 

Goijco) 

Nominal  MIMO  transfer  function 

G{jco) 

Uncertain  MIMO  transfer  function 

kp^ico) 

Nyquist  robust  stability  margin 

kjco) 

Multivariable  robust  stability  margin 

M{jco) 

Augmented  system  matrix 

WAs) 

SISO  weight 

Greek  Letters 

CO 

Frequency 

^Gijco)) 

'i'th  eigen  value 

Ais) 

Uncertainty  description 

Pc(®) 

Critical  perturbation  radius 

ll{(0) 

Structured  singular  value 

a{.) 

Maximum  singular  value 

Pi') 

Spectral  radius 

r{co) 

Eigen  template 

%{co) 

Critical  eigen  template 
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■  critical  line 


Re  gjco) 


Figure  2.1.  An  irregularly  shaped  uncertainty  template  at  frequency  cOi 
(shaded  area)  and  its  critical  perturbation  radius  Pc(ct)i).  The  critical 
template  %{cOi)  (solid  line)  is  the  subset  of  template  points  lying  on  the 
line  segment  with  end  points  go(^i)  ^^^^  gg[(0])  +  Pc{cOi)d{jcOi). 
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6^2  (M 


Im  X{G,{jco)) 


Re  X{gXjco)) 


h{G,{jo))) 


Figure  2.2.  Nyquist  plot  of  the  eigenvalues  of  a  2x2  MIMO  system 
showing  two  irregularly  shaped  eigen  templates  'T^[co^)  and  ^{cO]) 
(shaded  areas)  at  frequency  C0  =  C0i. 
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Figure  2.3.  Nyquist  plot  of  the  MIMO  Example.  The  contours  represent 
the  optimal  D-scaling  bounds  for  the  eigentemplates  T,  [co^ )  and  T2(ft>i )  at 
the  frequency  ©1=1.21  .  The  points  Zi(cOi)  and  22(^1)  li^  °^  the 
intersection  of  the  boundary  of  their  respective  eigentemplates  with  the 
critical  direction 


CHAPTER  3 

A  NEW  PERSPECTIVE  ON  COMPUTING  ROBUST  STABILITY  MARGINS  FOR 
COMPLEX  PARAMETRIC  UNCERTAINTIES 

3.1     Introduction 


Considerable  attention  has  been  given  to  the  problem  of  assessing  robust  stability 
and  computing  stability  margins  for  SISO  systems  with  parametric  uncertainties.  It  is 
now  well  known  that  affine  parametric  uncertainty  problems  result  in  a  rank-one  '" 
formulation.  While  the  general  robust  stability  margin  computation  problem  is  known  to 
be  NP-hard  [Braatz  et  ai,  1994],  its  rank-one  counterpart  results  in  a  convex  optimization 
problem  that  can  be  readily  computed,  and  can  be  often  derived  as  an  explicit  analytical 
expression.  This  aspect  of  the  problem  is  amply  documented  in  the  literature  [Qiu  and 
Davison,  1989]-[Hinrichsen  and  Pritchard,  1992]  [Chen  et  ai,  1994a]-[Chen  et  ai, 
1994b].  Of  particular  interest  are  references  [Chen  et  al,  1994b]  and  [Young,  1994], 
which  establish  connections  between  frequency-based  stability  conditions  and 
polynomial  conditions  in  the  spirit  of  Kharitonov's  theorem  [Kharitonov,  1979]. 

This  chapter  focuses  on  the  complex  parametric  uncertainty  problem  from  the 
perspective  of  the  recently  proposed  critical  direction  paradigm  [Latchman  and  Crisalle, 
1995]  [Latchman  et  ai,  1997].  The  approach  followed  in  critical  direction  theory 
consists  of  first  mapping  the  parameter  space  uncertainties  into  the  Nyquist  plane  in  the 
form  of  value  sets  or  uncertainty  templates,  and  then  invoking  specific  directionality  and 
phase  properties  to  infer  stability.  The  method  is  based  on  recognizing  that  at  any  given 
frequency  on  the  Nyquist  plane,  the  only  uncertainties  of  relevance  to  the  stability 
analysis  are  those  which  lie  on  a  well  defined  critical  direction.  In  fact,  from  a  stability 
perspective  all  points  on  the  uncertainty  template  that  do  not  lie  in  the  critical  direction 
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can  be  ignored.  This  approach  makes  it  possible  to  solve  problems  with  highly 
directional  uncertainties,  such  as  SISO  elliptical  uncertainties,  for  which,  at  this  time, 
approaches  based  on  the  more  traditional  structured  singular  value  methods  are  not 
obvious. 

This  chapter  addresses  the  case  of  element-wise  complex  parametric  uncertainties 
where  each  individual  uncertainty  lies  in  a  highly  structured  domain  that  is  convex  along 
any  line  from  the  origin.  This  is  a  mild  convexity  condition,  since  the  actual  domain 
could  be  highly  non-convex.  Exact  analytical  results  are  derived  for  geometrically  simple 
uncertainty  sets  such  as  ellipses  and  rectangles.  For  the  special  case  of  disk-bounded 
affine  uncertainties,  the  approach  recovers  in  a  transparent  manner  relevant  results  by 
previous  methods  [Hinrichsen  and  Pritchard,  1992][Chen  et  al,  1995b]. 

Section  3.2  of  the  chapter  provides  background  information  and  mathematical 
preliminaries,  presents  the  uncertainty  descriptions  considered,  and  gives  a  brief  review 
of  the  main  elements  of  the  critical  direction  theory.  The  main  results  are  presented  in 
Section  3.3,  and  connections  with  the  classical  approach  based  on  M-D  structures  is 
explored  in  Section  3.4.  Concluding  comments  are  given  in  Section  3.5. 

3.2    Background  and  Preliminaries 

Consider  a  SISO  system  with  an  open-loop  transfer  function 

n{jco,q)  ,-  IX 

8ijC0,q)  =  --. r  (3-1) 

dijCO,q) 

where  n{j(0,q)  and  dijco,q)  are  complex  polynomials  that  depend  on  an  uncertainty 
vector 

q  =  [q,{jC0),q2{j0j),-,q,n{jC0)f^Q'^C'"  (3.2a) 

whose  components  qi^iJQ))  belong  to  complex  domains  2^.,  i.e., 

qk{jC0)GQi^(zC,  k  =  l,2,---,m  (325) 
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and  hence  the  uncertainty  domain  is  the  cartesian  product  space  (2  =  2i  x  ^2  ^ ' ' '  ^  2™  • 
Figure  3.1  illustrates  a  case  of  a  radially  convex  domain  Qj^.  Note  that  the  entire  domain 
is  non-convex. 

The  uncertainty  class  considered  in  this  s  are  described  by  element-wise 
uncertainty  sets  Qj^,  k=  I,  2,  ...,  m,  that  are  closed  under  contraction  (i.e.,  are  radially 
convex),  and  are  unconstrained  in  phase.  Letting  (92^,  denote  the  boundary  of  domain  (2^, 
then  the  boundary  elements  dqj^  (jco)  g  dQi^  are  of  the  form 

dq^=r^((0,ei,)e'^',  k   =    l,2,  ■••,  m  (3.2c) 

where  0  <  9^  <  2;r  is  an  unconstrained  phase  angle,  and  7^(tt),0^)  =  [^/.(y©)]  >  0  is  the 

frequency-  and  phase-dependent  magnitude  of  the  boundary  element. 

A 

-''^k  r(t)     ^ 


y(t) 


Figure   3.1      Radially  convex  Figure  3.2  Unity  negative 

uncertainty    domain    2^     for  feedback  configuration, 

complex  component  qj^ 

The  elements  of  the  uncertainty  vector  q  are  assumed  to  appear  in  an  affine 
fashion  in  (3.1),  i.e.,  the  uncertain  model  is  of  the  form 


n(jCO,0)+  lni^U(O)qkU0^) 


8Uo),q) 


k=l 


diJQ),0)+  ^di^UG))qkU0}) 


k=i  (3.3) 

Under  the  unity  negative-feedback  arrangement  shown  in  Figure  3.2,  the 
characteristic  polynomial  for  system  (3.3)  is  given  by  the  polynomial  sum 


pijCO,q)  =  nijCO,q)  +  dijco,q) 


(3.4) 
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It  follows  that  the  robust  stability  of  the  closed  loop  with  respect  to  the  uncertainty  set  Q 

can  be  evaluated  by  studying  the  uncertain  polynomial  (3.4).     It  is  convenient  to 
decompose  (3.4)  into  its  nominal  and  uncertain  parts   p(ja),q)  =  Pg(jco)  +  5{jcO,q), 

where 

p^ijco)  =  n(jco,0)  +  d(jco,0)  p5g) 

and 


dijco,q):=  'ZPk{jC0)qk{J0^)  (3-5b) 

k=l 


where 


The  variable    6  Ajco)  is  simply  the  frequency-dependent  phase  of  polynomial  pj^{jco). 
Hence, 

m 
p{jC0,q)^p^ijCO)+  ^Pi^{jO))qi^{jCO)  (3.7) 

k=l 

Definition.  The  Parametric  Robust  Stability  Margin  is  the  entity  defined  by  the 
minimization  expression 

a  (o))=  min  \a\  0  €  p(jco,aq)    for  some  qeO] 

» 

The  parametric  robust  stability  margin  is  a  non-negative  real  scalar  that  can  be 
interpreted  as  the  minimal  magnification  a''{co)>l  or  contraction  a'{co)<l  of  the 
uncertainty  set  Q  that  brings  the  closed-loop  system  to  the  edge  of  stability. 
Geometrically,  the  parametric  robust  stability  margin  represents  the  minimum  tolerable 
blow-up  factor.  Note  that  the  parametric  robust  stability  margin  is  defined  for  each 
frequency. 
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Specific  Parameter  Uncertainty  Descriptions  Considered 

Three  types  of  uncertainty  descriptions  of  the  form  (3.2c)  are  considered,  namely 
the  (/)  circular,  (ii)  elliptical,  and  (Hi)    rectangular  element-wise  complex  domains 


illustrated  in  Figure  3.3. 
(a)  Im  cji. 

A 


^Req,^ 


(b)         lin  q 


Figure  3.3   Three  types  of  element-wise  complex  parametric  uncertainty 
regions  with  unconstrained  phase:  (a)  circular,  (b)  elliptical,  and  (c) 
rectangular  domains. 
Circular  Uncertainties 

As  depicted  in  Figure  3.3a,  this  is  the  familiar  case  of  circular  (disk-bounded) 
uncertainties  of  radius  A^.((y)  where 


r,{co,9,)  =  A,ico)  \fe, 


(3.9) 


Elliptical  Uncertainties 


Let  A^,(co)>0,  Bj.{a))>Q,  and  (f)i^{(o)  respectively  represent  the  frequency- 
dependent  semi-major  axis,  semi-minor  axis,  and  the  orientation  with  respect  to  the  real 
axis  of  the  an  ellipse  centered  at  the  origin  of  the  complex  plane  as  shown  in  Figure  3.3b. 
Then  the  boundary  of  each  elliptical  domain  is  given  by  the  map 

^^,  (0))  =  A,  (a)).^'^^('"W(0,-0,(a)))  +  ;5,(ffl)g^'^^-('"W(  0,-0,(0))) 
from  which  it  follows  that 


rk(o^A)=  A,{cDfco,\e,-<p,{co))  +  B,{cD)hm\9,-<p,{cD)y-    (3.10) 
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Obviously,  the  circular  uncertainty  case  (3.9)  is  recovered  by  (3.10)  after  setting 
A^{(o)  =  B,{a)). 

Rectangular  Uncertainties 

Let  Aj^{(o)>0  and  Bi^{co)>0  respectively  represent  the  frequency-dependent  half- 
width  and  the  half-height  of  a  rectangle  centered  at  the  origin  of  the  complex  plane  as 
shown  in  Figure  3c.  Consider  the  orientation-phase  definition 

0j«):=arctan  ^^4^  (3.11a) 

and  the  associated  phase  sets 

0bk-=l'Pk^^-'PM^+h^^^-(t>k)  (3-1  ic) 

The  shaded  and  the  plain  areas  in  the  Figure  3.3c  corresponds  to  0;,^.  and  0^^, 
respectively.  Then  an  analytical  expression  for  the  magnitude  of  any  point  on  the 
boundary  of  the  rectangular  regions  can  be  compactly  written  in  the  form 


Yki^A) 


for  9j^  e  0^1^ 


|cos(0^)| 

for  d/^  e  0^^ 


|sin(0,)| 


(3.12) 


3.3.    The  Critical  Direction  for  a  Characteristic  Polynomial 

The  critical  direction  theory  proposed  for  rational  systems  in  [Latchman  et  ai,  1997]  can 
be  readily  modified  for  the  case  of  polynomials  through  the  adoption  of  the  definitions 
given  below. 

(/)      The  Critical  Direction 

4y^):=_^4z4:.,A.M  (3.13) 

\PoUm 
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is  the  unit  vector  wliich  defines  the  direction  from  the  nominal  point 
Po{j(o)  towards  the  origin.  The  critical  direction  is  uniquely  identified  by  the 
critical  phase  angle  0^.{co). 
(ii)  The  Uncertainty  Template  (or  value  set) 

r{(D):={p{jco,q)eC\pijCO,q)  =  p^{jco)  +diJQ),q),  qsQ]       (3.14) 
(Hi)  The  Critical  Template 

%{<^)  '=  {MM^)  e  '2"((y)  I  pU(0,q)  =  PoU<^)  +rd{j(0),  for  some  r  e  9t+}(3.15) 
(iv)  The  Critical  Perturbation  Radius 

Pcico)  '=  max  [r\z  =  Po{J^)  +^d{i(o)  e  %{co)^  (3.16) 

(v)  The  Nyquist  Robust  Stability  Margin 

k^{co):=j^j^  (3.17) 

Note  that  the  Nyquist  robust  stability  margin  ki^{Q})  defined  in  (3.17)  is  a  metric 

that  characterizes  the  distance  on  the  Nyquist  plane  to  the  point  of  instability  {i.e.,  the 
point  l+jO).  In  a  complementary  fashion,  the  parametric  robust  stability  margin  a  {(o) 
defined  in  (3.8)  characterizes  the  "distance"  on  the  parameter  space  Q  to  the  limit  of 
stability. 

3.4.    Main  Results 

This  section  develops  a  technique  for  deriving  analytical  expressions  for  the 
parametric  robust  stability  margin.  The  approach  makes  use  of  the  critical-direction 
theory  elements  presented  in  Section  3.2.  The  developments  are  specialized  here  for  the 
case  of  uncertain  polynomial  systems. 

Lemma  3.1.    Consider  the  uncertain  characteristic  polynomial  (3.4),  and  assume 
that  nominal  stability  is  attained  under  unity  feedback,  and  that  the  critical  templates 
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T(&))  are  convex  at  all  frequencies.     Then  the  uncertain  system  remains  stable  for  all 
parameter  uncertainties  q  eQ  if  and  only  if 

kf^{co)<   1    Vffl  (3.18) 

Proof.  The  proof  follows  from  an  application  of  the  zero  exclusion  principle 
and  is  analogous  to  the  development  given  in  [Latchman  et  al,  1997]  for  transfer 
functions,  and  details  are  omitted  for  brevity.  V 

Note  that  Lemma  3.1  is  quite  general  and  holds  for  multiaffine,  and  complex 
and/or  real  uncertainties.  The  major  issue  is  the  calculation  of  the  critical  radius  Pc(ft)) 
needed  for  the  characterization  of  kj^{a))  in  (3.17);  however,  as  shown  below  in 
Theorem  3.1,  the  computation  of  the  critical  radius  for  complex  affine  parametric 
uncertainties  is  particularly  straightforward. 

Theorem  3.1.  Consider  the  uncertain  system  (3.3)  with  the  uncertainty  description 
(3.2a)-(3.2c),  and  its  corresponding  closed-loop  uncertain  characteristic  polynomial 
(3.7).  Then  the  critical  perturbation  radius  is  given  by  the  expression 

Pc{(^)=  S  \Pk{j(o)\  7k[(0,eM-QpS'^))  (3.19) 

Proof.  The  points  p{j(o,q)  on  the  value  set  'T{(X))  that  belong  to  the  critical 
template  %{(o)  must  satisfy  the  condition 

p{jco,q)  =  p^ijco)  +  5(jco,q)  =  p^ija)  +  rd{j(o)  (3.20) 

or 

ip,{jco)q,{jco)  =  reJ'^^^'^  (3.21) 

Due  to  the  radial  convexity  of  the  uncertainty  set,  it  suffices  to  consider  only  boundary 
uncertainties  (9^^ (ft));  then,  utilizing  (3.2c)  the  preceding  equality  can  be  rewritten  in  the 

form 
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ZkM7<.(«,e.)^'''''*''"  "'''""  =  '•  (3.22) 

k=l 

From  the  definition  (3.16)  it  follows  that  Pc{(o)  is  equal  to  the  maximal  value  of  r  in 
(3.22)  over  all  possible  uncertainty  phases  0  <  0/,  <2k,  i.e., 

p,{(0)=    max    Y\Pk{j(o)\rk{coA)  ^^       "  '  (3-23) 

O<0i.<2;r/;=i 
Clearly,  the  maximum  is  obtained  for  0^,  =6^{ct))-0p  (co),  leading  to  the  expression 

(3.19).  V 

Theorem  3.2.    Consider  the  uncertain  system  (3.3)  with  the  uncertainty  description 
(3.2a)-(3.2c),  and  the  corresponding  critical  perturbation  radius  Pc(ft)).     Then  the 

parametric  robust  stability  margin  is  given  by 

a\co)  =  ^  (3.24) 


Proof.  From  equation  (3.4) 


pijCO,q)  =  p^ijO))+Y,Pi^{jCo)qi^{jO)),    qeQ  (3.25) 


and 


pij(a,aq)  =  p^(ja>)  +  a^p,^{j(D)  qi,{jco),   qeQ  (3.26) 

k=i 

Note  that  the  map  (3.25)  yields  the  uncertainty  template  'T{(j))  defined  in  (3.14).  As  for 
the  question  of  stability,  note  that  only  those  uncertainties  that  give  rise  to  points 
pijO),q)  that  are  aligned  with  the  critical  direction  are  relevant.  Then,  using  the 
definition  (3.13)  for  the  critical  direction,  all  the  points  of  relevance  for  the  stability 
analysis  must  be  of  the  form 

pijOJ,q)  =  Poijco)  +  rd{jco)  (3.27) 

Equating  (3.25)  and  (3.27)  yields 
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I^PkU^)  qk{J(^)^rd{j(o)  (3.28) 

■Note  that  if  q^^eQ  is  such  that  p{j(i),q)  is  aligned  with  the  critical  direction,  then  for 
any  scalar  a>0,  (3.26)  implies  that  p{jco,aq)  also  is  aligned  with  the  critical  direction. 
Hence,  making  use  again  of  the  definition  of  the  critical  direction  given  in  (3.13),  and 
from  (3.25)  and  (3.27)  it  is  concluded  that 

p(ja),aq)^ p^(jco)  +  a  rd{jco)  =  d{j(o)(ar-\p^(jco)\)  (3.29) 

from  which  if  clearly  follows  that 

pijco,aq)  =  0  «  «  =  MM  (3.30) 

r 

Since  by  definition  the  maximum  possible  value  of  r  is  Pc(co)  ,it  follows  that  a^ico)  is 
obtained  when  r  =  p^{co)  and  hence  a" {co)  =  \pg{ja>)\/ p^{co) .  V 

Combining  Theorem  3.1  and  Theorem  3.2,  it  follows  that  the  parametric  robust 
stability  margin  (3.8)  can  be  written  in  the  form 

«*(«)  =  ^r -^ (3.3.) 

k=\  ^  >^      I 

Hence,  (3.31)  gives  an  analytical  expression  for  the  parametric  robust  stability  margin  for 
all  complex  element-wise  parametric  uncertainties  that  are  closed  under  contraction  of  the 
form  (3.2a)-(3.2c)  for  systems  with  the  affine  structure  (3.3). 

3.5.    Parametric  Robust  Stability  Margins  for  Highly  Structured  Uncertainties 

The  connection  between  the  Nyquist  robust  stability  margin  k^{(0)  and  the 
parametric  robust  stability  margin  a  {(o)  is  the  reciprocal  relationship 

a{a))  =  — ^  (3.32) 

which  follows  from  Theorem  3.2  and  the  definition  (3.17).  Invoking  Lemma  3.1,  it  then 
concluded  that  the  closed-loop  system  is  robustly  stable  if  and  only  if  a  (o))  <  1,  Vft). 
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This  result  holds  for  all  element-wise  complex  parametric  uncertainty  descriptions  that 
are  closed  under  contraction.  This  section  focuses  on  obtaining  analytical  expressions 
for  parametric  robustness  margin  for  the  three  specific  types  of  uncertainty  descriptions 
depicted  in  Figure  3.3 

Theorem  3.3  Consider  the  uncertain  system  (3.3)  with  the  characteristic 
polynomial  {3.7),  and  the  circular,  elliptical,  and  rectangular  complex  parametric 
uncertainty  descriptions  whose  element-wise  magnitudes  are  respectively  given  by  (3.9), 
(3.10)  and  (3.12).  The  robust  stability  margins  for  each  of  these  descriptions  are  given 
by  the  analytical  expressions: 

(3.i)  For  the  case  the  circular  uncertainty  description 

\PoiJ(0)\ 


a  (ft))  = 


S  4(^)  k(j«)l 


(3.33a) 


(ii)     For  the  case  the  elliptical  uncertainty  description 


a  (&»)  =  ■ 


PoO'^) 


2„;„2 


X   A^,(a;)-cos2  0c(«)-^pA(fi^)-'^;t(«) h-^-tCft^)  sin^  0cM-^p/.(«)-'/'/t(«)    "  k(j«)| 


(Hi)    For  the  case  the  rectangular  uncertainty  description 

boO'«)| 


a' {(o)  = 


m 
I 


\Mpk{m 


^ak  + 


bMpA'M 


(3.33b) 


-(3.33c) 


sm[e,{co)-e,((D) 


'bk 


COs(0^(G))-0/,^(ffl) 

where  5^^.  =  1  if  Oj^  e  0^,^,,  otherwise  d^j^  =  0,  cu%d  where  Si^^-  =1  if  Oj^  g  Oj^j^, 
otherwise  Sj^i^  =  0 . 
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Proof.  Equations  (3.28b)-(3.28a)  simply  follow  by  substituting  the  expressions 
for  the  phase-dependent  uncertainty  magnitudes  (3.9),  (3.10),  and  (3.12)  into  the 
analytical  expression  (3.27)  for  the  parametric  robust  stabihty  margin.  V 

The  generalized  structured  singular  value  approach  [Chen  et  al,  1995b]  has 
previously  been  used  to  derive  the  result  (3.33a)  for  the  case  of  circular  uncertainties.  In 
principle,  the  general  formalism  presented  in  [Chen  et  al,  1995b]  can  also  be  used  to 
derive  the  results  (3.33b)  and  (3.33c);  however,  the  critical  direction  theory  extracts  exact 
analytical  results  through  very  simple  and  straightforward  geometric  arguments. 

3.6     Connections  with  the  Classical  M-A  Formalism 

In  this  section  expression  (3.31)  for  the  parametric  robust  stability  margin  is 
rederived  utilizing  the  framework  of  the  classical  M-A  formalism  where  robust  stability 
is  characterized  by  a  determinantal  constraint.  The  approach  follows  well-established 
developments,  except  that  we  simplify  the  derivations  by  invoking  a  critical-direction 
concept,  namely,  that  a  specific  orientation  in  the  Nyquist  plane  is  the  only  direction  of 
relevance  for  stability  analysis 

Lemma  3.2.  The  unity  negative-feedback  configuration  of  the  SISO  uncertain  system 
(3.3)  with  affine  uncertainty  elements  can  be  transformed  into  an  equivalent  MIMO 
M-A  structure  via  a  linear  fractional  transformation  [Zhou  et  al,  1996]  where  the 
uncertainty  matrix  A{s)  is  diagonal  and  the  structural  matrix  M  =  [Mij^)  is  of  rank- one 
and  has  rational  elements  given  by  Mij^{s)  =  — - — ,  /  =1,2,  ■■■,  m;  k  =1,2,  ■■■,  m. 

Proof.  The  proof  is  an  exercise  in  block-diagram  algebra.  Details  are  given  in 
the  Appendix.  V 

For  the  MIMO  M-A  structure  in  question  it  is  typical  to  define  the  parametric 
robust  stability  margin  in  terms  of  a  determinantal  condition,  as  follows 

a^ico)  :=     min  la  I  det(/  -I-  aM(jcy)A(j(y))  =  0  |  (3.34) 
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which  is  formally  analogous  to  the  definition  of  the  parametric  robust  stability  margin  for 
a  SISO  system  given  in  (3.8).  In  fact,  the  theorem  below  shows  that  the  definitions  (3.8) 
and  (3.34)  lead  to  identical  expressions  for  the  parametric  robust  stability  margin  for  the 
case  where  the  uncertainties  appears  affinely  in  the  model. 

Theorem  4.  Consider  the  system  (3.3)  with  its  associated  affine  uncertainty 
description  (3.2)  and  suppose  that  the  nominal  system  is  stable  under  unity  negative - 
feedback,  and  that  the  nominal  and  the  uncertain  systems  have  the  same  number  of  open- 
loop  unstable  poles.  Then  the  parametric  robust  stability  margin  at  each  frequency  is 
given  by 

a  {(o)  = ■ ^^^^^^^ ^ (3.35) 


k=\ 


Ydd,{a)-d    {(o)-K]  \pj,{]CO 


Proof.  From  Lemma  3.1,  the  uncertain  system  with  unity  negative  feedback  can 
be  transformed  into  an  equivalent  MIMO M-  A  staicture.  Using  standard  arguments,  the 
robust  stability  condition  for  the  transformed  MIMO  system  is  given  by  the  determinantal 
condition 

det(/ -  M{](o)^[]Co))^Q   Vco  (3.36) 

Introducing  the  parametric  robust  stability  margin  a  e  9?"^  interpreted  as  a  blow-up  factor 

for     the     uncertainties,     the     determinantal     condition     is     rewritten     as 
det(/  -  aM(jft))A(j£o))?iO.    Since  from  Lemma  3.L  matrix  M  is  of  rank  one,  it  is 

straightforward  to  establish  the  inequality 


det(/  -  aM{io))A{iQ)))=  1 


+  ■ 


a 


k=\ 


(3.37) 


Since  the  aim  is  to  find  the  minimum  scalar  ae%'^  that  makes  (3.37)  identically  zero, 
equating  (3.37)  to  zero  and  solving  for  a  yields 
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a  = 


p  fj«)|e^^>""'^) 


'i\pAJCoU[coA)e^''^'"'^ 


(3.38) 


k= 
Clearly,  the  minimum  real  non-negative  value  of  (3.37)  will  be  realized  when  all  terms  in 

the  summation  in  the  denominator  (3.37)  are  aligned  with  respect  to  each  other  and  with 
respect  to  the  numerator,  i.e.,  using  the  notation  p  [jco]  =  p  [jco]  e    ''"  and  recalling  the 

fact  that  the  uncertainty  phase  is  0^  unrestricted,  the  minimum  is  attained  when 

dj^  =  Op^^  (CO)  -  9^^  {co)  +  7i,  k  =  l,2,--;m  (3.39) 

From  the  definition  (3.13)  it  follows  that^^  (co)  =  9^{co) - 7Z ,    hence  the  minimizer  is 

realized  at 

di^=9,{co)-dp^{co),  k^i,2,--,m  (3.40) 

Substituting  (3.40)  into  (3.38)  gives  the  minimal  value  a  -  a'{(x))  which  is  of  the  form 
(3.30).  V 

3.7     Conclusions 

In  this  chapter  we  demonstrated  the  significance  and  ease  of  application  of  the 
critical  direction  theory  to  the  complex  affine  uncertainty  problem.  The  development  is 
simple  and  intuitive  ,  and  for  the  special  cases  of  interest  recovers  the  exact  results  by  the 
previous  approaches.  The  case  of  multiaffine,  complex  and/real  uncertainties  may  also  be 
handled  using  the  critical  direction  arguments.  However,  for  all  these  cases  the 
computation  of  the  critical  radius  remains  the  major  theoretical  challenge;  fortunately,  for 
a  number  of  uncertainty  descriptions,  such  as  the  ones  considered  in  this  chapter  ,  it  is 
possible  to  obtain  analytical  expressions  for  the  critical  radius. 
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Nomenclature  for  Chapter  3 

s  Laplace  variable 

q  Uncertainty  vector 

n{ja),q)  Uncertain  numerator  polynomial 

d[jco,q)  Uncertain  denominator  polynomial 

p{j(U,q)  Uncertain  characteristic  polynomial 

d^{jco)  Critical  direction  in  the  complex  plane 

k^{o})  Nyquist  robust  stability  margin 

kijj{co)  Multivariable  robust  stability  margin 

Greek  Letters 

a''{co)  Blow  up  factor 

^  Frequency 

Pc  (ft))  Critical  perturbation  radius 

fi{co)  Structured  singular  value 

a{.)  Maximum  singular  value 

p(.)  Spectral  radius 

T(a))  Eigen  template 

^{co)  Critical  eigen  template 


CHAPTER  4 

VARIABLE  STRUCTURE  CONTROL  DESIGN  FOR  REDUCED  CHATTER  IN 

UNCERTAIN  STATE  DELAY  SYSTEMS 

4.1   Introduction 


The  study  of  uncertain  state-delay  and  input-delay  systems  has  received  much 
attention  by  researchers  in  the  past  decade  [cf.  Zaveri  and  Jamshidi,  1987  and  Hale  and 
Lunel,  1993  for  background  and  extensive  references].  One  reason  for  this  interest  is  that 
time  delay  is  encountered  in  various  engineering  systems  and  can  be  the  cause  of 
instability.  Another  reason  is  that  practical  control  systems  unavoidably  include  uncertainty 
or  disturbances  due  to  modeling  or  measurement  errors  and  other  approximations.  Several 
authors  have  investigated  the  problem  of  stabilization  of  uncertain  state-delay  systems.  The 
main  strategies  for  stabilization  of  delay  systems  include  designing  state  or  output  feedback 
controllers  (using  pole  assignment,  Lyapunov  or  LQ  theory),  and  variable  structure 
controllers  (VSC).  Cheres  et  al.,  [1989]  construct  a  min-max  controller  from  the 
knowledge  of  the  upper  bound  on  the  delay.  Shen  et  al,  [1991]  propose  a  memoryless 
linear  state  feedback  based  on  Riccati  equation  approach  [also  see  Wang  and  Lin,  1988]. 
Phoojaruenchanachai  and  Furuta  [1992]  also  construct  a  memoryless  state  feedback,  but 
consider  a  larger  class  of  systems.  Trinh  and  Aldeen,  [1994]  consider  interconnected 
systems  and  propose  memoryless  state  feedback  controllers  to  stabilize  the  uncertain  state- 
delay  systems.  Niculescu  et  al,  [1996]  design  memoryless  state  feedback  controllers  to 
stabilize  an  uncertain  state  delay  system  with  constrained  input  via  ricatti  approach. 

Recently,  sliding  mode  control  has  been  employed  as  a  tool  for  stabilization  of 
uncertain  state-delay  systems.  The  sliding  mode  control  possesses  well-known  features 
that  make  it  very  attractive  for  control  systems.  These  include  fast  response,  insensitivity 
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to  parameter  variations,  and  decoupling  design  procedure,  among  others.  Recently  Shyu 
and  Yan  [1993]  and  Oucheriah  [1995]  have  used  VSC  to  guarantee  the  stability  of  uncertain 
delay  systems,  deriving  sufficient  conditions  that  depend  on  the  size  of  the  delay. 
However,  the  Shyu- Yan  method  suffers  from  the  disadvantage  that  the  control  synthesis 
procedure  involves  unknown  matrices;  hence  the  control  law  can  not  be  implemented. 
Oucheriah  proposes  a  dynamic  switching-surface  control  scheme  based  on  pole- 
assignment.  Both  of  the  above  approaches  involve  strong  discontinuous  control  across  the 
switching  surface  to  overcome  the  effects  of  the  uncertainty,  and  hence  suffer  from  severe 
chattering.  Chattering  is  undesirable  because  it  might  excite  the  unmodeled  high-frequency 
components  of  the  system,  and  it  may  lead  to  premature  wear  or  tear  of  the  actuators.  The 
recent  technique  of  perturbation  compensation  has  been  shown  to  provide  effective  chatter 
suppression  for  delay-free  systems  [Kim,  1992;  Chan,  1996].  This  method  exploits  the 
specific  structure  of  VSC  to  estimate  and  compensate  for  the  effect  of  perturbation,  and 
thus  results  in  reduction  or  removal  of  chattering. 

This  chapter  proposes  a  VSC  design  methodology  to  robustly  stabilize  an  uncertain 
state-delay  system  with  nonlinear  and  possibly  mismatched  uncertainties,  utilizing  the 
concept  of  perturbation  compensation.  Sufficient  robust  stability  conditions  are  derived 
which  are  independent  of  the  size  of  the  delay;  furthermore  no  matching  conditions  on  the 
uncertainty  are  assumed.  The  main  advantages  of  the  proposed  approach  are  that  (i)  a 
realizable  control  law  is  obtained  which  does  not  need  the  bounds  on  the  system 
uncertainty,  (ii)  chattering  reduction/elimination  is  achieved,  (Hi)  delay-independent 
sufficient  robust  stabihty  conditions  which  can  accommodate  uncertainty  in  the  delay  are 
derived,  and  (iv)  guidelines  are  given  for  the  design  of  the  control  matrix  involved  in  the 
specification  of  the  switching  function.  Finally  the  results  are  illustrated  with  an  open-loop 
unstable  example. 
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4.2   Sliding  Mode  Control  Design 

Consider  an  uncertain  state-delay  system  subject  to  uncertainty  and/or  external 
perturbations  governed  by  the  equations 

xit)  =  Axit)  +  A^x{t  -h)  +  Bu{t)  +  f{x{t),t)  +  /^  {xit  -  h),  t) 
x{d)^(p{6)    ,     de[-h,0] 

where  x(t)  e  i?"  is  the  state  vector  which  is  assumed  measurable  and  has  the  initial  state 
^(o)  =  a:o,  u(t)  6  R""  is  the  input  vector;  A  e  i?"'"  and  B  e  R'"""  are  constant  matrices; 

^t)  is  a  continuous  vector-valued  initial  function;  and  h>  0  is  the  time  delay.  The  vector 
functions  f{x{t),t)e  R"    and  f^i(x(t - h),t)e  /?"  represent  nonhnear  perturbations  that 

depend  on  the  current  state  x(t)  and  on  the  delayed  state  x(t-h)  of  the  system, 
respectively.  It  is  assumed  that  the  modeling  uncertainties  satisfy  the  bounds 

\\f{xit),t)\\     <    k  \\x{t)\\ 
\\fa{x(t-h),t)\\<    kj\\xit-h)\\ 

where  k  and  k^j  are  known  positive  real  constants,  and  ||  ■  ||  is  any  vector  norm  or 

corresponding  induced  matrix  norm.  Notice  that  no  matching  assumptions  are  used.  It 

should  be  noted  that  these  bounds  are  required  not  for  controller  design,  but  for  robust 

stability  analysis. 

Two  definitions  are  useful  for  the  ensuing  developments.  The  matrix  measure 

[Vidyasagar,  1993]  is  defined  as  the  function 

ft :  C"'^"  ^  R  :  ^(A)  =  lim  "  " 

£^0+  £ 

The  matrix  measure  is  also  known  as  the  logarithmic  norm.  For  an  excellent  discussion 

on  logarithmic  norms,  see  [Strom,  1975]. 

Also,  a  class  of  matrix  pairs  (A,  B)  is  introduced  for  the  stabilization  problem, 

defined  through  the  set  membership 

U  :=  {(A,fi)  |B  F  G  i?"'''"  :  fi{A  +  BF)  <  o}  (4.3) 
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where  11  clearly  represents  the  class  of  systems  for  which  there  exists  a  state-feedback  law 
M(r)  =  FxiO  such  that  //(A  +  BF)  <  0 . 

The  objective  is  to  design  a  control  law  that  assures  reachability  of  the  sliding  mode 
within  a  finite  time  and  asymptotic  stabiUty  of  the  system  in  the  sliding  mode,  even  in  the 
presence  of  uncertainties  and/or  external  disturbances.  The  design  involves  two  major 
phases:  The  first  phase  is  the  selection  of  a  switching  surface  manifold  that  has  desirable 
dynamics,  and  the  second  phase  is  the  determination  of  an  appropriate  control  law  that 
guarantees  the  existence  and  reachability  of  the  sliding  mode  and  the  effective  rejection  of 
disturbances. 
Switching  Function  and  Control  Design 

This  section  discusses  the  selection  of  an  appropriate  switching  function.  In  the 
sliding  mode  the  system  must  not  only  be  stable  but  also  be  robust  to  any  kind  of  bounded 
perturbations  as  specified  in  (4.3).  The  procedure  is  as  follows.  First,  select  matrices  F 
and  Fj  such  that  matrix  A  =  A- BF  is  a.  stable  matrix,  and  that  matrix  Ad  =A^- BF^ 
satisfies  \Ad\  <  \\A^\\  .  At  this  point  we  anticipate  that  the  role  of  F^  is  to  reduce  the 
effect  of  delay,  since  the  stability  conditions  derived  later  in  this  chapter  involve  |At/||. 

Using  the  matrices  described  above  the  system  equation  (4.1)  can  be  rewritten  in  the  form 

i(0  =  (A  +BF  )xit)  +  {Ad  +BFa  )x{t-h)  +  Bu{t)  + f{x{t),t)  + f^{x{t-h),t){AA) 

An  integral  switching  function  [Fernandez  and  Hedrik,  1987]  is  adopted  in  the  modified 
form  proposed  by  Shyu  and  Yan  [1993]  namely 

t      _  _ 

s{t)  =  Cxit)  -  l{CAx(r)  +  CA^x(T:-h)]dt 

0  (4.5) 

where  matrix  C  e  R'^"  is  a  control  design  matrix  that  is  chosen  such  that  the  product  CB 

is  nonsingular.  The  reason  for  choosing  such  a  switching  function  instead  of  a  traditional 

combination  of  the  states  is  that,  this  choice  results  in  a  desirable  switching  dynamics.  This 

will  be  cleai-  after  the  design  of  the  controller. 
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The  goal  is  to  find  a  VSC  controller  that  drives  the  system  to  the  switching  surface 
and  ensures  a  sliding  mode  condition.  The  controller  structure  adopted  is  of  the  form 

uit)  =  u^^it)  +  u,.(t)+  Up{t) 
which  contains  three  terms.  The  first  term  in  turn  is  of  the  form 

Ueqit)  =  -{Fx{t)  +  Fax{t-h))  (4.6) 

and  thus  represents  the  model-based  equivalent  control  for  the  nominal  system  in  the 
absence  of  perturbations,  obtained  by  setting  s{t)  =  0  [DeCarlo  et  ai,  1988].  The  second 
term 

H,(0    =  -  {CB)-^  {Ks{t)  +d  sgn[s{t)])  (4.7) 

is  the  reaching  control  law  where  the  symmetric  positive-definite  matrix  K  and  the  scalar 
J  >  0  are  design  parameters.  The  final  term  is  intended  to  cancel  the  effect  of  perturbations 
on  the  dynamics  of  the  switching  function  [Chan,  1996]  and  is  given  by 

Up{t)    =  -  (CBfTpit)  (4.8) 

with 

Tp{t)    =i(0  -  CB(u,{t)  +Up{t-T))  (4.9) 

where  Tp(t)  is  designed  to  reject  the  effects  of  the  perturbation  signal.  In  order  to  avoid  an 
algebraic  loop  in  calculation,  a  deliberately  delayed  value  u  Jt  -  T)  is  used.   The  value  of 

the  computational  time  delay  T  is  chosen  to  be  small.  As  long  as  the  variations  in  the 
system  dynamics  are  slower  than  the  computational  delay,  this  will  be  good  approximation. 
However,  it  should  be  admitted  that  regardless  of  the  value  of  T,   Up{t)^Up{t-T). 

However,  as  long  as  the  estimation  error  is  less  than  the  actual  perturbation,  we  only  gain 
by  using  perturbation  compensation.  Also  a  finite  difference  scheme  is  used  to  estimate 
s{t).  This  could  also  be  a  source  of  errors  if  the  measurements  are  noisy. 

The  ideal  VSC  approach  assumes  that  it  is  possible  to  switch  infinitely  fast  across 
the  switching  surface.    The  phenomenon  of  non-ideal  but  fast  switching  is  labeled  as 
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chattering  [DeCarlo  et  al,  1988].  Interestingly,  the  name  comes  from  noise  generated  by 
the  relays.  Chattering  is  undesirable  for  two  reasons:  (i)  It  might  excite  the  unmodeled  high 
frequency  plant  dynamics  which  could  lead  to  unforeseen  dynamics.  Note  that  in  certain 
applications  it  is  not  of  concern  as  the  frequency  of  switching  is  far  above  the  structural 
frequencies  of  the  mechanical  system  [DeCarlo  et  al,  1988].  However,  in  most 
applications  this  will  be  a  major  concern,  (ii)  It  might  lead  to  the  premature  wear  and  tear  of 
the  actuators.  It  should  be  noted  here  that  it  is  the  discontinuous  control,  which  in  turn 
leads  to  chattering,  enables  one  to  guarantee  the  robustness  features  of  sliding  mode 
control.  In  a  sense,  we  gain  robustness  but  pay  for  it  with  chattering. 

To  alleviate  chattering,  many  approaches  have  been  proposed  in  the  literature. 
Here,  we  give  a  brief  review  of  the  past  efforts.  Although,  they  have  not  been  used  for 
sliding  mode  control  of  state-delayed  systems,  they  are  nevertheless  useful  as  a  concept. 
One  of  them  is  the  "boundary  layer  approach"  [see  Slottine  and  Sastry,  1983]  which  uses 
the  continuous  control  law  and  time  varying  sliding  surfaces.  Although  this  succeeds  in 
eliminating  chattering,  it  is  no  longer  possible  to  assure  the  asymptotic  stability  of  the 
system.  Rather,  one  has  to  be  content  with  uniform  ultimate  boundedness  of  the 
trajectories  within  a  neighborhood  of  the  origin  depending  on  the  boundary  layer. 
Ambrosino  et  al,  [1984]  cascaded  a  low-pass  filter  after  the  VSC  to  attenuate  the  high 
frequency  switching  component.  In  this  case  chattering  is  alleviated  at  the  expense  of 
robustness.  Espana  et  al,  [1984]  propose  a  zone  in  the  vicinity  of  the  switching  surface  in 
which  the  feedback  gain  is  adjusted  in  order  to  avoid  chattering.  Morgan  and  Ozguner 
[1985]  propose  to  decrease  the  control  gain  in  order  to  minimize  chattering,  but  the 
disadvantage  is  the  increased  time  that  it  takes  to  achieve  sliding  mode.  Chang  ef  al, 
[1990]  proposed  an  adaptive  alleviation  algorithm  which  requires  recursive  prediction  at 
each  step.  Their  procedure  involves  computing  the  derivative  of  the  states,  but  they  assume 
that  they  are  either  available  or  can  be  computed  without  errors.  Elmali  and  Olgac  [1992] 
do  an  on-Une  estimation  of  the  effect  of  the  perturbations  based  on  the  measurements  of  the 
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state  dynamics.  They  assume  that  the  error  in  estimation  is  bounded  within  a  scalar 
multiple  of  the  estimation.  Also  see  an  experimental  verification  of  their  approach  in  Elmali 
and  Olgac  [1996].  Chan  [1996]  proposed  an  perturbation  compensation  approach  utihzing 
the  specific  structure  of  the  variable  structure  control.  The  approach  seeks  to  estimate  the 
effect  of  perturbations  on  the  switching  dynamics  with  the  assumption  that  the  relative  error 
in  the  estimation  of  the  switching  dynamics  is  bounded  and  known.  One  of  the  major 
advantages  of  the  perturbation  compensation  schemes  is  that  one  need  not  have  a  priori 
bounds  on  the  uncertainties  for  controller  design. 

Note  that  with  the  application  of  the  control  law  (4.6)  -  (4.9)  and  the  choice  of  the 
switching  function  (4.5),  the  derivative  function  s{t)  is  given  by 

s{t)  =  CB  [u,{t)  +Up{t))  +  C{f{xit),t)  +  fMt-h),t)) 

s{t)=   -{Ks{t)+d  sgn[s{t)])  +  [c{f{x{t),t)  +  fa{x{t-h),t))-Tp] 

It  is  clear  that  s{t)  is  driven  by  the  estimation  error,  and  the  purpose  of  control  is  to 

counter  the  effect  of  perturbations  on  the  switching  surface.   The  usual  approach  is  to 

overpower  the  uncertainty  by  discontinuous  control.  For  example,  in  this  case  the  value  of 
d  should  be  chosen  such  that  d>\C\[k  \\xit)\\+  k^i\\xit-h)\\).    In  fact,  a  choice  of  this 

type  is  made  by  [Oucheriah,  1995].  For  proving  the  reaching  condition,  this  is  indeed  a 
clever  selection  of  parameters.  However,  there  are  some  significant  problems  with  this 
choice.  Firstly,  we  are  assuming  it  is  always  the  worst  case  that  is  affecting  the  system. 
Secondly,  there  is  no  method  to  estimate  even  conservatively,  the  norm  of  the  states  when 
the  system  oscillates  near  the  switching  surface.  Hence,  while  we  are  guaranteed  to  suffer 
from  chattering,  we  have  no  available  worst  case  upper  bound. 

In  contrast,  the  effect  of  Up{t)  is  to  ease  the  burden  on  the  discontinuous  control  by 

canceling  the  effect  of  perturbations  on  the  switching  dynamics.  The  discontinuous  control 
is  now  utilized,  to  combat  the  residual  error  in  estimation  rather  than  the  whole 
perturbation.  Irrespective  of  the  actual  value  of  the  perturbation,  as  long  as  the  estimation 
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scheme  works  well,  this  scheme  will  outperform  the  conventional  VSC.  One  point  needs 
to  be  clarified  here.  Estimating  s{t)  is  in  fact,  equivalent  to  saying  that  we  are  estimating 
the  perturbation.  Our  contention  here  is  that  s{t)  can  be  estimated  from  the  measurable 
states  of  the  system  via  (4.5).  It  goes  without  saying  that  we  are  introducing  a  new 
parameter  here,  namely  the  bound  on  the  estimation  error.  Apart  from  the  measurement 
error  in  the  states,  there  is  also  the  danger  of  noise.  Since  s{t)  depends  on  x(t),  it  is 
assumed  that  proper  filtering  strategies  will  be  employed  in  practice.  In  the  worst  case,  one 
has  to  cut  off  the  perturbation  compensation  and  has  to  live  with  the  chattering.  By  worst 
case,  we  mean  that  the  computed  perturbation  signal  is  greater  than  the  maximum  norm  of 
the  allowed  perturbation,  i.e., 

Tp  {t)\    <  \c{f{xit),  t)  +  fMt-h),  0)1  <  \C\(k\x{t)\  +  k^ [c{t  -  h)\) 

However,  this  decision  can  only  be  made  depending  on  the  application  at  hand.  For 
theoretical  development,  in  order  to  prove  the  reaching  condition  within  a  finite  time,  we 
need  to  assume  certain  error  bounds  on  the  derivative  estimation.  We  assume  that  the 
relative  error  in  the  estimation  of  s{t)  is  bounded. 


<5 


with  the  positive  scalar  bound  S  being  known.  An  estimate  of  d  can  be  obtained  via 
simulation  studies  with  appropriately  assumed  perturbations.  Then  with  the  choice 
d{t)  >  &g^f (0  +  £,  e  >  0,  and  K  the  reaching  condition  can  be  satisfied  and  in  addition, 

the  transient  behavior  can  be  shaped. 

Theorem  4.1.  The  uncertain  delay  system  (4.1)  with  model-uncertainty  and 
bounded  measurement  uncertainty  given  by  (4.2)  and  (4.10)  achieves  the  desired  sliding 
mode  within  a  finite  time  if  the  control  law  u(t)  -  u^Jt)  +  u,.{t)+  uJt)  is  applied, 

where  the  individual  components  are  given  by  (4.6  )-(4.9). 
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Furthermore,  the  sliding  mode  is  attained  within  t  =  t^,  where 


t,.< 


f       2       ^ 


In 


V^minWy 

and  the  upper  hound  on  the  norm  of  the  switching  surface  is  given  by 

^min(^) 


KOI  ^ 


Ko)|+ 


V2e 


^m.n(K), 


■42£ 


-^minC^) 


1  T 

Proof:    Choose  a  positive  definite  function    V{t)  =   -  s{t)  s{t).       Then 
differentiating  V{t)  along  the  trajectories  of  the  system  and  using  (4.5)-(4.10)  yields: 


V{t)  =  s{tYs{t)^s{t)' 


-Ks{t)  -d^  +  {c{f{x(t),t)  +  Mxit-h),t))-Tp{t) 


(4.11) 


In  the  rest  of  the  proof  the  time  argument  is  dropped  for  convenience  of  notation.  Noting 
that    As=   C[f{xit),t)  +  f^{x{t-h),t))-Tp{t)  and  using  the  bound  on  the  estimation 

error  (4.10)  the  equation  (4.1 1)  can  be  used  to  derive  a  bound  on  V{t)  as  follows: 

V{t)=  -  s^ Ks  -  ds'^ sign{s)  +  s'^ As 
We  note  that  since  s{t)  is  a  vector, s  sign{s)  =  \s\^  >  \s\2.  Hence, 

V{t)<  -s^Ks-d\s\  +\s\\As\ 
Since  d  =  (5|sg^;(0|  +  e,  and  |zls|  <  5|sg^,(f)| 

Vit)<  {-^^.^(K))V-s^ 
The  above  is  a  scalar  differential  equation  and  the  solution  is  given  by 

1         ^min(*^). 


no 


< 


no)+ 


^min(^) 


^™n(^) 
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<t)\ 


< 


[s(0)|  +  .   ^^ 


1        ^min(-^).  f^ 

— ^ — '        V2e 


^min(^) 


Note  that  first  term  is  always  positive  and  decreasing,  and  the  second  term  is  a  positive 
constant  and  hence,  their  difference  will  attain  zero  within  a  finite  time.  It  is 
straightforward  to  compute  the  upper  bound  on  the  reaching  time.  A  detailed  proof  is  given 
in  Appendix  F.  Thus  the  system  is  guaranteed  to  reach  the  sliding  mode  within  a  finite  time 
since  the  reaching  condition  s{t)^s{t)  <  0  for  s{t)  ^   0  Is  satisfied.  V 

Remarks: 

1 .  The  rate  at  which  the  switching  surface  is  reached  is  determined  by  the 
controller  parameters  K  and  d .  Therefore  by  adjusting  these  parameters  the 
designer  may  select  an  appropriate  reaching-transient  behavior. 

2.  Chattering  in  the  input  variable  is  unavoidable  for  any  positive  value  of  d, 
because  u^{t)  ^  ±d  when  the  system  oscillates  around  s  =  0.  To  reduce 
chattering,  the  parameter  d  must  be  chosen  as  small  as  possible.  In  some  cases 
s{t)  is  measured  directly,  therefore  d  can  be  set  equal  to  zero  without  any  loss. 
In  other  cases  igstCO  can  be  compared  with  simulated  s{t)  (for  a  given  set  of 
perturbations)  to  aid  in  choosing  a  suitable  the  parameter  d .  Alternatively,  one 
can  choose  to  set  (i  =  0  and  specify  a  large  K  to  overpower  the  estimation 
errors  if  any  are  present.  However  in  the  last  case  the  penalty  to  be  paid  is  high 
initial  controller  gains. 

4.3   Asymptotic  Stability  with  Perturbation  Compensation 

Once  in  the  sliding  mode  the  system  equations  reduce  to 

i(0  =  Ax(t)  +  Adx(t  -h)  +  f{xit),  t)  +  fj  {xit  -  h),  t)  -  B{CB)^  Tp  (0  (4_  j  2) 
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which  is  obtained  from  (4.4)  after  substituting  (4.6)-(4.8).  The  intent  now  is  to  derive 
conditions  that  establish  the  asymptotic  stability  of  system  (4.12).  First,  we  note  that  if 
perturbation  compensation  is  not  employed,  then  (4.12)  reduces  to 

x{t)  =  Ax(t)  +  Adx(t-h)  +  f(x{t),t)  +  f^i{xiT-h),t)  ^^  ^^^ 

With  perfect  perturbation  compensation,  the  system  equations  become 

i(0  =  Ax(t)  +  Adx{t  -h)+  y(f{x{t),  t)  +  /^  {xit  -  h),  0)  (4.14) 

where  y^ll-BiCSy^C 

There  are  basically  two  approaches  in  the  literature  for  deriving  conditions  for  the 
asymptotic  stability  of  these  systems.  One  is  to  follow  Lyapunov-type  of  arguments  [Hale 
and  Lunel,  1993],  and  the  other  is  to  use  norm-based  inequaUties  [Mori  et  a/.,  1981].  Both 
approaches  result  in  sufficient-only  stability  conditions. 

Lemma  4.1.     The  delay  system  (4.13)  with  nonlinear  uncertainties  (4.2)  is 
asymptotically  stable  if  the  following  condition  is  satisfied 

ld[A)  +  \\A^\\  +  k  +  k^i<0  (4.15) 

Proof.    See  Appendix  B.  V 

Theorem  4.2    The  delay  system  (4.14)  with  nonlinear  uncertainties  (4.2)  is 
asymptotically  stable  if  the  following  condition  is  satisfied 

lj{A)  +  \\A^\\  +  k\\'}\\  +  kM<0  (4.16) 

Proof.    The  proof  is  idendcal  to  the  proof  of  the  previous  lemma.  V 
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Remarks 


It  is  possible  that  different  choices  of  a  given  norm  and  the  corresponding 
matrix  measure  give  different  conclusions  about  stability.  This  is  consistent 
with  the  fact  that  the  conditions  derived  are  only  sufficient. 
In  order  to  improve  the  sharpness  of  the  stability  condition,  a  similarity 
transformation  z{t)  =  T~  x{t)  can  be  introduced  so  that  the  system  equation 
(4.14)  becomes 

z{t)  =  T-^ATz{t)  +  T-^AdTz{.t-h)+i{f{x{t),t)  +  fMt-h),t))    (4.17) 
where  7    =  T~^ [I  -  B{CB)  ' C ) .  In  this  case  the  stability  condition  reads 


^(r"Ur)+ 


r"'A^r 


^Kmii  +m 


<o 


(4.18) 


The  usual  approach  is  to  find  a  transformation  matrix  T  to  diagonalize 

4.   Under  perfect  estimation  and  the  assumption  of  matched  uncertainties  [Shyu 
and  Yan,  1993]  the  perturbation  error  is  (f  +  f^)  + BUp  =  0.  The  stability 


condition  then  reduces  to /i(r    AT 


+ 


T~^AjT 


<   0.   Of  course,  there  is  a 


5. 


penalty  for  imperfect  estimation  and  this  is  reflected  in  the  reduction  of  the  size 
of  the  perturbations  the  system  can  withstand  and  still  remain  stable. 
The  conditions  (4.15)  and  (4.16)  are  delay  independent.  This  could  be 
advantageous  in  cases  where  delay-dependent  stability  conditions  cannot  assure 
stability.  The  example  in  the  chapter  illustrates  this  fact.  However,  since  all 
conditions  reported  to  date  in  the  literature  (including  this  one)  are  sufficient 
only,  no  general  claims  can  be  made.  It  is  also  possible  to  derive  delay- 
dependent  conditions  for  the  problem  considered  in  the  chapter. 
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Optimal  Choice  of  C 

The  stability  condition  (4.18)  and  the  condition  derived  in  Shyu  and  Yan  [1993] 
involves  the  matrix      /'    =   T~^I -BiCBJ^c).    Then  the  matrix  C  is  chosen    to 


minimize 


T-^(l  -  B{CB)~^C) 


.    There  appears    to  be  no  documented  systematic 

F 


method  for  optimizing  the  left  hand  side  of  (4.18)  (or  similar  stability  conditions  for  other 
approaches)  over  all  allowable  matrices  T  and  C.  To  get  the  tightest  possible  stability 
condition  for  the  result  in  Shyu  and  Yan  [1993]  the  matrix  C  must  be  chosen  to  minimize 
the  norm  of  j .  We  propose  a  new  systematic  approach  to  design  matrix  C. 

The  following  definition  and  inequalities  concerning  Frobenius  norm  are  useful. 

IIf' 

value,  respectively.  Then  [Golub  and  VanLoan,  1989] 


Let  A   E  9t'"'",  and  let  ||  L,  and  af  represent  the  Frobenius  norm  and  'i'th  singular 


Jm    n  2 

II 


The  following  relationships  are  useful 

T,    +  •••  +C7„;    p   =   min|m,n| 

(4.19) 


'f     =    (Ji    +  ■■■  +(7„;    p   =   min{m,n} 


Ii2     -    \l"-\\F    -     ^'MK^II2 


Lemma    4.2.    Let  N  e  9t'"'"  and  M  e  9t'"".    Then  the  best  approximate 
solution  to  the  problem 


min     \\  M  -  NX 


is  Zopt    -  {n'^n]    N^M. 


Proof:    See  Barnett  [1971]. 

The  definition  of  X^^^  implies  that  for  all  Z,  X   g  ?K' 


either 


M  -NX  I    > 


M  -  NX,^, 
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or 


M  -  NX 


M  -  NX, 


opt 


and  \\X\\  > 


X. 


opt 


It  is  instructive  to  note  that  (n^n)    n'^  is  ttie  Moore-Penrose  generalized  inverse  of 
matrix  A^. 


rmn 

C  E  9t""" 


Theorem  4.3.    The  best  approximate  choice  of  matrix  C  for  the  following  problem 

T~^(l  -  B{CBy^c) 
B  T     T     1,  where  a  is  a  scalar. 

Proof.    Define  the  auxiliary  matrix  Z  =  {CB)~  C.  The  minimization  problem 

reduces  to 

T~\l  -  BX) 


mm 

X  e  Si""" 


-1 


Substituting  M  =   T     and  N  =  T    B  in  Lemma  4.2  the  best  approximate  solution  is 
found  to  be 

z,p,  =  [b^t-^t-'b^b'^t-^t-'  =  (Coptfi)~'c  ^^ 

It  is  obvious  that  Qpt  =  B^T~^T~^ .  Note  that  multiplying  by  a  scalar  a  does  not  change 


the  value  of  the  functional  being  minimized,  hence  Qpt  =  ayB^T     T    j 
Remarks: 


V 


1 .  Note  that  Theorem  4.3  gives  the  optimal  solution  in  the  sense  of  the  Frobenius 

norm.  Using  (4.19),  it  is  easy  to  see  that  for  the  case  where  B   e  9?"^  ,  the 
choice  Copt  -  oc[b^T~^T~^\  is  optimal  in  the  sense  of  the  2-norm  as  well. 

2.  The  Matrix  C  enters  in  the  last  three  terms  of  the  stability  condition  (4.18).  An 
alternative  method  for  finding  the  optimal  matrix  C  is  to  carry  out  a  numerical 
optimization.  However,  noting  that  the  contribution  from  the  last  term  will  be 
small  for  a  well  performing  estimation  strategy,  the  last  term  can  be  ignored  for 
analysis  purposes.  Then,  from  triangular  inequality  it  is  easy  to  verify  that  the 
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result  of  Theorem  4.3  provides  a  sub  optiraal  upper  bound  in  the  sense  of  the 
Frobenius  norm. 
Implementation   Issues 

1 .  The  estimation  error  depends  on  the  approximation  scheme  used.  Although 
there  is  no  a  priori  information  on  the  estimation  error,  different  estimation 
strategies  can  be  compared  by  assuming  appropriate  perturbations  and  tracking 
the  estimation  error. 

2.  Also,  note  that  the  main  challenge  in  the  proposed  approach  is  that  of 
approximating  s{t)  by  s^^A^).  In  general,  a  finite  difference  estimate  is  used. 

However,  this  may  cause  a  problem  near  the  switching  surface  where  some 
chattering  combined  with  a  small  step  length  may  contribute  to  large  estimation 
errors.  Also,  there  is  the  unavoidable  nuisance  of  noise,  and  the  danger  of 
differentiating  it.  Hence  it  is  necessary  to  employ  filtering.  Although  the 
switching  dynamics  with  perturbation  compensation  is  complex,  one  can  still 
perform  a  low  pass  filtering  on  s{t)  for  the  purpose  of  estimating  Sg^({t)  by 

finite  difference  [Elmali  and  Olgac,  1996].  Note  that  the  actual  value  of  s{t)  can 
be  used  elsewhere  in  the  controller.  Hence,  this  does  not  destroy  the  guarantee 
of  asymptotic  stability.  The  cutoff  frequency  for  the  filter  can  be  chosen  based 
on  the  maximum  frequency  component  that  we  wish  to  process  in  the  state 
dynamics.  For  example,  if  temperature  in  a  reactor  is  the  state  variable  that  is 
used  in  the  construction  of  the  switching  function,  a  cut-off  frequency  of  100 
Hz  is  acceptable,  whereas  if  pressure  in  a  flowing  line  is  the  variable  then  one 
needs  to  go  for  higher  frequencies.  With  the  advances  in  digital  signal 
processing  these  decisions  can  be  made  scientifically  depending  on  the 
application  at  hand. 
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4.4   Example 


Consider  the  uncertain  time-delay  system  studied  by  Shyu  and  Yan  [1993] 


• 

"  1 

1.5" 

"0 

-1" 

X    = 

x(r)  + 

0.3 

-2_ 

0 

0_ 

x{t-h)  + 


10" 
1 


u  +  fix(t),t)  +  fjix(t - h),t),    t>0 


The  delay  is  /z  =2,  and  the  unknown  non-linear  perturbations  are  of  the  form 

f{x{t),t)  =  0.2x{t)sin[x{t))  and  f^i{x{t - h),t)  =  OAx{t - h)sm[x{t - h)) 

The  design  objective  is  to  asymptotically  stabilize  the  uncertain  system.    The  initial 
condition  is  specified  as  jc:(r)  =  [1,3]     -h  <   t  <  0.  Note  that  for  the  delay  considered, 

the  sufficient  stabihty  condition  of  Shyu  and  Yan  [1993]  is  violated. 

Also  note  that  the  homogeneous  part  of  the  system  is  unstable  because  A  has  a 
positive  eigenvalue  A  =  1.1432.  Choosing  the  same  control  matrices  investigated  in 
Shyu  and  Yan  (1993),  namely  F  =  [0.3    0.15]  ,  F^=  [O    -O.l]  leads  to 


A    = 


-2 
0 


0 

2.15 


0     0  ■ 

0    -0.1 


The  particular  choice  f  adopted  makes  A  diagonal  (and  hence  T  =  I),    and  the 

In  this 


choice  of  i^/  ensures   A^    <  |Aj||.   As  per  Theorem  4.3,  C^^^    =   a  [B  T~ 

design  we  choose  K  =  2,  and  d  =  0.  Then  choosing  a  =  0.1    leads  to  C  =  [1    0.1].    If  a 
positive  value  for  d  is  chosen,  it  will  result  in  some  chattering. 

From  (4.16)  we  can  find  the  maximum  tolerable  estimation  error.  For  the  1-norm 
and  infinity  norm  we  find  that  asymptotic  stability  is  assured  if  e  <  1.4445 ,  while  for  the 
2-norm  asymptotic  stability  is  assured  if  £<  1.6085.  The  delay  introduced  in  the 
computation  of  m„(0  is  0.01  seconds. 

The  simulation  results  are  presented  in  Figures  4.1  -  4.4.  Figure  4.1  shows  the 
evolution  of  state  norm  as  a  function  of  time,  reflecting  the  asymptotic  stability  of  the 
system.  The  switching  function  and  its  evolution  towards  the  sliding  mode  is  shown  in 
Figure  4.2.  The  shding  mode  is  reached  rather  quickly,  and  most  importantly,  a  chattering- 
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free  behavior  is  observed.  Figure  4.3  shows  that  the  input  variable  is  also  smooth  and  free 
of  chattering.  The  jump  in  the  input  around  t  =  2  reflects  the  fact  that  control  variable  is 
responding  to  the  delayed  states  coming  into  effect.  Figure  4.4  shows  the  evolution  of  the 
estimation  error  as  a  function  of  time,  which  also  reaches  the  value  of  zero  asymptotically. 
From  Figure  4.4,  it  is  obvious  that  the  estimation  error  in  the  sliding  mode  is  much  less 
than  1.6085,  thus  implying  robust  stability  via  equation  (4.16).  In  this  example,  there  is 
no  chattering,  and  the  system  is  asymptotically  stable. 

4.5   Conclusions 

This  chapter  considered  the  robust  stabilization  of  uncertain  time-delay  systems 
with  mismatched  nonlinear  uncertainties.  A  practically  implementable  control  design 
methodology  is  proposed  and  an  illustrative  design  example  is  given.  The  present  method 
can  significantly  reduce  or  completely  eliminate  chattering  as  compared  to  previous 
approaches. 
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Nomenclature  for  Chapter  4 


A 

State  matrix 

A, 

State  delay  matrix 

B 

Control  matrix 

C 

Switching  function  matrix 

d 

Magnitude  of  discontinuous  control 

f 

State  dependent  nonlinear  perturbation 

fd 

Delayed  state  dependent  nonlineai"  perturbation 

F 

State  feedback  matrix 

Fd 

Delayed  state  feedback  matrix 

h 

delay 

K 

Positive  definite  matrix 

k 

Bound  on  state  dependent  uncertainty 

Ki 

Bound  on  state  delay  dependent  uncertainty 

s{t) 

Switching  function 

T 

Transformation  matrix 

Tp 

Estimated  perturbation  signal 

u{t) 

Control  input 

x{t) 

State 

z{t) 

Integral  equation  with  deviating  argument 

Greek  Letters 

6  Bound  on  the  derivative  estimation  error 

£  Derivative  estimation  error  in  the  sHding  mode 

£  Bound  on  the  derivative  estimation  eiTor  in  the  sliding  mode 

a  Scalar  multiplier 

0  time  index 
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^0)     Initial  condition 
jj.{A)    Matrix  Measure  of  A 


CHAPTER  5 
SLIDING  MODE  CONTROL  FOR  UNCERTAIN  INPUT  DELAY  SYSTEMS 

5.1   Introduction 

Time  delay  exists  in  various  branches  of  engineering.  From  biological,  optical, 
electric  networks,  to  chemical  reaction  systems,  time  delay  occurs  and  affects  the  stability 
and  performance  of  a  wide  variety  of  systems  [See  Schell  and  Ross,  1986,  Inamdar  et  al, 
1991].  Input  delays  occur  as  one  of  the  common  sources  of  time  delay,  such  as  in 
chemical  processes,  transportation  lags  and  measurement  delays,  etc.  It  is  well  known  that 
the  existence  of  time  delay  degrades  the  control  performances  and  makes  the  closed-loop 
stabilization  difficult.  Hence,  the  study  of  time  delay  systems  has  been  paid  considerable 
attention  over  the  past  years.  The  description  of  time-delay  systems  lead  to  differential- 
difference  equations,  which  require  the  past  values  of  the  system  variables.  There  are  no 
known  methods  to  get  necessary  and  sufficient  conditions  for  the  robust  stability  of  such 
systems.  There  are  many  approaches  that  yield  sufficient  conditions  with  varying  degrees 
of  sophistication. 

The  case  of  perfectly  known  {i.e.,  no  uncertainty)  time-delay  system  is  well  studied. 
It  is  known  that  a  stabilizable  time-delay  system  [see  Kamen  et  al.,  1985,  Logemann, 
1986]  can  always  be  stabilized  by  a  finite  dimensional  controller.  The  study  of  state 
uncertain  delay  systems  has  been  paid  much  attention.  However,  with  regard  to  the  control 
of  uncertain  input  delay  systems,  only  a  few  studies  have  been  reported  in  the  literature. 
Based  on  Riccati-equation  approach,  a  robust  controller  is  derived  in  Kojima  and  Ishjiama 
[1995].  Yan  et  al,  [1997]  propose  a  functional  observer  and  state  feedback  mechanism 
using  the  factorization  approach.  For  a  LQG/LTR  method  see  Lee  et  al,  [1988].  Wu  and 
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Chou  [1996]  propose  a  control  algorithm  from  input-output  feedback  linearization  by 
means  of  a  parametrized  coordinate  transformation.  Recently,  an  LMI  based  approach  was 
given  in  Niculescu  et  al,  [1997]  for  delay-dependent  closed  loop  stability  of  input-delay 
systems.  However,  the  authors  allowed  uncertainty  only  in  the  input  delay. 

One  of  the  common  methods  used  in  chemical  engineering  applications,  for  processes 
with  time-delay,  is  to  use  the  idea  of  the  Smith-predictor  scheme  to  cancel  the  effect  of  time 
delay  [Stephanopoulos,  1984].  However,  if  the  model  is  imperfect  (uncertain)  then  Smith- 
predictor  is  known  to  give  poor  performance.  Some  approaches  have  been  developed  for 
time-delay  compensation  based  on  prediction  strategies  [See  Henson  and  Seborg,  1994]. 

In  this  chapter  we  propose  an  approach  based  on  the  Sliding  Mode  Control  [Hung  et 
al,  1993;  Utkin,  1977,  Decarlo  et  al,  1988].  The  sliding  mode  control  approach 
possesses  many  advantages,  e.g.,  fast  response,  good  transient  performances,  robustness 
to  the  variation  in  plant  parameters  and  external  disturbances.  The  authors  in  [Shyu  and 
Yan,  1993;  Oucheriah,  1995;  Luo  and  Sen,  1993,  Basker  et  al,  1997]  study  the  sliding 
mode  control  for  uncertain  systems  with  state  delay.  Also  see  Luo  et  al,  [1997]  for  a 
sliding  mode  approach  for  uncertain  time  delay  systems  with  internal  and  external  point 
delays  via  various  types  of  feedback. 

This  chapter  proposes  a  SMC  controller  to  robustly  stabilize  an  uncertain  input-delay 
system  with  linear  and  possibly  mismatching  uncertainties.  Sufficient  robust  stability 
conditions  are  derived.  No  matching  conditions  on  uncertainties  ai-e  assumed. 

The  rest  of  the  chapter  is  organized  as  follows:  In  Section  5.2,  a  transformation  is 
adapted  to  convert  the  original  system  with  bounded  input  delay  into  one  without  any 
delay.  The  stabilization  of  the  modified  system  implies  the  stability  of  the  original  system 
as  long  as  the  control  law  is  bounded  [see  Kwon  and  Kim,  1980;  Artstein,  1982].  Section 
5.3  gives  the  guideline  to  design  the  sliding  mode  controller  which  consists  of  the 
switching  function  and  the  control  law.  In  Section  5.4,  the  analysis  of  the  perturbed 
system  is  carried  out,  and  the  reaching  condition  and  the  asymptotic  stability  of  the  sliding 
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mode  are  proven.  In  Section  5.5  some  open  issues  in  the  sliding  mode  literature  are  raised 
and  answers  are  provided.  Section  5.6  considers  constraints  on  the  controller  input. 
Section  5.7  gives  an  illustrative  example  with  an  open  loop  unstable  plant.  Conclusions 
and  future  work  are  discussed  in  Section  5.8. 

5.2   Preliminaries  and  Problem  Formulation 

Consider  the  uncertain  input  delay  system 

x{t)  =  Ax{t)  +  Bu{t)  +  B^u{t-h)  +  f{x,t)  /^  ^N 

where  jc(0  g  jR"  is  the  state,  «(•)  e  J?'"  is  the  control,  A,  B  and  B^  are  constant  matrices 
with  appropriate  dimensions,  h  isa.  known  constant  time  delay.  The  total  plant  uncertainty 
is  bounded  by 

\fix,t)\<P\x{t)\  (52) 

where  |.|  denotes  a  vector-norm  and  ||H  denotes  the  associated  induced  matrix-norm. 
Without  loss  of  generahty,  the  vector  2-norm  ||-||2  is  used.  The  positive  scalar  bound  ^3  is 

assumed  to  be  known.  In  rest  of  the  development,  without  loss  of  generality,  we  assume 
that  m<n.  Note  that  we  do  not  assume  any  matching  condition  on  the  uncertainty.  The 
method  described  here  can  be  trivially  extended  to  the  case  of  multiple  known  delays. 
Assumption  1:  The  pair  \A,  J?]  is  controllable  and  all  the  states  are  measurable. 

Transformation  into  a  delay  free  system 

We  consider  a  transformation  to  convert  the  original  system  (5.2)  into  an  equivalent 
system  without  delay.  This  reduction  enables  the  use  of  known  stabilizing  methods  for 
linear  systems  without  input  delay.  See  Artstein  [1982],  for  extensive  theoretical 
discussions  on  the  vahdity  and  the  appHcations  of  this  type  of  reduction  and  also  see 
Fiagbedzi  and  Pearson  [1986],  for  further  generalizations. 
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Define  the  transfonnation 

zit)  =  xit)  +  \l_j,e'^^'-''~'^B,iuir)dT  (5.3) 

Assumption  2.  The  solution  zit)  is  not  affected  by  a  change  in  u{t)  on  a  set  oj  zero 
Lebesgue  measure  [Artstein,  1982]. 

For  an  arbitrary  bounded  control  uit),  the  system  (5.2)  is  transformed  to  (see 
Appendix  D  for  details) 

z{t)  =  Azit)  +  Bu(t)  +  f(x,t)  (5.4) 

Where  B  =  B  +  e~^''B^i .  It  is  assumed  that  B  has  full  rank  m .  Then,  via  PBH  rank  test 

[Kailath,  1980]  it  is  trivial  to  show  that  the  system  (5.4)  is  also  controllable.    Note  that  the 
system  described  by  equation  (5.4)  is  an  uncertain  linear  system  without  input  delay. 

5.3    The  Design  of  the  Sliding  Mode  Controller 

In  the  absence  of  perturbations,  from  (5.4),  the  nominal  system  is  given  by 

z{t)  =  Azit)  +  Buit)  (5.5) 

In  the  following  we  will  design  the  switching  function  parameters  on  the  basis  of  system 
(5.5).  The  design  consists  of  two  phases.  The  first  phase  is  the  construction  of  the 
switching  function,  so  that,  the  nominal  system  restricted  to  the  switching  surface  is 
asymptotically  stable.  The  second  phase  involves  the  development  of  a  control  law,  which 
must  satisfy  a  sufficient  condition  for  the  existence  and  reachability  of  a  sliding  mode 
within  a  finite  time  period. 
Design  of  the  switching  function 

Since  system  in  (5.5)  is  a  typical  Unear  system,  choose  the  switching  function  to  be 

sit)  =  Cz(t)  (5.6) 

where  C  is  an  m  x  n  matrix  to  be  determined.  The  idea  is  to  design  the  switching  function 
so  that  the  nominal  system  is  asymptotically  stable  in  the  sliding  mode  [Woodham  and 


73 


Zinober,  1993].  Note  that  other  types  of  switching  function  such  as  an  integral  switching 
function  can  also  be  used  to  facilitate  desirable  switching  dynamics.  Since  the  matrix  B 
has  full  rank  m,  it  is  possible  to  reduce  (n-m)  of  the  rows  to  zero,  i.e.,  there  exists  a 
nonsingular  matrix  T  such  that 


TB  = 


0 


(5.7) 


where  Bj  is  an  m  x  m  nonsingular  matrix.  In  fact,  we  remark  that  such  a  T  has  a  simple 
form  for  m  =  1  and  is  given  by 


T  = 


0       - 


yB{n)j 


fn 


1 

0 


B{n-\) 
B{n) 


The  transformed  states  are  given  by 


and  (5.5)  becomes 


w{t)  =  Tz{t) 


(5.8) 


w{t)  =  TAT~^wit)  +  TBu{t) 
while  the  switching  function  becomes 


(5.9) 


s{t)  =  CT-'w{t)  =  Qw{t) 
We  introduce  the  following  partitions  for  later  use. 


H'(0  = 


w{t)' 

W[{t) 

(1) 


[n  -  m) 
(m) 


Q-    [Q     Qi]im) 

{n-m)     (m) 


T=[Tj     T2]n         TAT~^  = 

{n-m)    {m) 


All      A  2 
{n-m)     (m) 


[n  —  m) 
(m) 


(5.10) 


(5.11) 


equation  (5.9)  can  be  rewritten  in  the  form 
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M'(r)  =  Aiiw(0  +  Ai2H'/(0  (5.11a) 

Wi{t)  =  A2iw(t)  +  A22Wi(t)  +  Bjuit)  (5.11b) 

Furthermore,  (5.10)  can  be  written  in  the  form 

sit)  =  Qwit)  +  Q,Wjit)  (5.12) 

where   Q    is   an    mx{n-m)    matrix,   and    Qi    is    an    mxm  nonsingular  matrix. 
Differentiating  (5.12)  and  substituting  (5.1 1)  in  to  the  dynamic  equation  yields 

Ht)  =  (Au  -A,2Qr'Q)Ht)  +  AnQr's{0  (5A3) 

i(0=[(eA„ +!2/A2i)-(2Ai2 +Q/A22)Gr'e]iv(o+(eAi2 +2/^22  )er"s(0+!2/^^^ 

(5.14) 
When  s{t)^0,  equations  (5.13)  and  (5.14)  describe  the  reaching  mode.  On  the  other 
hand,  when  the  system  is  in  sliding  mode,  i.e.,  s{t)  -  0  and  s{t)  =  0,  (5.13)  yields 

^t)  =  [An-A,2Qr'Q)w{t)  (5.15) 

Wi{t)^-Qr^Qw{t)  (5j5a) 

Since  the  pair  (Aj  i ,  A12)  is  completely  controllable  (see  the  Appendix  E  for  a  proof),  there 
exists  an  m  X  (n  -  m)  matrix  K  that  satisfies 

X{A,,-A,2K)  =  r  (516) 

where  A''  represents  the  desired  poles  in  the  sliding  mode.  Since  Ajj  and  A]  2  are  known, 
a  set  of  linear  equations  can  be  solved  to  yield  K.  Comparing  (5.15)  and  (5. 16),  we  have 

Q^QiK  (5.17) 

Choosing  Qj  =  I  yields 

e=[e  Qi]=iK  I]  (5.18) 
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Substituting  (5.18)  into  (5.10)  results  in 

C  =  QT  =  lK  /,„]r  (5.19) 

Therefore,  the  switching  function  has  the  form 

sit)  =  Czit)  =  [K    I]Tz(t)  (5.20) 

It  is  instructive  to  recast  the  dynamics  of  the  system  in  terms  of  our  choice  for  C,  i.e.,  K. 
By  using  (5. 1 1),  it  is  easy  to  show  the  following: 

W(t)  =  {A]^-Ai2K)w{t)  +  Ai2s{t)  (521) 

s{t)  =  [k{A^  1  -  Ai2^)  +  (A21  -  A22K)]w{t)  +  {KA^2+  ^22  >(0  +  '^/"(O  (5.22) 

Wi(,t)  =  s(t)-Kw(t)  (523) 

We  also  remark  that  in  the  sliding  mode,  Wj{t)  =  -Kw{t)  and  hence,  the  asymptotic 
stability  of  (5.15)  ensures  that  lim  Wj{t)  -^  0. 

Remarks: 

(1)  The  above  also  illustrates  an  interesting  system-theoretic  point.  We  are  seeking 
to  control  an  n-dimensional  system  via  a  m-dimensional  input.  From  (5.21)- 
(5.23),  it  is  obvious  that  the  input  exphcitly  affects  only  a  m-dimensional  subset 
of  this  reahzation  and  the  remaining  (n-m)  dimensional  subsystem  is  affected 
via  the  feedback  of  the  states  as  determined  by  the  system.  This  information  is 
not  new.  In  fact,  by  recasting  the  system  in  the  controllability  realization,  this 
fact  can  be  trivially  seen.  But,  in  contrast  to  the  controllability  canonical  form, 
(5.21)-(5.23)  also  guide  in  choosing  a  controller  to  exploit  the  fact  that  the 
matrix  Bj  is  of  rank  'm'.  More  about  this  will  follow  in  later  sections. 

(2)  Note  that  with  the  choice  of  C,  (C5)"^  =  Bf^  and  since  B;  is  nonsingular, 
{CB)~   always  exists. 
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Design  of  the  Control  Law 

The  control  law  chosen  is  a  combination  of  the  traditional  state  feedback  and  a  time 
varying  discontinuous  control  element 

u{t)  =  -(CBy^\CAz{t)  +  k\sit)f&gn[sit)]  +  J(r)sgn[s(r)]]  (5.24) 

where  the  parameters  are  chosen  such  that  k>0,  0  <  a  <  1.0  and  d{t)  =  l5\C\\x{t)\  +  £, 
with  e  >  0.  When  s(t)  is  small,  a  faster  approach  rate  than  the  linear  approach  algorithm 
is  realized.  The  state  feedback  is  designed  to  cancel  the  effect  of  the  state  matrix  A  in  the 
switching  dynamics.  Once  the  system  is  in  the  vicinity  of  the  switching  surface,  the 
magnitude  of  the  discontinuous  control  is  linearly  proportional  to  the  norm  of  the  states. 
Since  the  states  x{t)  and  z{t)  are  available,  realizing  the  above  control  law  is  straight 
forward.  Using  the  partitions  (5. 1 1),  it  is  also  easy  to  verify  that  (5.24)  is  equivalent  to 

u  =  -{Biy\[K{Au-Ai2K)  +  {A2i-A22K)]w{t)  +  {KA^2+^22Ht)} 

+ 

-(fi;f'[/:|s(Orsgn[s(0]  +  ^Wsgn[s(0]]  (5.25) 

Comparing  (5.25)  and  (5.22),  it  is  obvious  that  we  exploit  the  fact  Bj  is  non-singular 
mxm  matrix,  and  obtain  u  by  equating  (5.22)  to  zero. 

With  the  chosen  controller,  the  switching  dynamics  (5.22)  is  of  the  form 

s{t)  =  -k\sit)fsgn[s{t)]  -  d{t)sgn[s{t)]  +  Cf{x{t),t)  (526) 

Remarks: 

(i)  From  (5.26),  it  is  obvious  that  the  discontinuous  control  is  used  to  overpower 
the  effect  of  uncertainty.  Since  the  only  information  we  have  about  the 
uncertainty  is  the  norm  bound,  we  can  not  use  any  lower  value  of  d{t)  and  be 
guaranteed  to  reach  the  shding  mode  within  a  finite  time.  From  (5.25),  it  is 
clear  that  there  will  always  be  chattering  for  this  choice  of  parameters. 
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(ii)     For  the  nominal  system  d{t)  can  be  chosen  arbitrarily  close  to  zero. 

(Hi)  The  parameter  a  can  be  used  in  a  time  varying  fashion  to  guarantee  a  faster  rate 
of  approach  to  the  switching  surface  than  the  normal  choice  a  =  1.0.  For 
example,  a  =  1.0  for  |s(r)|>1.0  and  a  =  0.5  for  |s(r)|  =  LO  would  assure  a 

faster  rate  of  approach  than  the  proportion  control  (i.e.,  a  =  1.0,  Vf  >  0  ). 
5.4   Analysis  of  the  Perturbed  System 

In  this  section,  the  robust  stability  of  the  uncertain  system  (5.1)  is  analyzed.  First,  the 
analysis  is  carried  out  for  system  (5.4),  in  order  to  give  conditions  for  reaching  the  sliding 
mode  and  the  asymptotic  stability  in  the  sliding  mode.  Since  system  (5.4)  is  related  to 
system  (5.1)  via  (5.3),  conditions  are  given  under  which  robust  stability  can  be  inferred  for 
the  original  system  (5.1). 

For  an  uncertain  system  to  be  robustly  stable  under  sliding  mode  control,  the 
following  conditions  have  to  be  met:  (/)  The  system  should  attain  the  sliding  mode  within 
a  finite  time;  that  is,  the  state  trajectories  should  reach  the  switching  surface  and  should  stay 
on  the  switching  surface  for  all  time  thereafter;  and  (//)  the  system  in  the  sliding  mode 
should  be  asymptotically  stable. 
Reaching  Condition  for  Sliding  Mode. 

Theorem  5.1:  The  uncertain  delay  system  (5.4)  with  model-uncertainty  and 
bounded  measurement  uncertainty  given  by  (5.2)  achieves  the  desired  sliding  mode  within 
a  finite  time  if  the  control  law  (5.25)  is  applied . 

Furthermore,  it  can  be  shown  that  the  value  of  a  chosen  1,  then  the  sliding  mode  is 
attained  within  t  =  t,.,  where 


(,<l^lln 


1+* 

£ 


^i7(oMoy^ 


V 


2 
and  the  upper  bound  on  the  norm  of  the  switching  surface  is  given  by 


(5.27) 
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<t)\< 


-^Ko)|+- 

V2'  ^  ^'     k 


— t     p 


(5.28) 


Proof:  We  choose  a  positive  definite  Lyapunov  function 


^IT, 


V[t)  =  ^[s'{t)s{t) 


(5.29) 


Our  aim  is  to  show  that  V{t)  <  0  for  the  choice  of  parameters  k  and  e.  The  development 
here  is  fairly  standard. 


utilizing  (5.26)  for  s{t) 


y{t)-s\t)i{t) 


=  s^{t)\k\sit)f&gn[s{t)]  -  d{t)sgn[s(t)]  +  Cf{x,t) 


(5.30) 


From  the  available  bounds  on  the  uncertainty 

\Cf{x,t)\  <  \C\\f{x,t)\  <  m\At)\  ^  d  (5.31) 

Using  (5.31)  in  (5.30) 

V{t)  <  -\sitf[ks'^{t)sgn[s{t)-\]  -  J(0[/(0sgn[5(0]]  +  /(OC/(x,0   (5.32) 
We  observe  that 

s^(r)sgn[5(0]  =  l5/(0sgn[5,(0] 
=  l|5,.(0|[sgn[s,(r)]' 

=  i:K(0|  =  k(0|i  (5.33) 

Using  (5.33)  in  (5.32) 

V{t)  <  -k'^{tf\s{t\  -  e\s{t)\  (5.34) 

Since  both  the  terms  on  the  right  hand  side  are  negative,  the  system  is  guaranteed  to  attain 
the  sliding  mode  in  finite  time.  Furthermore,  for  the  choice  of  a  =  1.0,  (5.34)  becomes 

V{t)<-k\s{tf-e\s{t%  (5.35) 

since  for  a  vector  |  •  L  <  |  •  L . 
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Using  (5.5.29) 

V{t)<-kV{t)-£V^'^{t) 
The  above  is  a  scalar  difference  equation,  and  from  Appendix  F, 

k 


(5.36) 


#W^ 


F(0)+- 


K 


m  < 


\m\ 


42  e 


-2'_V2£ 


Equating  the  right  hand  side  to  zero  we  get  the  result. 
Corollary: 


(5.37) 

(5.38) 
V 


For  the  choice  0  <  a  <  1.0,  a  tighter  result  can  be  obtained.  Following  the  same  steps  as  in 
the  proof  of  the  above  theorem  (see  Appendix  F),  it  can  be  shown  that  for  a  control  law 
(5.24)  with  d{t)  =  P\C\x{t)\,  the  sUding  mode  is  achieved  within  t  =  t,.,  where 

-|0.5(l-a) 


t< 


'     k{l-a) 


1     r 


s'  (0)5(0) 


and 


\s{t)\  <  V2 


{I -a) 


:Wo)] 


0..'5(l-a) 


-kt 


2 


0<t<L 


This  is  advantageous  because,  we  avoid  an  extra  discontinuous  control  element,  e . 
Asymptotic  Stability  in  the  Sliding  Mode 

Note  that  when  the  system  is  in  the  sliding  mode,  the  control  law  reduces  to  a  linear 
state  feedback.  In  this  chapter,  we  follow  a  Lyapunov-Razhamikhin  type  approach  to 
derive  sufficient  stability  condition. 

The  dynamic  equation  in  sliding  mode  is  given  by 

^(t)^{A^,  -A,,K)w{t)  +  TJ{x,t)  (5.39) 

Since  in  the  sliding  mode,  Wi{t)  =  -Kwit),  it  suffices  to  show  the  asymptotic  stability  of 
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(5.39).  The  first  task  is  to  convert  the  uncertainty  description  in  tenns  of  wit) . 


'my 

w,it)_ 


wit) 
-Kw,it) 


=  T 


n—m 

-K 


z{t)  =  T 
Next  we  seek  an  expression  relating  xit)  and  zit).   From  (5.5.3) 


wit) 


(5.40) 


\xit)\<\zit)\+\[_,^e^^'-''-'^B,iuiT)dT 


(5.41) 


where  uit)  =  -iCB)'' CAzit). 
Thus 


\xit)\<\zit)\  + 


max 

-h<e<Q 


AG 


B^^iCB)-'CA\\JziT)\dT 


(5.42) 


The  main  hurdle  at  this  point  is  to  express  |z(t)|  in  terms  of  \z{t)\.    To  obviate  this 
difficulty,  we  take  Lyapunov— Razhamikhin  type  approach.  First  of  all  we  show  that  for 
the  case  of  matched  uncertainties,  proving  asymptotic  stability  is  straight  forward. 
Definition:  Matching  condition  [Hung  et  al.,  1993].  If  the  uncertainties  satisfy 

fix,t)  =  Bfix,t) 
where  fix,  t)  e  9t'"^  ,  then  the  matching  condition  is  satisfied. 
This  means  that  the  perturbation  is  in  the  image  of  B  ■ 

Tfix,t)  =  TBfix,t) 

0 

T^Bjix,t) 

Hence,  in  the  sliding  mode,  the  system  equation  becomes 

m)  =  [Au-^X2Qr'Q)w{t)  (5.43) 

Note  that  (5.43)  is  of  the  same  form  as  the  equation  describing  sliding  mode  for  the 
nominal  system  (5.15).  This  implies  that  the  uncertainties  will  have  no  effect  on  the 
system  in  the  sliding  mode.  This  is  one  of  the  major  advantages  of  the  sliding  mode 
design.  The  state  feedback  is  simply  used  to  place  the  poles  of  the  nominal  system  in  the 
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sliding  mode.  However,  the  assumption  of  matched  uncertainties  is  a  highly  restrictive 
one,  as  for  example  in  the  m  =  1  case,  we  are  constraining  the  perturbation  to  be  a  scaled 
version  of  B .  One  should  note  that  the  situation  represented  in  (5.43)  is  an  ideal  case. 
That  is,  if  the  controller  is  able  to  exactly  cancel  the  perturbation,  then  the  sliding  mode  is 
independent  of  the  perturbation.  This  obviously  suggests  that  if  one  could  estimate  and 
cancel  the  perturbation  (or)  the  effect  of  the  perturbation  on  the  sliding  mode,  then  the 
system  in  sliding  mode  can  be  made  to  be  independent  of  the  matching  perturbations. 
These  kinds  of  approaches  are  commonly  called  perturbation  compensation  and  have  been 
applied  to  linear  state  delay  systems  [Chan  et  al.,  1996].  Even  if  mismatched  uncertainties 
are  present,  perturbation  compensation  can  be  employed  to  reduce  the  magnitude  of  the 
required  discontinuous  control. 

In  the  mismatched  case,  it  is  not  possible  to  completely  cancel  the  perturbation  and 
there  will  be  a  residual  effect  of  the  perturbations  in  the  sliding  mode  and  the  state  feedback 
should  be  such  that  the  overall  system  is  asymptotically  stable.  Theorem  5.2  gives  a 
sufficient  condition  for  the  mismatched  case. 

Before  we  state  and  prove  Theorem  5.2,  we  explain  about  a  difficulty  that  commonly 
arises  in  proving  asymptotic  stability  in  state-delayed  systems  and  input-delayed  systems 
(as  the  control  usually  involves  state  feedback,  the  input-delay  problem  has  relevance  and 
resemblance  to  state-delayed  problem).  There  are  two  classical  approaches  in  the  literature 
for  proving  asymptotic  stability  for  delayed  systems.  Over  the  yeai's  there  have  been  many 
incremental  improvements,  but  the  message  is  that  all  methods  are  conservative.  This  is 
because  the  differential  equations  that  determine  the  dynamics  of  the  system  have  two 
difficult  elements.  One,  the  perturbations.  Two,  the  delayed  states.  We  know  only  the 
norm  bound  on  the  perturbation,  and  while  using  this  sufficiency  creeps  in.  The  problem 
with  delayed  states  is  more  interesting.  In  using  Lyapunov-Razhamikhin  type  approaches, 
a  restriction  is  imposed  on  the  dynamics  of  the  system. 
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Consider  a  Lyapunov  candidate 


where  P  =  P^  >  0 . 


V{t)  =  w'^it)Pwit) 


At  this  point,  we  introduce  the  usual  Razhamikhin  type  restriction  on  the  system.  We 
restrict  that 

V{9)<SV{t)    t-h<9<t,V5>l 

It  is  important  to  understand,  how  restrictive  the  above  assumption  is.  The  reason  for 
using  such  an  assumption  is  that  it  allows  one  to  bound  the  norm  of  the  delayed  states  of 
the  system  in  terms  of  the  norm  of  the  present  states  of  the  system.  Obviously  any 
condition  based  on  this  will  only  be  sufficient.  Since  the  following  condition  holds  for  any 
positive-definite  matrix  P  =  P^   >   0,  and  w  e  9?""'" 

it  follows  that 

\w{e)\  <  5a\w{t)\ 
where  a 


_    l^m.AP) 


^rainiP)  ' 

The  above  condition  is  not  very  restrictive  as  it  is  true  for  some  5^  >  5  >  1,  for  all 
stable  systems  whose  state  norm  continuously  decreases  (and  it  is  easy  to  show  that  for 
unstable  systems,  the  above  is  true  for  any  positive  d  ).  The  actual  values  of  5  and  a 
(and  hence  P)  quantify  the  decay  rate  of  the  system.  This  condition,  however  could  be 
violated  for  a  short  durations,  for  example,  when  a  delayed  disturbance  hits  the  system 
(i.e,  the  norm  of  the  states  could  increase).  However,  the  controller  will  act  against  the 
disturbance  and  restore  the  performance.  Hence,  temporary  violation  of  this  constraint  is 
not  fatal,  and  will  not  affect  the  stability  conclusion. 
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Theorem  5.2:  The  system  (5.4)  in  sliding  mode  with  the  mismatched  uncertainties 
(5.2)  is  asymptotically  stable  if  the  following  condition  is  satisfied 


where 


iS 


n-m 

-K 


\  +  ha 


max 

-/z<0<O 


AB 


BACBY^CA 


<  ^"''"f^,)      (5.44) 


a  = 


^max(^) 


_  '"max 


^min(^) 


and  P >0  is  the  solution  of 


P[A,,-A,^K)  +  (Au-AnKfp  =  -Q,  Q>0 


(5.45) 


Proof.  Under  the  similarity  transformation  (5.8)  the  system  (5.4)  transforms  to 
(5.11).  If  the  uncertainties  are  present,  using  (5.11)  equation  (5.11a)-(5.11b)  adopts  the 
form 


w{t)^A^iw{t)  +  Ai2Wi{t)  +  TJ{x,t) 


(5.46) 


Wiit)  =  A2iw(0  +  A22Wi{t)  +  Biu{t)  +  T2f{x,t)  (5.47) 


Consider  a  Lyapunov  candidate 

V{t)-^w^{t)Pw(t)  (5.48) 

where  P  =  P^  >  0,  P  =  P^  >0,P  e  gt^"-'")^^"-'") . 
Taking  the  derivative  of  (5.48) 

no  =  w^it)[p{Ai  1  -  Ai2K)  +  (Ai I  -  A,2/i:)^Pp(r)  +  2w^(t)PTJ{x,t)(5.49) 

At  this  point,  we  make  the  usual  Razhamikhin  type  restriction  on  the  system.  We  restrict 
that 


V{d)<5V{t)    t-h<0<t,yd>l 


(5.50) 


it  follows  that. 
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\wid)\<5  Hr^4^\w(t)\ 


^mm(^) 


(5.51) 


Denoting  a  - 


'   '"^'^^'^^  using  this  in  (5.42) 


^min(^) 


k(0|  < 


n-m 

-K 


\  +  h5a 


max 

-/!<0<O 


AQ 


BACBY'CA 


k(0|      (5.52) 


Therefore, 


\f{x,t)\<l5 


-K 


1  +  h5a 


max 

-/i<£><0 


AO 


BACBrCA 


\w{t)\  (5.53) 


We  combine  all  the  terms  to  get 


A=iS 


'  n-rn 

-K 


\  +  h5a 


max 

-/!<0<O 


AB 


D^-^ 


BACBY'CA 


\f{x,t)\<p,\w{,t)\ 


(5.54) 


hence 


V{t)  =  w^{t)\p{A^  1  -  A12K)  +  (A|  1  -  A,  2Kfp]wit)  +  2/3i  ||pH7|||>v  (Op(5.55) 
Choose  P>0,   Q>0,  such  that 


P{A,,-A,2K)  +  {A,,-A,2KyP  =  -Q 


tttT 


<-w'  {,t)Qw{t)  +  2j5MTM{t)\ 


(5.56) 


no  ^  [-^min(e)  +  2A  Il^ll7i|lll^(0p  (5.57) 

Hence  a  sufficient  condition  for  asymptotic  stability  in  the  sliding  mode  is  given  by 

^nun(e) 


A< 


2IIPI  7; 


I.e., 


-K 


l  +  h5a 


max 
-h<e<o 


A0 


BACBY'CA 


^^mML     (5.58) 

211^1  rJ 
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Note  that  if  the  above  expression  is  true  for  5  =  1  then  there  exists  a  5>l,  sufficiently 
small  that  it  will  still  hold  true.  Therefore,  a  sufficient  condition  for  stability  is  given  by  the 
above  expression,  with  5  =  1.  V 

Remarks: 

(i)      In  the  sliding  mode,  the  derivative  of  the  Lyapunov  function  is  given 

nt)<[-Kin{Q)+m\p\tTi\\iw{tf  <o. 

Denoting  e  :=  -{-^miniQ)  +  2/?i ||P||||ri ||) ,  we  can  obtain 


:i— /.^|2 


V{t)  <  -£\witf 
£ 


<- 


V{t) 

(5.59) 


.KiniP) 

This  explains  why  we  can  get  only  asymptotic  stability  guarantee,  that  is,  we 
can  only  state  that  the  norm  of  the  states  will  reach  zero  as  time  goes  to  infinity. 
In  the  next  section  we  will  elaborate  more  about  this  issue, 
(ii)  For  a  given  system  description,  once  the  poles  of  the  closed-loop  nominal 
system  are  chosen,  everything  except  P  and  hence  Q  is  known.  Hence  P  is 
an  optimization  parameter.  We  remark  that  the  best  choice  is  to  make  sure  that 
the  condition  number  of  P  is  minimized  while  satisfying  (5.44)-(5.45).  For 
large  systems,  this  constrained  optimization  problem  can  be  cast  as  an  convex 
LMI  problem  and  easily  solved  through  available  commercial  software 
[Niculescu  er  (3/.,  1997]. 

Theorem  5.3.   If  the  uncertain  system  (5.4)  is  asymptotically  stable,  it  implies  the 
asymptotic  stability  of  system  (5.1). 

Proof.     From  the  transformation  (5.3),  we  have 

.Ae 


\xit)\<\z{t)\  +  h 


max 
-h<e<o 


e 


■^hKI  (5-60) 


From  Theorem  5.2,  the  transformed  system  (5.4)  is  asymptotically  stable.    Hence, 
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lim  \zit)\  -^  0.  Since  u{t)  =  -(CB)  '  CAzit),  in  the  sliding  mode,  lim  |«(0|  -^  0     From 


t^cx, 


t^oo 


(5.60)  the  result  follows.  V 

Implementation  Issues 

Due  to  inherent  limitations  such  as  the  switching  delay,  system  inertia,  etc.,  the  sliding 
motion  does  not  take  place  in  any  real  system.  Instead,  the  trajectories  chatter  in  a 
neighborhood  of  the  switching  surface.  The  chattering  phenomenon  is  undesirable  because 
it  may  adversely  affect  the  actuating  mechanism  and  excite  the  high-frequency  unmodeled 
plant  dynamics,  leading  to  instability.  The  discussion  about  chattering  and  the  ways  to 
alleviate  this  problem,  as  outlined  in  Chapter  4,  applies  to  the  present  case  of  stabilizing 
systems  with  input-delay  as  well. 

5.5    Some  Open  Issues  in  Sliding  Mode  Literature 

Although  the  sliding  mode  control  approaches  have  been  used  for  a  long  time  for  a 
variety  of  systems,  there  are  nevertheless  some  open  issues.  We  raise  the  following  issues 
of  interest: 

(1)  Is  it  possible  to  guarantee  that  the  norm  of  the  states  continuously  decrease 
during  the  reaching  phase?.  If  so,  what  is  the  non-trivial  upper  bound,  on  the 
norm  of  the  states? 

(2)  Is  it  possible  to  get  an  upper  bound  the  norm  of  the  states  when  the  system  is 
in  the  reaching  phase  (other  than  of  course,  the  trivial  estimate  \x{0)\ )  ? 

(3)  We  are  able  to  prove  that  the  sliding  mode  will  be  attained  within  a  finite  time. 
At  the  same  time,  the  best  we  can  prove  for  the  system  in  the  sliding  mode  is 
asymptotic  stability.  We  also  seek  to  explain  why  this  is  so. 

We  also  remark  that  these  questions  apply  to  all  sliding  mode  control  schemes  (for 
linear  uncertain,  Uncertain  linear  systems  with  state  delay  and  so  on..)  and  not  just  for  the 
input-delay  case  that  we  consider  here.  Note  that  the  usual  approach  in  the  literature  for 
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sliding  mode  systems  is  to  prove  the  reaching  condition  within  a  finite  time  (as  we  did  via 
Theorem  5.1),  and  then  prove  asymptotic  stability  for  the  reduced  system  in  the  sliding 
mode  (as  in  Theorem  5.2).  However,  there  are  problems  with  this  approach.  For 
example,  one  does  not  have  an  estimate  of  the  norm  of  the  states,  when  the  system  reaches 
the  sliding  mode.  Hence,  we  do  not  know  the  proper  initial  condition  for  the  scalar 
Lyapunov  differential  equation  given  by  (5.59).  One  only  hopes  that  it  is  less  than  the 
norm  corresponding  to  the  initial  condition.  Even  assuming  so,  it  is  still  highly 
conservative.  What  we  need  is  a  way  of  estimating  a  tight  upper-bound  on  the  norm  of  the 
states,  during  the  reaching  phase.  If  we  obtain  such  an  expression,  then  we  will  be  able  to 
answer  questions  (Q1)-(Q2). 

Before  we  seek  to  get  an  analytical  expression  to  the  above  questions,  it  is  instructive 
to  ponder  over  the  significance  of  Theorem  5.1.  What  (5.27)-(5.28)  imply  is  that  the  norm 
of  a  linear  combination  of  the  states  continuously  decreases  and  attains  zero  within  a  finite 
time.  But,  a  linear  combination  of  two  time-dependent  functions  (here  states  of  the  system) 
may  be  decreasing  towards  zero,  while  the  individual  elements  may  be  diverging  in 
opposite  direction,  making  their  sum  (or  difference)  approach  zero.  Hence,  it  seems  like 
that  we  do  not  even  have  a  guarantee  that  the  norm  of  the  states  of  the  system  will 
continuously  decrease,  when  the  norm  of  the  switching  function  continuously  decreases 
and  attains  the  sliding  mode.  In  other  words,  the  central  question  here  is,  can  one  obtain 
sliding  mode  when  the  system  is  unstable  . 

However,  upon  reflection,  one  simple  fact  emerges.  Consider  the  linear  uncertain 
system  (5.4)  and  bounded  perturbation  (5.2).  The  solution  will  be  driven  by  the  interaction 
of  the  control  and  perturbations  and  the  rate  at  which  the  states  converge  (or)  diverge  will 
depend  on  the  state  transition  matrix,  and  whether  or  not  it  is  stable  enough  to  overcome 
the  perturbations.  Ultimately  the  form  of  the  solution  will  contain  complex  exponentials, 
apart  from  other  terms.  Since,  the  perturbation  is  bound  by  the  norm  of  the  states,  an 
unique  solution  does  exist.  It  is  not  our  intention  to  derive  it.  In  fact,  it  is  impossible  to 
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derive  it,  unless  the  form  of  the  perturbation  is  exactly  known.  Now,  Theorem  5.1  implies 
that  the  norm  of  a  linear  combination  of  the  complex  exponential  terms  is  continuously 
decreasing.  If  the  system  is  unstable,  it  is  the  real  part  of  the  complex  exponentials  that 
will  dominate  the  system  dynamics.  In  that  case,  the  resulting  linear  combination  of  the 
complex  exponentials  will  diverge.   Hence,  it  is  not  possible  for  \s{t)\  to  continuously 

decrease  and  attain  zero  within  a  finite  time,  if  the  system  is  unstable. 

This  however,  does  not  necessarily  mean  that  the  norm  of  the  states  will  continuously 
decrease.  For  example,  for  a  system  with  state  delay,  when  the  delayed  states  first  come 
into  action,  there  could  be  a  temporary  increase  in  the  norm.  Alternately,  there  could  be  a 
non-linear  uncertainty  that  satisfies  all  our  constraints  and  still  has  a  spike  leading  to  an 
increase  in  the  norm  of  the  states.  Of  course,  the  controller  will  seek  to  react  to  this 
situation  and  bring  the  system  under  control.  Therefore,  although  mathematically  not 
precise,  it  is  nevertheless  acceptable  to  assume  in  general,  that  the  norm  of  the  states 
decrease  in  the  reaching  phase.  This  conclusion  is  important  to  assess  the  behavior  of  the 
state  norm  in  the  reaching  phase.  For  example,  without  such  a  guarantee,  it  is  impossible 
to  bound  the  discontinuous  control  signal  d{t).  From  (5.24),  d{t)  =  P\C\\x{t)\  +  £.  Hence, 

the  value  of  the  discontinuous  control  is  linked  to  the  norm  of  the  states.  In  proving 
asymptotic  stability  in  the  sUding  mode,  we  obviate  this  difficulty,  as  &gn[s{t))  =  0.   With 

the  above  conclusion,  however,  one  could  make  a  confident  assumption  on  the  expected 
value  of  d{t),  much  in  the  same  vein  of  the  Razhamikhin  type  assumption,  that  seeks  to 
restrict  the  norm  of  the  states. 

The  following  lemma  is  needed  for  the  ensuing  developments. 

Lemma   5.1 
Consider  the  following  scalar  differential  equation 

x{t)<-ax{t)  +  {btx^{-ctl2)  +  e  exp[-S2t/2)- d)x^'^{t)  c^  f,]\ 

with  the  initial  condition  Jt(0)  and  the  following  constraints: 
(1)    The  constants  {a,b,c,d,e)  >  0  and  a^c  and  a^  82 
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(2)    x>0,  Vt  >Qand4x>Q,  \ft  >0 
Then, 


4^)^ 

V4o)  +  - 

_                          6 
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Proof: 
Remark: 

See  Appendix  G 

fcAJJ 

~'n"2Vj 

a 

(5.62) 


In  (5.62)  the  first  three  terms  are  positive,  while  the  third  term  is  a  negative  constant. 

Furthermore,  the  first  term  continuously  decreases  from  the  initial  value,  while  the  second 

and  terms  always  have  a  maxima.  As  a  parallel  we  remark  that  the  second  and  third  terms 

in  (5.62)  look  like  the  response  of  an  under  damped  second  order  system  to  unit  impulse 

input. 

Upper  Bound  on  the  Norm  of  the  States  During  the  Reaching  Phase 

In  this  section,  the  aim  is  to  find  an  upper  bound  on  the  norm  of  the  states,  for  the  stable 
closed-loop  system.  The  analysis  for  the  nominal  system  is  straightforward  and  is 
presented  in  the  sequel.  However  for  the  perturbed  system,  the  analysis  is  quite  involved. 
In  the  reaching  phase,  the  dynamic  equations  ai'e  given  by 


nt)  =  (A,  I  -  A^2K)yv{t)  +  TJ{x,t)  +  Ai2s(r) 


w  ,{t)  =  s(,t)  -  Kw  it) 


(5.63) 


(5.64) 


In  complete  parallel  to  the  proof  of  Theorem  5.2,  consider  a  Lyapunov  candidate 


V{t)  =  w'^{t)Pw{t) 


/\  ^  '■ 


where  P  =  P^  >  0,  P=^  P^  >  0,  P  e  9t(«-'«)x("-"0  and  impose  the  restriction 


V{0)<SV{t)    t-h<6<t,\f5>i 


(5.65) 


(5.66a) 


It  is  straightforward  to  show  from  (5.66a)  that 


\w{e)\  <  5 


\     -mm  \'  I  J 


W 


(01 


Choose  P>0,   Q>0,  such  that 
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(5.66b) 


(5.67) 


w^it) 


(5.68) 


Taking  the  derivative  of  (5.65)  yields 

V{t)  =  vv ^(0[p(Ai  1  - A^oK)  +  (A,,  - A^^Kfp 
+  2w'^{t)PTJ{x,t)  +  2h'^(OPAi2s(0 
Converting  the  Uncertainty  in  Terms  of  w{t) 

Again,  the  challenge  is  to  convert  the  uncertainty  description  in  terms  of  w{t).     Since 
|/(A:(r),?)|  <  P  |x(r)|,  we  run  into  a  problem,  since  the  stated  aim  is  to  find  an  upper  bound 

on  \x{t)\.    In  the  ensuing  development,  we  will  pinpoint  this  difficulty  and  propose  an 

additional  assumption  on  the  perturbation  that  will  enable  us  to  obviate  this  hurdle..  From 

Theorem  5. 1,  we  have 


For  simplicity,  we  denote 


Hence  \s{t)\<be  ^  -d. 
In  the  reaching  mode, 


k(ON 


|s(0)|  + 


42  £ 


2    _ 


V2e 


b  = 


Ko)|+ 


42  £ 


c  =  k,    d  = 


2e 


z(0  =  r-V(0  =  T 


■w{ty 
w,{t) 


z{t)  =  T 


'n-m 

-K 


w{t)  +  T 


-1 


w{t) 
s{t)  -  Kw{t) 

S{t) 


(5.69) 
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MOI 


< 


n-m 

-K 


\w{t)\+  %    1^(01 


where  T2  is  a  matrix  of  dimension  nxm,  partitioned  out  of  T    . 
Since  we  used  the  transformation  (5.3),  from  (5.42) 

\xit)\<\zit)\  +  \lli^e^^'-'-'^BMt)dT 
where  the  controller  is  given  by 

uit)  =  -(CB )"'  [CAz{t)  +  k\sit)\sgn[sit)]  +  d{t)sgn[sit)]\ 
since  |sgn[s(0]|  ^  -Jm ,  we  take  an  upper  bound 


1^^(01  ^ 


\x(t)\<\z{t)\  + 


.,^\-l 


(CB)    CAz{t)+-^  [CB) 


kbe  -  -kd  +  d{t) 


max 

-h<e<o 


A0 


B,i(CBrCA    1;_Jz(t)|Jt 


+  vm 


max 

-h<e<Q 


AO 


B., 


iCB)-^    Jl 


kbe  2    -kd  +  dij) 


dt 


Substituting  for  \z{t)\  in  terms  of  |M'(r)|  and  \s{t)\,  and  using  Razhamikhin  type  assumption 
for  the  delayed  states,  via  (5.66b),  results  in 


|x(0| 


< 


n—in 

-K 


\w{t 


T, 


+  h5a 


7, 


+ 


+  'Jm 


max 

-/z<0<O 


max 

.-/!<0<O 


n-m 

K 

Ad 


AO 


max 

-h<e<Q 


y\  T  A 

he   2    -d 

V  J 

AQ 


BjiCBy^CA   \w{t)\ 
B,iiCB)-'CA    j^ 


h-h 


f      c        ^ 

be  ^    -d 


dx 


(5.70a) 


\\B, 


(CB)-'     t 


V  J 

c 


kbe  2    -kd  +  d{r) 


dx 


In  the  above  equation,  the  only  unknown  quantity  is  d{t).   According  to  our  control  law, 
d{t)  -  j3  ||C|||x(f)|  +  e,  and  our  stated  aim  is  to  find  an  upper  bound  on  the  norm  of  the 
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states  (note  that  this  difficulty  does  not  arise  for  the  nominal  system).  This  difficulty  can 

be  overcome,  if  an  additional  assumption  can  be  made  on  the  uncertainty  description.  We 
assume,  that  the  uncertainty  vector  |/(jc(r),f)|  has  finite  energy.  Hence,  it  can  be  bounded 

by  a  decaying  exponential  and  can  be  represented  as  d{t)  =  a-^  Qxp^-S^t),  ^2  >  0.    It  is 

assumed  that  the  actual  value  of  82  is  known. 

To  simplify  notation,  we  denote 

„A0 


ai  := 


max 

-/!<0<O 


BJCBT'CA 


a-i  '.=  k-^Jm 


max 
-/i<6i<0 


AO 


B, 


(CB) 


0-2  :=  (aj  +  a2  )j5dh 
\f{x,t)\<P\x{t)\ 


|/(x,r)|<A|vv(0|  +  /3(ai+a2) 


+  /3(|r2  l^exp-^ 


^2b       (ch^ 
— exp 
c 


-1 


V  ^  y 


exp 


-ct 


■^3 


../+/? 


.       V2 

fexp(52^)  ~  ^]^^v{~^i^) 


where  /3i  =  /3 


"■n-in 

-K 


l  +  hSa 


max 

-/i<e<o 


Ad 


BACBT'CA 


To  further  simphfy  notation,  we  denote 


m  =  /3(ai  +^2) 


2b 

—exp 
c 


^1-1 
V2 


A 


+  I5b 


p  =  a2+  Pd 


—                 ( -ct\ 
\f{x,t)\<P^\w{t)\  +  mtx^  — p  +  n&x^{-52t) 

V  2  y 


(5.70b) 


From  (5.67),  (5.68)  and  (5.70b)  we  have 
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V{t)<-w^ {t)Qw{t)  +  2\w{t)\   P   (||r||  |/(;c,0|  +  ||Ai2||  \s{t)\)       (5.71a) 


^(0    <  -A^in(e)  +2AK0I   P  \\m\witf+2\wit)\P   IIA12II    boxv 


—  -a 

2  J       ) 


+  2^(01 


^   .       f-ct'^ 
mexp 

V  ^  y 


(5.71b) 


-  p  +  nexp(-520 


Denoting 


«  =  --^min(2)+2AK0| 


ll^lll/^maxl^ 


Z?  =  2 


|il|  m+  Aj9  o 


'^max  I  " 


K^.[P) 


d  =  2\T^\p 


^max  I  " 


VV^^inC,. 


^max(^) 
^min(^) 


equation  (5.71b)  can  be  written  as 


V{t)    <-aAn,axU'    k(0|   +  ^exp 


-ct 


-d 


y^mx^,[p) 


+ 


^exp^— ^ 


i^(OUA^in(^) 


^  ^_^A  /_^2? 


V(r)    <-a:y+   texp  +e  exp 

V  V    2    y* 


d 


(5.72) 


We  remark  that  (5.72)  is  in  a  form  similar  to  (5.61)  and  hence  it  is  straightforward  to  get  a 
bound  on  the  norm  of  the  states  during  the  reaching  phase. 
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Implication  of  the  Result 

Using  Lemma  5.1,  for  the  system  (5.72),  the  norm  of  the  states  can  be  bounded.  The 
maximum  value  of  the  upper  bound  and  the  time  at  which  it  occurs  can  also  be  easily 
computed,  since  the  constants  are  known.  This  explicitly  answers  questions  (1)  and  (2). 
Upper  Bound  on  the  Norm  of  the  States  for  the  Nominal  System 

For  the  nominal  system,  from  (5.71a),  it  follows  that  the  following  equation  holds. 

V{t)<-w'^it)Qw(t)  +  2\wit)\  \p\  |Ai2|  \s{t)\ 
Hence, 


V{t)<-^min{Q)Htf+Mt)\ 


^exp|  -Y 


-rt\  '^  1 


■d 


m  < 


KJpmt  + 


il2 


Vin(^) 


bexp 


f-rt\       ^\ 


-Ct 


■d 


V        V  ^  y       J 


^r.in(^)KO| 


y(0<- 


'^minM^ 


^max(^) 


no+ 


M2 


VV^--(^/y 


-d 


^exp 

V         \  ^  J       J 


no 


Then,  from  Lemma  5.1,  it  follows  that  an  upper  bound  on  the  norm  of  the  states  can  be 
established  for  the  nominal  system,  without  resorting  to  any  additional  approximations. 
Caveat 

It  is  easy  to  show  that  the  solution  of  (5.72)  goes  to  zero  within  a  finite  time.  Does 
that  mean  that  we  are  proving  that  the  norm  of  the  states  are  guaranteed  to  go  to  zero  within 
a  finite  time.  The  temptation  to  reach  such  a  conclusion  must  be  avoided.  This  is  because, 
(5.72)  is  valid  only  until  the  time  sliding  mode  is  attained.  After  that  it  is  the  differential 
equation  (5.39)  and  the  stability  condition  (5.58)  that  hold  true. 
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Why  only  Asymptotic  Stability  Guarantees? 

Now  we  turn  our  attention  to  the  third  issue  that  we  raised.  Why  is  it  not  that  we  are 
able  to  prove  that  we  can  drive  a  m  dimensional  manifold  of  the  system  to  origin  within  a 
finite  time,  whereas  it  is  impossible  to  get  such  a  resuh  for  the  remaining  (n-m) 
dimensional  manifold  of  the  system.  The  answer  lies  in  the  fact  that  we  are  limited  by  the 
rank  of  the  matrix  B.  The  following  theorem  is  a  step  towards  understanding  the  inherent 
system  limitation. 

Theorem  5.4.    Consider  the  following  nominal  system 

x{t)  =  Ax{t)  +  Bu{t) , 

(1)  If  the  control  matrix  B  is  of  full  rank,  the  system  trajectories  can  be  driven 
to  the  origin  within  a  finite  time  by  the  discontinuous  control  law 

u{t)  =  -{By\Ax{t)  +  a^x{t)  +  a2?>%n{x{t))\  a,  > 0,  «2>0 

(2)  If  the  control  matrix  B  is  not  of  full  rank,  then  it  is  not  possible  to  drive  the 
system  trajectories  to  the  origin  within  a  finite  time  by  a  discontinuous  control  law. 

Proof:  First  consider  the  case,  where  B  is  of  full  rank.  Then,  applying  the  above 
discontinuous  control  law  results  in 

x(t)  =  -tt]  x(0  -  a2  sgn{x(t)) 
Introduce  a  simple  Lyapunov  function  of  the  form 

V{t)^^x'^{t)x{t) 

V{t)  =  x^  it)x{t) 

=  -2ai  V{t)  -  a2 x^(r)sgn(x(0) 
Since  x'^ {t)sgn{xit))  -  \x{t)l ,  the  above  equation  can  be  rewritten  as 


V{t)  <  -2^1  V{t)  -  V2a2  V^ 
Following  Appendix  F,  it  is  straightforward  to  show  that  the  trajectories  of  the  system  will 
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origin  within  a  finite  time.  Furthermore,  it  is  also  possible  to  get  an  explicit  upper  bound 
for  the  norm  of  the  states.  This  completes  the  first  part  of  the  proof. 

Now  consider  the  case  where  B  is  not  full  rank.  Consider  the  discontinuous  control 
law 

u{t)  =  -K^xit)  -  K2  sgn(x(0) 
The  system  dynamics  is  given  by 

x{t)  =  Ax{t)  -  BK2  sgn(x(0)  ^  A={A-  BK) 
Consider  the  usual  Lyapunov  function 


V{t)^x^{t)Px{t) 


V{t)  =  x\t) 


A^P  +  PA 


x{t)  -  2x^{t)PBK2  sgn{x(t)) 


Denoting  Q  =  ■ 


A'  P  +  PA 


V{t)  =  -x^{t)Qx{t)  -  2x'^{t)PBK2  sgn(x(0) 
=  -^n^n{Qy{t)x{t)  -  2x\t)PBK2  sgn(x(r)) 


< 


V{t)  -  2x^{t)PBK2  sgn(x(r)) 


In  order  to  prove  that  the  system  will  attain  the  origin  within  a  finite  time,  the  second  term, 
i-2x^ {t)PBK2Sgn{x{t))\  should  be  negative-definite.  We  will  show  by  a  simple  generic 

second  order  counter-example  that  this  term  is  indeed,  sign-indefinite. 

-2x'^{t)PBK2  sgn(x(0)  =  -2x'^(0^sgn(x(0) 
We  note  that  any  singular  matrix  of  dimension  two,  can  be  written  as 


'a     b' 

a 

b' 

p= 

= 

_c     d_ 

cm 

ab 

(This   assertion  can  be  proved  easily.      Since    P    is    assumed   to   be    singular, 
detP  =  ad-cb  =  0.  Since  3a  €  9t,  such  that  c  =  aa,  it  naturally  follows  that  d  =  ab) 
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-x^{t)PBK2  sgn(x(0)  =  -x^{t) 


a       b 
aa    ab 


sgn(x(0) 


Since 


sgn(A:(r))  =  [sgn{x^(t)) ; sgn(^2(0)]  , 

Xi  sgn(x,-(0)  =  |-^/(0|  and 

xi  sgn(xy(r))  =  \xiit)\  sgn{xi{t))  sgn[x jit) 

-x^{t)PBK2 sgn{xit))  =  -{a  +  b)  [|xi(r)|  +  a\xo{t)^  [l  +  sgn(xi (r))sgn(x2(0)] 

In  the  above  equation,  [l  +  sgn(xi(r))sgn(x2(0)]  ^  0  (in  fact,  it  can  take  the  values 

0,1,2). 

-X   {t)PBK2  sgn(jc(r))  >  0  under  the  following  two  conditions, 

{a  +  b)>Q,  a  <  0,  \x^it)\  <  \x2{t)\  for  some  t  >  0 

(a  +  6)<0,  a>0 
and  -x   {t)PBK2  sgn(x(0)  ^  0  under  the  following  two  conditions: 

{a  +  b)>0,  a>0 

{a  +  b)<0,  a<0,  |x,(0|<|x2(0|  for  some  z^>0 
Hence,  -2x   {t)Ps,gn{x(t))  can  not  be  guaranteed  to  be  negative,  an  essential  pre-requisite 
to  prove  finite  time  stability.  The  reason  for  this  inability  is  of  course,  the  fact  that  B  is  not 
full  rank.  V 


In  fact,  we  earlier  exploited  the  spirit  behind  Theorem  5.4  to  prove  the  reaching 
condition  for  the  system  (5.26).  To  see  this  more  clearly,  let  us  revisit  the  equations 
(5.21)-(5.23)  for  the  uncertain  system. 

W{t)  =  A,  ,w  (r)  +  A^2Wj{t)  +  TJ{x,t) 


w,  it)  =  Aj  |W  (r)  +  A22H'/  (r)  +  Bjuit)  +  T2f{x,  t) 
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Since  the  matrix  Bj ,  is  known  to  be  full  rank,  any  linear  combination  of  m  states  of  the 
system  can  be  driven  to  and  maintained  at  zero.  Hence,  we  construct  a  m -dimensional 
manifold  ('switching  function')  as 

s{t)  =  Kwit)  +  Wiit) 
In  order  to  drive  and  maintain  these  m  states  to  zero  {i.e.,  sliding  mode),  we  require  that 

sit)  =  0,  s{t)  =  0 

s{t)  =  [K(Au  -Ai2^)  +  (A2,  -A22K)]w{t)  +  {KA^2  +A22XO 
+  Biu{t)  +  {KTi+T2)f{x,t) 

Since  Bj  is  of  rank  m ,  we  exploit  this  and  design  a  controller 

uit)  =  -{Biy\K(A^ ,  -  A^2K)  +  (A2,  -  A22K)]w{t)  +  [KA^  +  A22XO 
+  {-B,y\{k\s(t)\  +  d{t))  sgn{s{t))] 

This  of  course,  assures  sliding  mode  and  the  reader  is  referred  to  the  proof  of  Theorem 

5.1,  for  further  details.  The  important  message  is  the  following:    When  the  number  of 

control  inputs  is  less  than  the  number  of  states,  i.e.,  in<n,  the  best  we  can  achieve  with 

the  given  physical  system  is  to  drive  a  chosen  m  dimensional  subspace  of  the  system  to  the 

origin  within  a  finite  time  and  maintain  it  there  for  ever.    Of  course,  we  can  assure 

asymptotic  stability  for  the  signals  corresponding  to  the  remaining  {n  -  m)  dimensional 

manifold,  if  the  original  system  is  controllable.  Since  we  have  an  uncertain  system,  there 

is  also  the  additional  challenge  of  having  to  overcome  the  effect  of  uncertainties. 

The  beauty  of  the  sliding  mode  control  is  that  we  achieve  these  twin  objectives  via 

design  of  our  controller  (5.24).  The  choice  to  be  made  then,  is  the  design  of  an  appropriate 

m  dimensional  subspace.    There  are  two  scenarios:  one,  we  are  free  to  chose  any  m- 

dimensional  subspace.    In  such  a  case,  we  choose  the  subspace  to  guarantee  that  the 

nominal  system  in  the  sliding  mode  is  asymptotically  stable.  This  determines  the  choice  of 
K,  and  hence  C  =  \K  /,„]r.  We  then  design  a  controller  to  drive  the  chosen  subspace 
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with  the  desired  speed  to  origin  (i.e.,  to  attain  sliding  mode).  Another  scenario  is  where 
the  switching  function  is  available  in  the  form  of  a  performance  objective.  For  example, 
consider  the  case  where  the  objective  is  to  reduce  the  norm  of  a  given  Unear  combination  of 
the  state  deviations.  In  this  case,  C  is  specified  and  K  is  obtained  via  the  appropriate 
partition  of  CT~  .  Then,  the  conditions  that  we  have  developed  in  this  chapter  aid  us  to 
ascertain  the  feasibility  of  meeting  such  a  requirement.  Of  course,  a  more  useful  constraint 
would  consider  the  combination  of  individual  state  norms,  rather  than  a  norm  of  a  linear 
combination  of  the  states.  This  aspect  merits  further  attention. 


5.6   Constraints  on  the  Control  Input 


When  the  initial  condition  is  known,  we  can  find  an  upper  bound  on  the  norm  of  the 
control  input.  This  is  important  as  in  many  cases  of  practical  interest,  there  are  hard 
constraints  on  the  norm  of  the  control  input  u{t)  (for  example,  in  chemical  engineering 
process  applications,  the  control  valve  as  hard  bounds  of  0%-100%  valve  opening.  There 
may  also  be  safety  considerations). 

First,  we  get  an  upper  bound  on  the  norm  of  the  control  input.  We  have  the  control 
law 


u{t)  =  -{CB)  ' \cAz{t)  +  k\sit)fsgn[sit)]  +  J(r)sgn[s(0] 


(5.73) 


For  the  sake  of  simplicity,  we  write 


where 


u,onti^)^-{CBy\cAz{t)] 
^discontit)  =  -{CBy'\k\s(t)\"sgn[s(m  +  d{t)sgn[sim 


(5.74) 
(5.75) 
(5.76) 


Using  (5.7),  (5.11)  and  (5.19),  (5.75)  can  be  written  as 
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u 


cont 


[t)  =  -[Bi)  \k  /m]  TAT~^ 


In-  m 
-K 


w{t)  + 


On-m 

s{t) 


(5.77) 


ii,antit)  =  -{BiT\K{Au-A^2K)  +  {A2i-A22K)]w{t)   -(fi;)-'[^i2 +^22>(0     (5.78) 
Define: 

Uj^,  :^-{Bi)-'[k{A^^-A^2K)  +  {A2i  -A22K)],  U,2  :=-{B,y\KA,2- A22]  (5.79) 

Then  it  follows  that 


\u{t)\<\\Uiiilw{t)\  +  \\Ui^2\Ht)\+  {BiY^  [k\sit)f+dit) 
From  Theorem  5. 1,  we  have  that 

V2e 


sgn[s(r)]      (5.80) 


<t)\< 


k(o)  +  - 


k 


k  r- 


and  hence,  max|s(r)|  =  |s(0)|.    Also  noting  that  \sgn[s{tyi<  4m ,  and  d{t)<d{Q),  and 
utilizing  the  previously  developed  expression  for  \vv^^^  {t)\ ,  we  have 


|«(0|<^i||Kax(0|+|felHo)|+V^  (^/)    [^k(0)r  +^(0)]     (5.81) 

By  focusing  attention  on  the  above  equations,  it  is  clear  that  there  is  a  price  to  be  paid  for 
placing  the  poles  of  the  nominal  system  in  the  sliding  mode  via  feedbaclc  control  matrix  K 
arbitrarily  far  from  the  imaginary  axis.  The  expression  (5.81)  gives  a  practical  way  of 
balancing  fast  response  to  allowable  values  of  control  signals.  Unfortunately,  due  to  the 
complexity  of  the  controller  and  the  stability  conditions,  it  is  not  possible  to  develop  a 
systematic  way  of  enforcing  a  given  controller  constraint.  However,  for  the  interested 
reader  we  mention  that  such  approaches  have  been  studied  for  simple  state  feedback 
control,  via  Linear  Matrix  Inequalities  based  approaches  [see  Niculescu  et  al,  1997]. 

In  some  applications,  it  is  more  advantageous  to  have  an  expression  for  the  upper 
bound  on  the  norm  of  the  individual  control  inputs.  For  example,  consider  the  well 
studied  case  of  control  of  exothermic  CSTR  where  temperature  and  concentration  are 
controlled  by  manipulating  the  inlet  concentration  and  temperature  of  the  feed ,  and  cooling 
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water  temperature.    Obviously,  there  are  different  constraints  on  tliese  manipulated 
variables.  We  will  get  an  expression  for  the  upper  bound  on  the  individual  control  inputs. 
Very  little  further  work  is  required,  and  we  will  exploit  our  previous  developments. 
For  any  matrix  R,  define:  /?;•:=  ith  row  of  R. 

Then  it  simply  follows  that 


An  additional  advantage  of  the  expressions  (5.81)-(5.82)  is  that,  given  a  priori  upper 
bounds  on  the  control  signals  (for  example,  a  pneumatic  valve  can  only  take  values 
between  0%  -  100%),  they  can  be  used  to  get  an  estimate  (though  conservative)  of  the 
region  of  attraction.  Towards  this  we  note  that  \z{0)\  =  \x{0)\  since  m(t)  =  0,  T  <  0,  and 
hence  \w{0)\  =  |r,z(0)|  =  \T^x{0)\.  Then  from  (5.81)-(5.82),  we  can  get  an  upper  bound  on 
\x{0)\. 

5.7   Extension  to  the  Case  of  Multiple  Delays 


The  control  scheme  presented  in  this  chapter  can  be  extended  in  a  strightforward  fashion  to 
the  case  of  multiple  known  delays.  Consider  the  uncertain  input-delay  system 


x(t)  =  Ax{t)  +  Bu{t)  +  'ZB^,u{t  -hi^)  +  f{x,t) 

k=\ 

To  convert  this  into  a  system  without  delay,  consider  the  transformation 


z(0  =  ^(0+I 

k=\'- 


l,/'^-'^-''B,uiT)dT 


Then  following  the  arguments  in  Appendix  D,  it  is  straightforward  to  show  that 


z{t)  =  Az{t)  +  Bu{t)  +  f{x,t) 


I 


-Ah, 


where  B  =  B+  Y,Bj^e      '' .  Of  course,  (5.85)  is  of  the  same  form  as  (5.4). 

k=\ 


(5.83) 


(5.84) 


(5.85) 
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5.8  Illustrative  Example 


Consider  the  uncertain  input  delay  system  described  by  (5.1)  where 


A  = 


0 
0 


1 

0  " 

'0' 

"0" 

1.1 

0.51 

,     B  = 

1 

.     B,= 

0 

0 

0.8  _ 

_1_ 

_1_ 

and   h  =  0.5 


The  uncertainty  is  constrained  to  be  such  that  \f{x,t)\<P\xit)\  where   j3  =  0.25.     In 
particular  the  following  uncertainty  is  chosen  for  the  purpose  of  simulation. 


f{x{t),t)  =  {0.2 +  smt) 


0.1       0.03    -0.02 
-0.08    0.05      0.0 
-0.03    0.12      0.15 


X2{t) 
^3(0 


The  initial  condition  is  specified  as  x(0)  =  [11  1]   . 

The  design  objective  is  to  asymptotically  stabilize  the  uncertain  system.    The 

homogeneous  part  of  the  system  is  unstable  because  A  has  a  positive  eigenvalue  A  =  0.8. 

After   the    transformation    (5.3),    the    system   equation    (5.4)    is    realized    with 
B  =  [0.0975    0.7147    1.6703]^ and  it  is  verified  that  rank[A,5']  =  3.  The  transformation 

matrix  T  is  given  by 


r=: 


1     0 

-0.0583' 

0     1 

-0.4279 

0    0 

1 

,  TAT'^  = 


-2       1 
0     -1.1 
0        0 


0.2645  ■ 
-0.3030 
0.8 


and  from  equation  (5.7) ,  Bj  =  1.6703 .  The  poles  in  sliding  mode  for  the  nominal  system 
are  assigned  at  -1  and  -2.  This  is  achieved  for  K  =  [0  0.3301]  .  From  (5.19),  it  follows 
that    C  =  [0    0.3301    0.8588].   The  controller  parameters  are  chosen  to  be  k  =  2  and 

e  =  10"  .  Finally,  we  set  a  =  0.5.  To  guarantee  asymptotic  stability  in  the  sliding  mode, 
the  Lyapunov  matrix  is  chosen  as  P  =  0.245/2x2  and  this  satisfies  (5.46)  (In  the  stability 

condition  (5.46),  LHS  =  0.6337  and  RHS  =  0.8231)  . 

The  simulation  results  are  presented  in  Figures  5.1-5.4.   Figure  5.3  shows  the  state 
norm  of  the  transformed  system  (5.4)  as  a  function  of  time,  confirming  the  asymptotic 
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stability  of  the  transformed  system.  Figure  5.1  shows  that  the  original  system  also  is 
asymptotically  stable  as  predicted  by  the  theoiy.  The  switching  function  and  its  evolution 
towards  the  sliding  mode  is  shown  in  Figure  5.2.  The  sliding  mode  is  reached  rather 
quickly  at  around  t  =  0.909.  Figure  5.4,  shows  the  evolution  of  controller  output  and 
some  chattering  is  observed  in  the  sliding  mode  (as  the  discontinuous  control  is  given  by 
d{t)  -  j5  ||C||  \x{t)\  +  e).    The  level  of  chattering  continuously  decreases  to  e  as  the  norm  of 

the  states  decay  to  zero.  Perturbation  compensation  will  further  minimize  this  problem. 

5.8   Conclusion 

This  chapter  deals  with  the  robust  stabilization  for  a  class  of  input  delay  systems  with 
uncertain  states  via  sliding  mode  control  approach.  By  a  simple  transformation,  the 
original  system  is  converted  into  a  system  without  any  input  delay  and  this  system  is 
controlled.  The  asymptotic  stability  of  the  transformed  system  guarantees  the  asymptotic 
stability  of  the  original  system.  By  a  proper  choice  of  the  switching  function  and  the 
controller  parameters,  the  reaching  condition  and  the  asymptotic  stability  are  assured.  The 
problem  of  chattering  is  discussed  and  promising  ways  to  alleviate  this  problem  are 
discussed.  Further  work  will  focus  on  applications  to  systems  with  input  saturation  and 
uncertainty  in  the  time  delay.  It  should  also  be  interesting  to  compare  the  performance  of 
the  sliding  mode  approach  with  other  conventional  methods  for  time-delay  compensation. 

Nomenclature  for  Chapter  5 


A 

State  matrix 

(dimension  nxn) 

B 

Control  matrix 

(dimension  nxm) 

Bd 

Delayed  control  matrix 

(dimension  nxm) 

C 

Switching  function  matrix 

(dimension  mxn) 

d{t) 

Controller  parameter 

(dimension  1x1) 

fix{t),t) 

State  dependent  perturbation 

(dimension  n  x  1) 
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h 
K 

s{t) 
T 

u{t) 

x{t) 

z{t) 

w{t) 

w{t) 

Wj{t) 

Vit) 

P 

P 


Delay 

mx{n-ni)  matrix 
Switching  function 
Transformation  matrix 
Control  input 
State 

State  of  the  transformed  system 
Transformed  state 
Transformed  state 
Transformed  state 
Lyapunov  Function 
Lyapunov  function  matrix 
Lyapunov  function  matrix 
Scalar/vector  norm 
Matrix  norm 


(dimension  1x1) 

(dimension  m  x  1) 
(dimension  n  xn) 
(dimension  mxl) 
(dimension  nxl) 
(dimension  nxl) 
(dimension  nxl) 
(dimension  {n  -  m)  x  1) 
(dimension  mxl) 
(dimension  1x1) 
(dimension  {n-m)x{n-m)) 
(dimension  {n  -  m)x{n-  m)) 
(dimension  1x1) 
(dimension  1x1) 


Greek  Letters 

a 

controller  tuning  parameter 

P 

Bound  on  state  dependent  uncertainty 

3 

Razhamikhin  constant 

s 

controller  tuning  parameter 

e 

time  index 

m 

Initial  condition 
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§0.4 

0.2 

0 

-0.2 
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z1 



z2 
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A 

~  '^V 

-    \  _  _  r~^'"-'  - 

0 


1 


2        3        4        5 

t  (sec) 
Figure  2.  Evolution  of  the  transformed 
states 
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t  (sec) 

Figure  5.3.    Evolution  of  the  switching    Figure  5.4  Evolution  of  the  control  input 
function 


CHAPTER  6 
FUTURE  WORK 


6.1     Robust  Stabilization  of  Uncertain  Input-Delay  Systems-A  Linear  Matrix 

Inequality  Approach 


The  problem  of  robustly  stabilizing  an  uncertain  input-delay  system  via  sliding  mode 
control  was  considered  in  Chapter  5.  In  this  section,  we  consider  a  system  with  time 
varying  input-delay  and  uncertain  states.  Then,  we  consider  the  following  questions  of 
interest:  (1)  What  is  the  maximum  delay  that  the  system  can  tolerate  and  still  be  robustly 
stable?,  (2)  What  is  the  corresponding  stabilizing  state  feedback?,  (3)  Given  a  constraint 
on  the  controller,  is  it  feasible  to  stabilize  the  system  and  (4)  What  is  the  maximum  L2 

gain  for  the  system?.  We  approach  this  problem  via  Linear  Matrix  Inequalities  (LMI) 
approach.  The  work  in  this  section  is  motivated  by  and  in  fact  is  an  extension  of  the 
previous  works  of  Niculescu  et  al,  [1997]. 

Preliminaries  and  problem  formulation 

Consider  the  following  uncertain  time  delay  state  space  system 

x{t)  =  Ax{t)  +  Bou{t)  +  Bju{t  -  Tit))  +  f{x{t),t)  (6. 1.  la) 

Where  the  uncertainty  is  bounded  by 

||/(x(0,0||</3ol|x(0||,  (6.i.b) 

the  time  delay  is  in  the  range 

0<t(0<t'  (6.1.1c) 

and  the  controller  is  given  by  the  state  feed  back 
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uit)  =  Kxit)  (6.1.2a) 

Given  this  system,  the  first  objective  is  to  find  a  stabilizing  static  feedback  controller 

gain,  K.  Another  goal  is  find  the  maximum  value  of  T   for  which  the  close  loop  system 

is  stable.  In  addition,  the  problem  can  be  formulated  taking  into  consideration  constraints 

on  the  control  effort,  that  is 


max||M(0||<//,    fi>0  (g_^_2b) 


t>0 

Assumption  6.1.1. 

The  pair  iA,BQ  +  Bi)  is  stabihzable. 

Assumption  1,  implies  the  closed  loop  system,  which  is  free  of  delays,  can  be  stabilized. 

Controller  Design 

Let  x{t),  t  >  Obe  the  solution  of  (6.1.1).  Since  x{t)  is  continuously  differentiable  for 
f  >  0,  the  following  equation  holds  (Hale  and  Lunel,  1993) 

0 

x(r  -  T(0) = a:(o  -  jxit+eyie 
-no 

0 

=  jc(0  -    jAxit  +  d)  +  BgRxit  +  d)  +  BjKxit  -  t(0  +  6)+  ,^  j  2) 

-T(0 

f(xit+e))de,  t>T 

Substituting  (6.1.3)  into  the  (6.1.2a),  to  obtain  the  delayed  input,  and  then  substituting 
this  into  the  system  equation  (6.1.1a),  letting  5  =  Bq+Bj  gives 

0 

x{t)  =  (A  +  BK)xit)  -  BjK   jAx{t  +  0)  +  B^Kxit  +  6)  + 

-^C)  (6.1.4) 

BiKx{t  -  T(t  +  e)  +  e)+ fit + d)d6 + fit) 

Which  is  valid  for  arbitrary  initial  data  in  the  range   -2t  ,0  .  It  can  be  shown  that  if  the 

zero  solution  of  (6.1.4)  is  asymptotically  stable  then  the  zero  solution  of  (6.1.1a)  is  also 
asymptotically  stable  [Hale  and  Lunel,  1993]. 
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Theorem  6.1.1 

Consider  the  uncertain  input-delay  system  given  by  (6.1.1a-6.1.1c)  with  known 

initial  condition  x(0)  under  the  control  action  given  by  (6.1.2a-6.1.2b)  with  tuning 
parameter  kQ(P). 

If  there  exists  symmetric  and  positive-definite  matrices  g  e  i*^"^"  and  iV  g  5^"^",  a 
matrix  W  e  5^™^",  and  positive  scalars  {^,-,  /  =  l,---,5}  and  {//,-,  /  =  1,2}  such  that  the 

following  LMI's,  (6.1.5-6.1.12),  hold: 

1  r 


TdT 


QA'  +AQ  +  BW+W'B'  +  p^I  +  N 

[/5,+l52+P,lko(P)Po]  +  p4] 

w^bJ 


+ 


w"^b{^ 


AW      I3,W 


— Q      0 


-P3Q 


<o 


(6.1.5) 


-P,Q  +  AQA^<0 


P2 


-Q 
0 

0 
0 

+ 

'    0 

W^B^^ 

B^W 

-Q. 

'-N 

_Q 

Q     ' 

<0 

<0 


(6.1.6) 
(6.1.7) 

(6.1.8) 


^^I<Q<^MP)I 


Pa 


'  -1 
.^(0) 

xiof 

-Q  . 

<0 

-Q 

w''' 
-/-i_ 

<0 

'-Q    0] 

'      0          BqW' 

_  0 

0 

+ 

W^fio' 

-Q. 

<0 


(6.1.9) 
(6.1.10) 

(6.1.11) 
(6.1.12) 


then  the  system  is  closed-loop  asymptotically  stable  via  a  finite  dimensional  controller  of 
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the  form  uit)  =  Kx{t),  for  all  input  delays  satisfying  0<Tit)<T'' .  Furthermore,  the 
corresponding  input  is  given  by 

uit)  =  WQ-^xit)^Kx(t)  (6ji3) 

Proof. :  The  proof  of  Theorem  1  is  based  on  obtaining  sufficient  stability  conditions  for 
the  system  given  by  (6.1.4).  This  system  still  contains  the  specified  norm  bounded 
perturbations  (6.1.1b),  and  given  constraint  on  the  control  input  (6.1.2b).  The  stability 
requirements  are  derived  using  the  Lyapunov-Razhamikhin  approach.  By  choosing  the 
following  quadratic  Lyapunov  function 

Vix{t))  =  xitfPx{t)  (6114) 

where  P  is  symmetric,  positive  definite  and  an  element  of  5(."  ,  all  that  is  need  is  to 
show  Vix(t))<0,  for  z^>Oto  assure  asymptotic  stability  under  the  perturbations  and 
constraints.  This  choice  of  Lyapunov  function  has  several  advantages.  One  advantage  is 
that  if  there  is  an  asymptotically  stabilizing  controller,  then 

e  =  [x(t)e!Il''\xitfPxit)<l]  (6.1.15) 

is  a  holdable  ellipsoid.  This  characterization  also  makes  it  possible  to  impose  LMI  based 
conditions  to  constrain  the  control  input. 

In  the  proof  of  Theorem  6.1.1,  it  is  useful  to  impose  another  restriction  on  the 
system.  Namely,  for  any  positive  number  5  >  1 ,  the  following  holds 

Vix{0))<mxit)),       t-2t<9<t  (6116) 

This  implies  that 

\\x(e)f<5k{P)\\xitf  (6.1.17) 

This  comes  from  the  fact  that  the  Lyapunov  function,  (6.1.14),  can  be  bounded  by 

A^i„(P)||;c(0|p  <  Vixit))  <  A^,,(P)||x(0|P  (6.1.18) 
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Considering  asymptotically  stable  systems,  it  is  desired  that 

||x(0)|p  < %(0|^      <5>1,      t-2r<e<t  (6.1.19) 

so  that  the  states  decay  to  zero.  Combining  (6.1.18)  and  (6.1.19)  gives 

^miniPMsf  <  5^m^iP)Htf  (6.1.20) 

or 

Wmf  <  S^^^hitf  =  dk(P)\\x(t)f  (6_i_21) 

Where  k(P)  is  the  condition  number  of  P  and  is  equal  to  -:^ for  symmetric  P. 

This  restriction  is  standard  in  literature. 

The  time  derivative  of  the  Lyapunov  function  (6.1.14),  utilizing  the  system  given  by 
(6.1.4)  is 


\T 


\T  n  ,    nT, 


V(xit))  =  x(t)  '[iA  +  BK)'P  +  P'(A  +  BK)\x{t)  + 

-2   ]xitfPBjKAxit  +  d)  +  xit)^PBjKBoKxit  +  6)  + 

-^(')  (6.1.22) 

xitf  P{BjKfx(t  -  T{t  +  9) +  9)  +xitfPBjKf(t  +  9)de  + 

2x{t)^Pf{t) 
It  is  desired  to  bound  the  individual  components  of  the  integral  term  of  (6.1.22)  in 
order  to  factor  out  the  uncertain  time  delay  and  eliminate  the  integration.   This  can  be 
done  by  using  the  following  inequality  for  any  real  constant  /?  >  0  and  any  symmetric 
and  positive  definite  matrix  S,  [Su,  1994] 

-2u^v<^~^u^S^^u  +  j5v'^Sv  (6.1.23) 

Although  the  upper  bound  of  (6.1.23)  is  conservative,  this  conservatism  can  be  reduced 
by  optimizing  over  /3.  For  the  first  term  of  the  integral,  consider  the  following 
substitutions 


Ill 


u'^  =xitf{PBjKA) 
v-xit  +  d) 
S  =  P 


(6.1.24a-6.1.24c) 


Which  gives 


]-2x{tfiPBjKA)xit  +  d)dd 


-T(f) 


0         T 

\-2u^vde 

-T(0 

< 

-T{t) 


(6.1.25) 


\Pfh{t)^  {PBiKA)P~\PBiKA)^  xit)  +  p^x{t  +  9fPx(t  +  6) 

-no 


Using  (6.1.14)  and  (6.1.16)  this  becomes 


^  T 

j-2x{tyiPBiKA)x(t  +  0)d9 

-no 

< 

xitf    jPf\PBjKA)P~\PBiKAf  +  l5^5Pdexit)  (6.1.26) 

-no 


Tit)x(tf[pf\PB^KA)p-\PB^KAf  +  /3,5P]x(0 

As  desired,  (6.1.26)  is  a  bound  on  the  integral  term  that  is  not  in  terms  of  delayed  states 
nor  requires  integration.  Similarly,  for  the  second  term  of  the  integral 


0 


j-2x(0 '  {PBjKBoK)x(t  +  e)de 

-no 

< 

r(t)xitf[l3f\PBiKBoK)P~\PBjKBgKf  +  I3^5p}:cit) 


(6.1.27) 
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For  the  third  term 


0 
I 

-T(0 


0  ,  s 

J  -2xitf\^p(BjKyjxit-Tit+e)+e)de 


< 


(6.1.28) 


T(r)x(0' 


P2-HP{BiKf]p-'ip{BjKf]    +p,5P 


xit) 


The  last  term  in  the  integral  is  slightly  different.  Using  (6.1.23)  it  is  bounded  by 

0  T 

l-2x{tyiPBiK)fit+d)de 

-T(0 

< 

0 

j/5f^x(tf{PBjK)p-\PBjKfxit)  +  l3^fit  +  dfPfit  +  d)dd 

-T(t) 


(6.1.29) 


Now,  the  goal  is  to  approximate  the  term  containing  the  delayed  uncertainty  with  an 
appropriate  bound.  From  (6.1.18)  and  the  original  bound  on  the  uncertainty,  (6.1.1b),  it  is 
true  that 

f{t+dfpf{t+e)<x^^^{P)f{t+dff{t+e)<pQi^,^{P)x{t+e)'x{t+e)     (6130) 

From  (6.1.20)  this  becomes 


/(r  +  0)'rp/(r  +  0)  < /3oA^,,(P)5^^naxg);,(r)T^(^) 


^min(-P) 


(6.1.31) 

Using  (6.1.18)  and  previously  mentioned  definition  of  the  condition  number  this  becomes 

i2 


\T 


f{t+6ypf{t+d)<%5 


2     /  nT 


A^i„(P)x(0'  x(0  <  ^Q8[k{P)]-x{tYPx{t)    (6_^32) 


The  condition  number  can  be  bounded  by 


k{P)<k{PQ) 


(6.1.33) 


such  that 


i2„,.nT 


f{t+dypf{t+e)<iiQ5[k{PQ)]x{typx{t)  (6^34) 

Combining  (6. 1 .34)  and  (6. 1 .29)  gives  the  bound  for  the  last  term  of  the  integral 
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T{t)x{ty 


\-2x{t)^{PBiK)f{t+e)de 

-T(0 


/33-'  iPBjK)p-^  (PBjK)^  +  p^/5od[k(Po)f  P 


(6.1.35) 


;c(0 


Combining  (6.1.26-6.1.28,  6.1.35)  and  (15)  gives 


Vixit))<xitf[MiTit),d)]xit)  +  2xitfPfit) 


(6.1.36) 


where 


M(Tit),S)  =  (A  +  BKfp  +  P^iA  +  BK)  +  T(r)[^f'  iPBjKA)P~\PBjKAf  + 

p2-'lP{BjKf)p-Hp{BjKff  +{5f\PBjK)p-\PBjKf  + 


l5f\PBjKBoK)p-\PBjKBoKf  +  A  +  /?2  +  p3pQ[k(Po)f  +  Pa 


8P 


(6.1.37) 
Asymptotic  stability  can  now  be  guarantied  using  the  Lyapunov-Razhamikhin 
approach  by  showing  that  the  time  derivative  of  the  Lyapunov  function,  (6.1.36),  is 
negative  definite  subject  to  the  uncertainty  constraint  (6.1.1b).  A  necessary  condition  for 
(6.1.36)  to  be  negative  definite  is  that  M(T(t),5)  be  negative  definite.  Using  the 
continuity  property  of  the  eigenvalues  of  M,  it  can  be  argued  that  if  MiT(t),5  =  1)  is 
negative  definite,  then  35  >  I,  sufficiently  small  such  that  Mir(t),d)  is  negative  definite. 
An  LMI  based  scheme  can  now  be  designed  to  guarantee  M(T(r),(5  =  1)  is  negative 
definite.  From  (6.1.36)  it  follows  that 


V{x(t))  <  x{tf[M(t(t),d  =  i)]x(t)  +  IxitfPfit) 
In  addition,  the  uncertainty  constraint  can  be  rewritten  as 

fitffit)-xitfipQl)x{t)<0 
From  (6.1.38),  V{x{t))  is  negative  definite  if 


xitf     fit) 


The  uncertainty  constraint  is  equivalent  to 


M    P 
P     0 


xit) 
f{t)_ 


<0 


(6.1.38) 


(6.1.39) 


(6.1.40) 
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x{tV     fit) 


-PqI  0 

0       I 


'x(ty 
fit) 


<0 


(6.1.41) 


(6.1.40-6.1.41)  can  be  combined  into  a  "single"  matrix  equation  through  the  s-procedure 
for  quadratic  functions  and  non-strict  inequalities  [see  Boyd  et  al.,  1994].  If  there  exists  a 
/?5  >  0,  then  (6.1.40-6.1.41)  are  equivalent  to  the  condition 


M  +  Ms'l        P 


<0 


(6.1.42) 


(6. 1 .42)  can  be  expanded  as 


Defining  a  matrix 


Mirit),5  =  l)  +  {5Ql5s'l  +  P,P^<0 


(6.1.43) 


j-i] 


M,:=MiTit),5  =  l)  +  PoP^'l  +  P,P'  (6_1_44) 

then  for  any  continuous  time  varying  delay,  t(0  ,  satisfying  0  <  T(r)  <  T ' ,  there  exists  an 
a3  >  0  such  that 

Ml  <  0  -»  Vixit),t)  <  -a4xit)f  (^_i_45) 

Lemma  6.1.1  gives  the  sufficient  conditions  for  guaranteeing  Mj  <  0.  Lemma  6.1.2 

gives  sufficient  conditions  to  ensure  that  the  constraint  on  the  controller,  Equation  (4b)  is 
satisfied.  Lemma  6.1.3  constrains  Q  such  that  Equation  (6.1.33)  holds. 

Lemma  6.1.1 

There  exists  a  symmetric  and  positive-definite  matrix  P  g  5^;^"^"  and  a  matrix 
2^g^mxn   g^^]^  j.j^^j.   j^^^Q^  If  1-j^ej.g  exist  the  symmetric  positive-definite  matrix 

gei^"''"  and  a  matrix  We  5^'"''"  satisfying  (6.1.5  -  6.1.10). 
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Proof:  From  (6.1.37),  Mj  is 

M^{T{t),5  =  1)  =  (A  +  BKfP  +  PiA  +  BK)  +  Pq^^^I  +  P^P^  + 


T(0 


Pf\PBiKA)P~\PBjKAy+p2~HP{BiK))p-Hp{BjK))    + 


2\„-l 


Pf^(PBjK)P~'(PBjKy  +  pf\PBjKB(,K)P  '{PBjKBoKY  + 

p,+P,+l5,Po[k(Po)f+p4]P 

(6.1.46) 
The  goal  is  to  determine  conditions  that  guarantee  Mi(T(r))  <  Oby  satisfying  Mi(t'')  <  0. 

This  condition  is  not  jointly  convex  in  P  and  K.  Fortunately,  a  change  of  variables 
results  in  an  equivalent  condition  that  is  an  LMI  [Bernussout  al.,  1989].  Multiplying 
(6.1.46)  by  Q  on  both  sides  and  substituting  P  =  <2"'  and  K  =  WQ~^  results  in 

QMiQ  =  QA^  +AQ  +  BW  +  W^b'^  +  /'\^Pf^BjWQ~^AQA^Q~^w'^Bj  + 

P2~^BjWQ'^BjWQ-^w"^bJq-^w^bJ  +  pf^  BjWQ-^w^bJ  + 

Pf'B^WQ-'BoWQ-'w^B^Q-'w^Bj  +  [^  +^2  +  PM^iPo)]^  +  P^JQ 

PoP;'q'  +  PsI 

(6.1.47) 
In  order  to  get  a  tractable  sufficient  condition,  the  following  constraints  are  imposed 
[Niculescu  et  al.,  1997] 


+ 


-I5^Q  +  AQA^<0 


-j32Q  +  BjWQ-'W^Bj<0 


(6.1.48) 


(6.1.49) 


In  addition,  the  following  constraints  are  also  imposed 

Note  that  (6.1.5  -  6.1.6,  6.1.12)  represent  (6.1.48  -  6.1.51)  in  LMI  form.    Substituting 
(6.1.48  -  6.1.51)  into  (6.1.47)  yields  the  sufficient  condition 


116 


QM^Q  <  [qA^  +AQ  +  BW  +  W^B^  +  P^I  +  n]  +  t''[3BjWQ~^w'^b'^  + 


Q 


<0 


(6.1.52) 


(6.1.1b)  represents  (6.1.52)  in  LMI  form.    Hence,  if  the  LMI's  given  by  (6.1.5  -  6.1.8, 
6.1.12)  are  satisfied,  it  is  guaranteed  that  Mj  <  0. 

Lemma  6.1.2 

Suppose  that  there  exists  matrices  We!?(.™^"  and  Q  =  Q^  >  0  e  Hi'^^^  such  that 
uit)  =  WQ~  x(t)  stabilizes  the  closed  loop  system  given  by  (6.1.1  -  6.1.3)  with  x{0) 

known,  then  the  following  two  statements  are  satisfied. 

1.  The  ellipsoid  represented  by  e  =  jx  6  ^"  x(t)  Q~  x(t)  <  1  >  is  invariant,  and 

2.  The  input  constraint  given  by  (6.1.4b)  is  enforced  at  all  times  r>0  if  the 
following  LMI's  are  satisfied 


-1 

x{Of' 

x(0) 

-Q  . 

-Q 

w^' 

w 

-^h_ 

<0 


<0 


(6.1.53) 


(6.1.54) 


Proof:  An  outline  of  the  proof  of  Lemma  6.1.1  is  given  in  Boyd  et  al,  [1994].  The 
details  are  given  here  for  clarity.  Since  the  feedback  u{t)  -  WQ~  x(t)  stabilizes  the 
closed  loop  system  the  following  is  true 

Vixit),t)  <  y(jc(0),0)  =»  xiOfQ'^xiO)  <  1  ^  x{tfQ~^xit)  <  1 
Hence,  e  is  invariant  if  the  LMI  given  by  (6.1.10)  holds.    In  order  to  constrain  the 
controller,  consider 


max||M(0||  =  rnax  WQ    x(t) 
t>0  t>Q 


< 


^ 


-  max 
e 

<  max 

e 


WQ~'x 


x^Q'^W^WQ'^x 


Ml 
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1/2, 


Using  the  transformation  i  =  jg      x  it  is  true  that 


X '  it)x(t)  <  1 
which  comes  from  the  invariance  of  e .  Using  this  transformation  gives 


(6.1.55) 


Hence, 


max||M(f)||  <  max 
f>0  e 


<A„ 


i-l/2wTTjr^-l/2 


i^(oe"'^^w^we"'^^x(o 


1/2 


g-l/2-^T^g-l/2 


g-l/2^T^g-l/2 


72 


x^(Oia) 


(6.1.56) 


1/2 


Q-"^W'WQ^"'  -ii'l<Q    ^    max||H(0||  < //,    A^  >  0       ^g^^^^) 


f>0 


This  condition  can  be  recast  as  an  LMI 


\w'^W-QI<0 
1^' 


(6.1.58^ 


Equation  (49)  is  a  generalized  eigenvalue  problem  if  jj.  is  not  known.  For  this  problem  it 
is  known,  and  (6.1.58)  can  be  rewritten  as 


-Q     W 


<0 


(6.1.59) 


Lemma  6.1.3 

The  symmetric  positive  definite  matrix  Q  has  condition  number  k(Q)<k(PQ)  if 

there  exists  jj.  >  Osuch  that 

fa<Q<^k{P,)I  (6  160) 

Proof:    Equation(6.1.60)  with  /:(/o)  unknown  is  an  eigenvalue  problem  [Boyd  et  ai, 
1994].  In  this  case  k(PQ)  is  a  tuning  parameter.  From  Equation  (6.1.60)  it  follows  that 


jjx^x  <  x^Qx  <  fik{PQ)x^x 


(6.1.61) 


In  addition 


X^,^(Q)x^x<x'^Qx  <  X^,^iQ)x''x 


(6.1.62) 
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Equation  (6.1.61  -  6.1.62)  imply 

K.AQ)^^k{Po)  (6.1.63) 

Therefore 

k{Q)  =  AmaxM  <  ^(p  ) 


The  goal  is  to  find  a  static  state  feedback  controller.  Note  that  if  the  static  controller 
is  not  feasible  through  these  LMI's,  then  LMIs  that  would  result  through  a  dynamic 
controller  will  also  be  infeasible. 

L2  Gain 

This  section  considers  the  problem  of  finding  a  state  feedback  gain,  K,  such  that  the 
L2  gain 

II  II  I'l^il 

sup  114=   sup  r^  (6_i_65) 

Il/Il2=l  11/112^0  ll/ll  ^  ^ 

of  the  closed  loop  system  is  less  than  a  specified  number  / .  Assuming  that  the  bound  on 
the  delay,  r  ,  is  known.  Again,  the  system  considered  here  is 

x{t)  =  Ax(t)  +  Bguit)  +  Bju{t  -  T(0)  +  f(x(t),  t) 
with  the  output  given  by 

z{t)^cx{t)  (g_^_g6) 

Theorem  6.1.2 

The  uncertain  input-delay  system  characterized  by  (6.1.1  -  6.1.4)  with  known  T  ,  has 
an  L2  gain  less  than  or  equal  to  7 ,  if  the  following  LMI's  are  satisfied 


QA^  +  AQ  +  BW  +  W'^  b'^  +  P^I  +  N  +  c'^c  +  N2 


+ 


W"^B^^ 


Q 


Q 

-1 


r'Po 


■ 

PiW 

AW 

0 

0 

<0 


<0 


In  addition  the  LMI's  given  by  (6.1.6  -  6.1.9,  6.1.12)  must  also  be  satisfied. 
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(6.1.67) 


(6.1.68) 


Proof:    These  LMI's  come  from  arguments  similar  to  those  used  in  proving  Theorem 
6.1.1.  Using  the  Lyapunov-Razhamikhin  approach  results  in  the  following  conditions 


V{x{t))  +  z'^z-7'-ff<0 


.T„T, 


(6.1.69) 

V{x{t))  +  x'c'cx-  y^PqX ' X  < 0  (g^  7Q^ 

As  mentioned,  the  remaining  steps  are  similar  to  those  show  for  Theorem  6.1.1,  and  are 
therefore  omitted. 

Future  work 

The  conditions  derived  in  this  section  are  sufficient-only.  They  give  a  systematic 
framework  for  getting  a  feasibility  answer  for  the  host  of  questions  raised.  However,  it  is 
imperative  to  try  numerical  experiments  to  see  the  relative  advantages  of  this  method 
versus  the  others  in  the  literature.  This  will  form  part  of  the  future  work. 

6.2     Automatic  Covariance  Resetting  for  On-Line  Recursive  Identification 


The  most  well  known  on-line  identification  technique  is  probably  the  recursive  least- 
squares  method  (RLS)  which  has  found  wide  application  for  the  parameter  estimation  of 
physical  systems.    The  reasons  for  its  wide  spread  use  are  its  ease  of  implementation. 
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relatively  low  computational  load,  and  well-known  statistical  properties.  However, 
problems  arise  when  RLS  technique  is  applied  to  time-varying  systems.  The  gain  of  the 
algorithm  decreases  monotonically  towards  zero  as  a  function  of  time.  Hence,  after  some 
time  the  algorithm  loses  its  capability  for  further  tracking  any  parameter  changes.  This  is 
a  well  known  phenomena  and  the  approaches  proposed  for  addressing  this  limitation  can 
be  classified  in  two  major  groups  [see,  Ljung,  1987],  namely,  techniques  that  use 
selective-forgetting  algorithms,  and  techniques  based  on  covariance  management. 

The  first  group  of  techniques  use  the  concept  of  a  forgetting  factor  that  selectively 
weights  the  data  depending  on  its  information  content  or  time  of  acquisition.  The  main 
utility  of  the  forgetting  factor  is  to  prevent  excessive  gain  decay,  and  in  addition,  to 
automatically  readjust  the  gain  to  a  high  value  whenever  parameter  changes  are  detected. 
There  are  many  techniques  in  the  literature.  For  extensive  review  and  further  discussions, 
see  Lung  and  Gunnarsson  [1990],  Ramombason  et  al.,  [1992],  and  standard  texts  such  as 
Goodwin  and  Sin  [1984].  The  techniques  range  from  simple  exponential  forgetting 
[Goodwin  and  Sin,  1984],  to  elaborate  schemes  such  as  the  directional  forgetting  strategy 
of  Kulhavy  and  Karny  [1984].  Particularly  good  performance  has  been  observed  by  the 
directional-forgetting  algorithm  of  Canudas  de  Wit  and  Carillo  [1990],  which  features  the 
concept  of  a  dead  zone  that  can  be  specified  as  a  function  of  the  measurement-noise 
variance.  Also  see.  Woo  et  al.,  [1995]  where,  a  directional  forgetting  matrix  is  employed. 
Unfortunately,  the  performance  of  estimation  algorithms  with  forgetting  factors  depends 
strongly  on  the  adequate  adjustment  of  a  number  of  heuristic  parameters. 

The  second  group  of  techniques  are  those  that  directly  manipulate  the  parameter- 
covariance  matrix  of  the  RLS  algorithm.  Such  techniques  are  effective  because  the  gain 
of  the  algorithm  and  the  parameter-covariance  matrix  are  directly  related.  Covariance 
management  is  typically  done  ad  hoc,  with  the  assistance  of  arbitrary  premises.  The  reset 
limit  as  well  as  the  new  covariance  value  must  be  arbitrarily  specified  by  the  user.  Often 
these  specifications  are  made  invoking  strong  heuristics,  without  regard  for  the  level  of 
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process  noise  or  type  of  parameter- variation  dynamics.  The  arbitrary  modifications  made 
to  the  covariance  matrix  make  it  lose  its  meaning  as  a  statistical  covariance.  Also,  these 
choices  are  not  obvious  in  practice. 

The  high  level  of  heuristic  considerations  involved  in  selective  forgetting  and 
covariance  management  is  undesirable.  Although  these  two  alternative  approaches  for 
modifying  the  RLS  algorithm  can  most  often  be  tuned  for  successful  performance,  they 
involve  heuristics  or  choices  which  are  not  obvious  to  make  in  practice.  This  is 
particularly  relevant  in  the  context  of  adaptive  controllers,  where  the  estimation  algorithm 
should  not  add  to  the  burden  of  the  work  required  to  tune  the  operation  of  the  controller. 
Therefore,  there  is  a  strong  motivation  to  minimize  the  number  of  ad  hoc  adjustable 
parameters  in  RLS  estimation  algorithms  designed  for  tracking  time-varying  systems. 

In  order  to  implement  an  on-line  identification  strategy  along  with  the  controller  in  a 
process,  it  is  imperative  to  minimize  the  number  of  ad-hoc  parameters.  With  this  goal  in 
mind,  studies  have  been  done  with  an  improved  RLS  algorithm  for  time-varying  systems 
that  automatically  adjusts  its  gain  matrix  in  order  to  track  parameter  changes  Isaksson 
[1987],  [Ramombason  et  ai,  1992].  The  new  algorithm  maintains  its  capability  for 
tracking  parameter  changes  by  estimating  on-line  the  parameter-covariance  matrix  as  well 
as  the  measurement-noise  variance.  Therefore,  the  scheme  has  an  embedded  covariance- 
estimation  algorithm  that  in  principle  can  characterize  the  variance  of  process  parameter. 

However,  there  are  still  some  open  issues.  For  example,  the  performance  of  the 
algorithm  is  suspect  when  there  is  process  noise.  Also,  there  are  no  guarantees  for 
convergence.  In  the  ensuing,  the  standard  problem  formulation  is  given  and  an 
interesting  insight  is  presented.  The  problem  is  the  same  as  considered  in  Isaksson 
[1987]  and  Ramombason  et  ai,  [1992] 

Let  the  process  be  characterized  by  the  time-varying  linear  regression 

yit)  =  fit)9{t)  +  e{t)  (62.1) 

where  the  output  y{t)  and  noise  e{t)  are  scalar  signals,  and  0  (r)G9t"  and  0(r)e9t    are 
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the  regression  and  parameter  vector,  respectively.  The  noise  e{t)  is  assumed  to  be  a 
zero-mean  Gaussian  white-noise  sequence  with  vaiiance 

Ee{tf  =  R{t)  (52.2) 

where  E  is  the  expectation  operator.  The  process  parameters  in  (1)  may  change  with 
time.  For  example,  many  physical  processes  undergo  step-wise  parameter  changes  at 
random  times.  Other  patterns  of  parameter  evolution  include  ramps,  sinusoidal  forms, 
progressive  drifts,  etc.  Often,  parameter  variations  are  modeled  as  a  Markovian  process, 
or  more  specifically,  as  the  random  walk 

dit)  =  e{t-l)  +  wit)  (62.3) 

where  w{t)  e  91  is  a  parameter-innovations  vector  whose  elements  are  zero-mean  white 
Gaussian  noise  sequences  such  that 


Ew{t)w''{t)  =  Q{t)  (62.4) 


and 


Ee{t)w{t)  (62.5) 

where  the  parameter-covariance  matrix  Q{t)  e  9t"^"  is  positive  semi-definite.  Although 
the  random-walk  model  (6.2.3)  may  not  be  a  realistic  representation  of  the  type  of 
parameter  drifts  experienced  by  a  particular  process.  The  random- walk  paradigm  is 
nevertheless  useful  for  theoretical  purposes. 

Least-squares  parameter  estimation  techniques  seek  to  find  the  estimate  0{t)  that 
minimizes  the  quadratic  functional 

m)  =  Ei{yii)-cp\i)e{i)) 

i=V  '  (6.2.6) 

The  optimal  solution  to  (6.2.6),  under  the  assumption  of  parameter  drift  according  to  the 

random-walk  (6.2.3),  is  given  by  the  BLUE  (Best  Linear  Unbiased  Estimator),  a  Kalman 

filter  that  strikes  the  best  compromise  between  parameter  tracking  and  noise  tolerance 
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[cf.  Ljung  and  Gunnarson,  1990]: 

For  the  case  of  time-varying  systems,  the  following  modified  formulation  of  the 
BLUE  is  used: 


e{t)  =  e{t-l)  +  L{t)e{t)  (6.2.7) 


,r 


£{t)  =  yit)-(p'{t)d{t)  (g_2.8) 

^  P{t-l)(p{t) 

^^     R{t)  +  (p^{t)P{t-\)(p{t)  (^■2-9) 


P{t-\)cp{t)(p''{t)P{t-\) 
R{t)  +  (p\t)P{t-\)(p(t) 


P{t)  -  P{t  -  1)  -    A...    -JtAP<f-^^rn(A  ^  ^^'^  (6-2. 10) 


where  the  alternative  covariance  matrix  Q{t)  and  variance  function  R{t),  obtained  by 
various  methods,  are  used  instead  of  their  tme  unknown  counterparts. 

The  approaches  in  Ramombason  et  al.,  [1992]  and  Isaksson  [1987]  consist  of 
estimating  Q{t)  and  R{t).  Towards  this,  consider  the  following 

£{t)  =  (p^{t)w{t)  +  e{t)  (6.2.11) 

where 

w{t)'=0{t)-d{t-\)  (6.2.12) 

Under  the  assumptions  that  the  random  sequence  is  ergodic,  the  regressor  vector  (p{t)  is 
deterministic  and  Ew{t)\v^ {t)  =  Eyv{t)w'^ {t),i\.  is  possible  to  arrive  at 

£{tf=(p'^m{t)(p{t)  +  R[t)  (6.2.13) 

Since  the  above  can  be  rearranged  in  the  form  of  a  hnear  regression,  it  is  then  possible  to 
estimate  the  elements  of  Q{t)  and  R{t).  This  is  a  constrained  optimization  problem, 
since  Q{t)  has  to  be  positive  semi-definite.  The  convergence  properties  of  the  available 
algorithms  have  not  been  proven.  Now,  we  find  that  a  simple  substitution. 


Q{t)  = 


(6.2.14) 


124 

satisfies  (6.2.13).  In  other  words,  this  solution  parametrizes  all  possible  diagonal 
solutions  in  terms  of  a  single  parameter  R{t).  If  R{t)  can  be  estimated  on-line  (or)  an 
approximate  value  is  known  this  then  can  be  used  to  estimate  an  approximate  value  of 
Q{t)  without  any  need  for  constrained  optimization.  Note  that  for  any  value  of  chosen 
R{t),  the  resulting  Q{t)  is  automatically  positive  semi-definite.  This  simple,  yet  elegant 
solution  appears  to  be  new  and  deserves  to  be  investigated  under  various  scenarios. 

6.3     Sliding  Mode  Control  of  Bilinear  Systems 

Another  idea  is  to  use  the  VSC  approach  for  the  stabilization  of  uncertain  bilinear 
systems.  Bilinear  systems  are  a  special  class  of  nonlinear  systems.  They  are  interesting 
as  well  as  important  since  various  processes  in  science  and  engineering  can  be  modeled 
as  bilinear  system.  Past  work  in  the  literature  has  included  optimal  control  strategies  for 
bilinear  systems  and  the  intent  of  the  future  research  is  to  compare  the  performance  of  the 
sliding  mode  control  technique  against  the  existing  approaches. 

Nomenclature  for  Chapter  6 


A 

Bo 

Bi 

e{t) 
E 

f(x{t),t) 
T(r) 

T 
K 


State  matrix 
Control  matrix 

Delayed  control  matrix 

Noise 

Expectation  operator 

State  dependent  perturbation 

Variable  delay 

Upper  bound  on  variable  delay 

State  feedback  matrix 


(dimension  nxn) 
(dimension  nxin) 

(dimension  nxm) 
(dimension  1x1) 

(dimension  nxl) 
(dimension  1x1) 
(dimension  1x1) 
(dimension  mxn) 
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KiP) 

Tuning  parameter                               ( 

^dimension  1x1) 

u{t) 

Control  input                                       ( 

; dimension  mxl) 

x{t) 

State                                                   ( 

; dimension  nxl) 

V{t) 

Lyapunov  function                             ( 

;dimension  1x1) 

P 

Lyapunov  function  matrix                  ( 

^dimension  nxn) 

Pit) 

Parameter  covariance  matrix              ( 

; dimension  nxn) 

N 

Positive-definite  matrix  variable         ( 

^dimension  nxn) 

Q 

Positive-definite  matrix  variable        ( 

^dimension  nXn) 

m 

Parameter-covariance  matrix              ( 

^dimension  nxn) 

m 

Estimated  covariance  matrix              ( 

^dimension  nXn) 

Rit) 

Variance                                            ( 

;dimension  1x1) 

R(t) 

Estimated  variance                             ( 

;dimension  1x1) 

W 

Matrix  variable                                   ( 

^dimension  mxn) 

wit) 

Parameter  innovations  matrix             ( 

; dimension  nxl) 

x(0) 

System  initial  condition                      ( 

; dimension  nxi) 

Greek  Letters 

^ 

Bound  on  state  dependent  uncertainty 

A 

Scalar  variable 

m 

Model  prediction  error 

jU 

Upper  bound  on  controller  effort 

Ml 

Scalar  variable 

5 

Razhamikhin  constant 

e 

time  index 

Bit) 

Parameter  vector 

¥it) 

Regression  vector 

APPENDIX  A 

In  this  section,  we  give  a  detailed  proof  of  Lemma  3.2.  However,  before 
proceeding,  we  review  certain  fundamental  aspects  of  tlie  M-A  paradigm.  For  complete 
details,  the  reader  is  referred  to  Chapter  10  in  Zhou  et  al,  [1996].  A  standard  tool  in 
robustnesss  analysis  is  to  exploit  the  concept  of  Linear  Fractional  Transformation  (LFT). 
The  following  diagram  shows  an  upper  LFT,  F^^{M,A) , 


Figure  A.  1  The  classical  M-A  structure 


where  A  is  the  norm  bounded  uncertainty  and  ^  is  a  complex  matrix  that  can  be 
partitioned  as 


M  = 


12 
M21       M22 


This  results  in. 


^z^ 


yyj 


Mil     ^12' 

^21  ^22 


ry,\ 


\fj 


From  the  above  diagram,  it  is  straightforward  to  show  that 

y  =  F,{M,A)  =  \M22  +  M2iA{l-MnAy'M,2 


\-i 
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In  the  above  equation,  M22  represents  the  nominal  mapping,  and  is  perturbed  by  A,  while 
Mjj,  Mj2,  and  M21  reflect  the  way  in  which  the  perturbations  affect  the  nominal  map. 
The  control  theoretic  insight  is  gained  by  obeserving  that  the  F^^{M,A)  is,  in  fact,  the 

closed  loop  transfer  function  between  the  input  and  the  output.  Obviously,  for  the  upper 
LFT  to  be  well  defined,  the  [l -M^^A)  must  be  invertible.   Naturally,  this  leads  to  the 

determinantal  condition  det(/  -MnA^j^O. 

In  the  following  analysis,  with  little  abuse  of  notation  we  will  refer  to  Mjj  as  M . 

Our  task  is  to  find  M,  that  is  the  transfer  function  between  z  and  w .   The  SISO  transfer 

function  (3.1)  with  affme  uncertainty  can  be  written  as  shown  in  (3.3). 

m 
n{jO),0)+  ln,^(jCO)q,^(jCO) 

8ij(0,q)  = ^ 

d{j(O,0)+  Y,dj^{jO))q^{j(0) 

k=\ 

The  aim  is  to  transform  the  system  (3.3)  under  unity  feedback  in  to  a  M-A  structure  by 
extracting  the  pai^ametric  uncertainties  and  including  them  in  a  diagonal  the  matrix 

4  =  diag{(?,,g2,--,  ^m} 
For  ease  of  notation,  henceforth,  we  will  drop  the  arguments  and  denote,  Uq  =  nijco,0) 
and  cIq  =  d{jco,0).  With  these  simplifications,  we  can  write 


y{s) ^ 


'>      do+ld,q, 


Cross  multiplication  yields, 


«0  ^=1 


k=\ 


Hj.       dj, 

-^u — -y 

Jg  dQ 


Conceptually,  this  can  be  represented  by  Figure  A. 2.  While  for  ease  of  exposition,  Figure 
A. 2  shows  only  one  uncertainty,  multiple  uncertainties  can  be  trivially  added  to  the 
diagram.  Figure  A. 3  shows  the  SISO  system  under  unity  feedback.  Furthermore,  in  order 
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to  extract  the  effect  of  uncertainty,  we  label  the  input  and  output  of  the  uncertainties  as  z^ 
and  w^,  respectively.  Then,  by  following  Figure  A. 3,  yields 

"0       k=\ 


( 


\ 


V"o  +^^0  7 


m 
k=\ 


Uq 


y^riQ+dQ  J 


Zk={f-y)^--Ty 


(Iq      dQ 


HL- 


i'h+dk) 


Substituting  for  y  reults  in 


Zk  = 


dQ      {nQ+dQ)dQ 


r  — 


7 — TTT  ^^^ 


Defining  M,/,  as  the  transfer  function  from  the  output  of  q]^  to  the  input  of  qi ,  it  is 
straightforward  to  conclude  that 

Using  definitions  (3.5)  and  (3.6)  leads  to  the  final  result 

Pkis) 


Mikis)  = 


Poi^) 


,     i  =  \,2,---,m;     k  =  l,2,- 
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^k 

1 

"0 

Figure  A.2  Block  diagram  of  the  open-Loop  SISO  Transfer  function 


Zk 

Ik 

w^ 

^ 

dt 
do 

r 

"/t 

1 

V     , 

'^r^ 

1 

do 

-< 

+ 

r  + 

~^\ 

/ 

"0 

J 

y 

-y 

Figure  A.3  Block  diagram  of  the  closed-Loop  SISO  Transfer  function 


APPENDIX  B 

Proof  of  Theorem  4.2 

From  (4.15)  we  have, 

xit)  =  Axit)  +  Adxit-h)  +  f{x{t),t)  +  f^l{xit-h),t)  ^g  j^ 

In  order  to  derive  sufficiency  conditions,  thie  perturbations  can  be  replaced  by  their 

upper  bounds  as  in  equation  (4.2).  Then,  following  steps  similar  to  those  in  the  proof  of 

Theorem  4.1  in  [Mori  et  ai,  1981]  and  using  the  Comparison  Theorem  [Hale  and  Lunel, 

1993],  it  can  be  shown  that  system  (4.12)  is  asymptotically  stable  if  (4.13)  holds.   The 

details  are  given  below: 

The  solution  x{t),\/t  >  0  can  be  written  as 

_  t        _ 

xit)  =  exp(Ar)jc(0)  +  j  exp[A{t  -  s)y{s)ds  (B  .2) 

0 

where  r{s)  is  given  by 

r{s)^Adxis-h)  +  f{xis),s)  +  f^{x(s-h),s)  (B.3) 

Taking  the  norm  on  both  sides  of  (B.2)  and  using  (B.3)  and  the  following  inequality 

exppfjl  <  exp(/i(A)r),  Vf>0  [Vidyasagar,  1993] 
yields 

\x{t)\  <  cxpUA)t)x{0)  +  lexp[ii(A){t  -  s))\r{s)\ds  (B.4) 

0 

Put  v{t)  =  \x{t)\,  we  can  rewrite  (B.4)  as 

v(0  <  exp[iLi(A)t)v{0)  +  j  Qxp[fl[A){t  -  s)jr{s)\ds  (B  .5) 
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with  v{t)  =  |(/>(r)|,  - /j  <  r  <  0,  and  v(0)  =  vq. 

Now,  consider  the  integral  equation  with  the  deviating  argument 

z{t)  =  exp(/i(I)?)v(0  +  jexp(//(A)(f -  s))  ([|a^|+  kAvis -  h)\  +  k\v{s)\]  ds{B.6) 

with  zit)  =  \(^{t)\,  -h<t<0. 

Of  course,  (B.6)  is  the  solution  of  the  scalar  differential  equation 

z{t)  =  pi{A)z{t)  +  [\Ad\+Ki]zit-h)  +  kz(t),  ziO)  =  vo  (B.7) 

The  solution  of  the  above  differential  equation  uniquely  exists  and  is  asymptotically  stable 
if  and  only  if 


-[li[A)  +  k)>[\Ad\+k 
fi[A)  +  k  +  \Ad\+k^  <0 


(B.8) 


If  the  condition  (B.8)  is  satisfied,  using  the  comparison  theorem  with  m=l  (See  Appendix 
C),  we  can  write  K^)I  =  ^W^^W. 


Hence,  lim  z{t)  ->  0  ^  \x{t)\  -^  0 . 

t->oo 


APPENDIX  C 
Comparison  Theorem 

Assume  a  vector  valued  function  m{t,s,z)  :JxJx  R'"  -^  R"\J  :=  [tQ,^]  has  the 

following  property. 
For  any  fixed  t,  s, 

Zi  <  Z2  =^  m{t,s,zi)  <  m{t,s,Z2) 
Let  z(t)  be  the  solution  to  the  inequality 

z{t)  <  z{0)  +  jm{t,s,z{s  -  h))ds 

'a 
t 

Then  the  maximal  solution  r{t)  of  w{t)  <  z(0)  +  \m[t,s,w{s  -  h))ds  satisfies 


Z{t)<r{t),t>t^ 
For  a  proof,  see  [Laksmikantham  and  Leela,  1969]. 
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APPENDIX  D 


Consider  the  system 

x{t)  =  Ax{t)  +  Bu{t)  +  B^iu{t  -h)  +  f{x{t),t) 
Introduce  the  transformation 

z{t)  =  xit)  +  \l,^e'^^'-''-'^B^iu{T)dT 
Taking  the  time  derivative  to  (D.2),  we  have 

m  =  x(t)  +  j^[j;_^^e^^'-"-'^B^iu(T)dT]  +  f{x{t),t] 
Using  Leibnitz  rule  for  differentiation  under  the  integral  sign,  we  have 


d_ 
dt 


(D.l) 


(D.2) 


(D.3) 


-  4  ff     ,A{t-h-T) 


Ah-, 


!l_f,e^^'-"-'^B^iu(z)dz  =  Aj;_i^e^^'-''-'>B,iuiT)dT  +  e-^^'B^uiz)  -  B^uit  -  h) 


Substituting  (D.4)  into  (D.3),  having 

z{t)  =  Az{t)  + 


B  +  e-'^'^B^ 


u{t)  +  f{x{t),t) 
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APPENDIX  E 


If  Ag'^H"  '\  B  G  9i"  '"  is  a  controllable  pair,  the  a  nonsingular  matrix  T  g  SK 


nxn 


can  be  found  such  that 


TAT 


Ai 

^u 

TB  = 

'  0' 

A2] 

A22_ 

A- 

then  the  system   Aj  j  ,  A12 


is  also  completely  controllable. 

Proof:  First  of  all,  we  note  that  controllability  (and  observability)  are  preserved 
under  similarity  transformation  (Kailath,  1980). 
Hence,  from  PBH  test,  given 


rank   si  -  A  ,B  \  =  n      \/s 


it  implies  that 


rank 


sI-TAT~\TB 


=  n      \/s 


(E.l) 


(E.2) 


That  is  no  left  eigenvector  of  TAT     is  orthogonal  to  TB  .  This  implies  that  there  can  not 
be  a  non  zero  row  vector  p  e  91"^   such  that 


si -TAT 


-1 


^0     V5 


(E.3) 


and 


pB^o    ys 


(E.4) 


We  partition  p  as 


P  =  [Pn-m      Pm]  (E-5) 

Using  (E.3)-(E.5)  in  conjunction  with  (5.11)  yields  an  equivalent  condition  for 
controllability. 
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The  system  [a,B]  is  controllable,  if  and  only  if  there  exists  no  non  zero  row 
vectors  p„_,„  e  9t""'" ,  /?,„  e  9t'"  such  that 


yPn-m      Pin  J 


5/  -  A, ,        -A,  o         0 


Ml 


12 


-A21  ^/-A22       i? 


;7„,fi=0      V5 


=  0      \/  s 


(E.6) 


(E.7) 


Note  that  since  B  is  nonsingular  square  matrix  the  only  solution  for  (E.7)  is  given  by 
jP„j  =  0  Therefore,  there  exists  no  row  vector  jP„_,„  e  9t^""'"^^  such  that 


/7„_,„(^/-A,i)  =  0      \f  s 


(E.8) 


i^n-wAz   =^       ^^ 


Of  course,  (E.8)  and  (E.9)  imply  that 


rank 


5/ -All  'A2 


=  n-m      Vs 


(E.9) 


(E.IO) 


Hence  the  dynamical  subsystem  An  ,  A12   is  completely  controllable. 


APPENDIX  F 


Consider  the  following  scalar  differential  equation 


/2/ 


x{t)<-ax{t)-bx'\t) 


with  known  x(0)  and  subject  to  the  constraints 
(1)     the  constants  a>0,  b>0 


(2)      x{t)  >  0  and  Jx{t)  >  0 


40^exp|  Y^ 


-a 


x{0)+- 


b 

a 


Proof: 


By  effecting  a  transformation 


x{t)  =  [y{t)f 


(F.l)  becomes 


2y{t)^<-ayHt)-by{t) 
dt 


Since  y{t)  >  0  from  the  constraints 


(F.l) 


(F.2) 


^<-f><"H-f 


(F.3) 


(F.3)  is  of  the  form 


dy 

—+py^q 

dt 


Hence, 


(F.4) 


^'(Oexpl  ^t\<  if Y  Rpf f  ^  \dt  +  C 
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4-&B'- 


<^exp(frKC 


yW^  — +Cexp[^-^? 


XO^expl^r 


-a 


J(0)  +  - 


a 


■^(O^expl  Y? 


x(0)+- 


(F.5) 
(F.6) 
(F.7) 

(F.8) 


The  transformation  from  (F.2)  to  (F.3)  is  valid  for  the  range  0  <  r  <  —  In 

a 


aJx{0) 


+  1 


For  all  practical  purposes,  the  upper  bound  can  be  taken  as  the  time  taken  for  the  system  to 

attain  zero. 

To  prove  the  assertion  in  the  corollary,  consider  the  scalar  differential  equation 

-|0.5(l+a) 


x{t)  <  -k[x{t)] 

with  known  x{0)  under  the  constraints: 
(1)     the  constant  k  >  0 


x(t)  >  0 


(2)     x{t)  >  0  and  Jx{t)  >  0 


(F.9) 


Integrating  (F.9),  it  is  straightforward  to  show 

0.5x(l-«) 


<-kt  +  c 


Substituting  for  the  constant  yields 

\l-a) 


x{t)< 


{I -a) 


[m] 


0.5(l-a) 


■kt 


2 


(F.IO) 


The  time  to  attain  x{t)  =  0  is  given  by 
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t< 

"     k{\-a) 


i/(0>(0) 


0.5(l~a) 


(F.ll) 


and  since 


x{t)  =  V{t)=^-s\Q)s{Q)  =  h^s{t)f 


m  ^  V2 


(1-a) 


(1-a) 


2 
1^ 


Q<t<t,.      (F.12) 


APPENDIX  G 
Consider  the  following  scalar  differential  equation 

x{t)  <  -ax{t)  +  [bexp{-ct/2)  +  e  exp[-d2t/2)  -  d)x^'^{t) 

with  the  initial  condition  x{Q)  and  the  following  constraints: 
(1)    The  constants  {a,b,c,d,e)  >  0  and  ai^c  and  a^  82 


(2)    x{t)  >  0  and  ^x{t)  >  0 

Then, 

v^< 

a_ 

f    a^         b 
K    2  J     [a-c) 

'       (    c\          ( 
exp  — t   -exp  - 
V    2  y          V 

a  V 
-  —  t 

2 ;. 

''{a -82) 

€ 

;xp        2  ^ 

( 
-exp  — 

V 

a  V 
-t 

2  ^ 

d 
a 

Proof:  By  effecting  a  transformation 


m-[y{t)\ 


(G.l)  becomes 


M) 


2y{t)^^^^  <  -ay^{t)  +  (^exp(-cr/2)  +  e  exp(-(52r/2)  -  d)y{t) 
Since  y{t)  >  0  from  constraints,  one  can  divide  both  sides  by  y{t)  to  get 
^^  <  (--\it)  +  [b/2exp{-ct/2)  +  e/2  exp{-d2t/2)  -  d/2) 


(G.3)  is  of  the  form 


dt 


+  py<q 


Hence,  the  solution  of  (G.3)  is  given  by 


(G.l) 


(G.2) 


(G.3) 
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3;(0exp|  —    <  j  bjltxA  —  I  +  e/2  exp 


^-5,t-^ 


-d/2 


^w'l^Wf '" 


^     e     ■^ 


a -do 


exp 


V    2    J 


expi  —  Wr  +  c 


d 

-  +  cexpl 

a  \  2  J 


(G.4) 


solving  for  the  constant  c,  yields 


b    \      f-ct)     d 
a- c J       \  2  J     a 


y{Q)  +  -- 


a      a  —  c     a  — 5' 


exp 


2,7 


-at 


It  is  instructive  to  rewrite  (G.5)  as 


,«.i„o).^ur^].f-i- 


ct  \  I  -at 

expl—   -exp  — 


+ 


^     .     ^ 


exp 


^-5,r^ 


-expl 


-a  A 


(G.5) 


(G.6) 


In  (G.6)  the  first  three  terms  are  positive,  while  the  third  term  is  a  negative  constant. 
Furthermore,  the  first  term  continuously  decreases  from  it's  initial  value,  while  the  second 
and  terms  always  have  a  maxima.  As  a  parallel  we  remark  that  the  second  and  third  terms 
in  (G.6)  look  like  the  response  of  an  under-damped  second  order  system  to  unit  impulse 
input. 


Since  y{t)  =  ■^x{t) ,  the  result  follows 
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